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IC  DESIGN 


1  Introduction 

The  design  of  high-performance  integrated  circuits  requires 
an  accurate*  inexpensive  way  in  which  to  assess  circuit  perfor¬ 
mance*  Discrete  "breadboard"  techniques  are  inadequate  because 
they  do  not  accurately  model  the  effects  of  parasitics*  thermal 
coupling*  and  device  matching*  which  are  critical  to  integrated- 
circuit  performance.  Prototype  circuit  construction  must  also  be 
ruled  out  because  of  the  high  cost  -  typically*  six  man-months 
and  $25000. 

As  a  result  of  these  considerations*  the  semiconductor  in¬ 
dustry  has  turned  to  the  integrated-circuit  (IC)  simulator  to 
evaluate  the  performance  of  circuits  during  the  design  stage. 
The  SPICE  program  was  developed  to  meet  the  needs  of  the  IC  group 
at  Berkeley.  The  analyses  which  it  performs  include  dc  operating 
point*  nonlinear  transient*  small-signal  frequency-domain*  noise* 
harmonic  distortion*  and  dc  sensitivities. 

This  report  describes  the  internal  design  of  the  SPICE2  pro¬ 
gram  (the  fundamental  theory  and  algorithms  are  described  in 
(D).  Chapter  2  describes  the  dynamic  memory  allocation  tech¬ 
nique  which  is  used  throughout  the  program.  Chapters  3-12  detail 
the  program  structure.  Finally*  the  Appendices  (Chapter  13)  con¬ 
tain  a  listing  of  the  User's  Guide*  specifications  for  all  the 
internal  data  structures*  descriptions  of  COI4MON  variables*  and  a 
program  listing. 
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2  Data  Management 

2.1  Description  of  Data 

As  execution  proceeds,  SPICE  requires  storage  for  many  dif¬ 
ferent  data  blocks*  Memory  must  be  allocated  to  hold  the  circuit 
description,  the  linearized  circuit  equations,  operating-point 
information,  charge-storage  information  for  transient  (large- 
signal)  analysis,  and  output  variables.  In  general,  both  the 
number  and  the  size  of  these  blocks  is  a  function  of  the  circuit 
to  be  analyzed. 

2.2  Previous  Method 

SPICEl  solved  this  storage  problem  by  allocating  a  fixed 
amount  of  storage  in  FORTRAN  COMMON-blocks  for  each  of  a  fixed 
number  of  data  blocks;  for  example,  space  for  400  nodes  and  25 
device  models  was  reserved.  This  approach  had  the  advantages  of 
simplicity  and  efficiency,  but  it  suffered  from  one  drawback: 
storage  reserved  for  one  type  of  data  could  not  be  used  for 
another  data  type*  Thus,  a  75-node  circuit  whose  description  re¬ 
quired  26  device  models  would  "overflow  memory",  even  though  the 
total  memory  available  was  sufficient. 

2.3  Design  Considerations 

There  were  four  principal  considerations  in  the  design  of 
the  memory  allocation  package  in  SP1CE2.  First  of  all,  the  im¬ 
plementation  language  had  to  be  FORTRAN  -  no  other  language  suit¬ 
able  for  scientific  computation  was  (or  is)  as  widely  available 
on  as  many  different  computers.  This  consideration  required  that 
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the  memory  to  be  managed  had  to  consist  of  a  single  FORTRAN  ar¬ 
ray;  this  array  is  called  "VALUE"  in  SPICE2. 

Secondly*  both  the  number  and  the  size  of  managed  data 
blocks  had  to  be  dynamically  variable*  This  consideration  elim¬ 
inated  the  following  type  of  "dynamic"  memory  allocation  (assum¬ 
ing  some  fixed  number  of  blocks): 

DIMENSION  A(IOOO) 

•  •  • 

•••  determine  sizes  of  blocks*  and  store 

•••  in  variaoles  I*  J*  and  K 

•  •  • 

CALL  D0IT(A(1 ) *  I *A (l+I ) * J *A (1  +  I+J ) * K ) 

•  •  • 

END 

SUBROUTINE  DOIT ( A  *  I  * B * J , C * K ) 

DIMENSION  A('l )  *B  (J  )  *C  (K) 

•  •  • 

•••  all  real  computation  performed  nere 

•  •  • 

RETURN 

END 

The  third  consideration  was  that  more  than  two  blocks  could 
be  varying  in  size  at  the  same  time*  At  first*  the  pattern  of 
memory  usage  in  SPICE2  permitted  memory  to  be  managed  as  a  last- 
in-first-out  (LIFO)  stack*  in  which  the  (only)  block  whose  size 
was  varying  at  any  given  moment  could  be  positioned  at  the  top 
(end)  of  allocated  memory.  Subsequent  enhancements  to  SPICE2  el¬ 
iminated  this  possibility. 

The  fourth  consideration  was  the  dynamic  region-size  capa¬ 
bility  which  some  computer  operating  systems  support  -  for  exam¬ 
ple*  the  SCOPE  (2)  operating  system  on  Control  Data  Corporation 
6000-series  computers.  Such  a  capability  offers  potential  reduc- 
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tions  in  both  response  time  and  analysis  cost  sufficient  to  jus- 
'tify  its  inclusion  in  the  memory  manager. 

2.4  Management  System  Description 

2.4.1  Definitions 

The  memory  manager  in  SPICE  controls  "tables"  using  "table 
pointers".  A  "table"  is  a  contiguous  block  of  memory;  a  "table 
pointer"  is  a  variable  which  serves  both  to  identify  a  block  (to 
serve  as  the  block  "name")  and  to  indicate  the  block  origin  in 
memory.  As  an  example>  consider  a  table  of  100  words  with  table 
pointer  "LTAB".  In  SPICE2/  the  global  allocation  array  is  named 
"VALUE";  the  contents  of  table  "LTAB"  would  oe  accessed  as 
VALUE (LTAB+1 )  through  VALUE { LTAB+100 ) ,  regardless  of  where  in 
VALUE  the  table  happened  to  be  at  any  given  moment. 

2.4.2  Memory  Manager  Calling  Sequences 

There  are  8  entry  points  in  the  memory  management  package: 
SETMEM,  GETMEM,  RELMEM,  EXTMEM,  CLRMEM,  SIZMEM,  PTRMEM ,  and 
CRUNCH.  The  rest  of  the  SPICE2  program  uses  only  these  entry 
points  for  all  memory  allocation  and  manipulation.  A  brief 
description  of  each  entry  follows: 


s,n.£„r,Y.(axqy:n]ga.tg) 

SETHEH(IFAMWA) 


GETMEM(P,S) 


RELMEH ( P , S ) 


description 

initializes  the  memory  manager;  the 
dynamically-managed  memory  starts  at 
(absolute)  address  IFAMWA 

makes  a  new  managed  table.  P  is  the 
table  pointer  for  this  new  table;  the 
size  of  the  new  table  is  S  words. 

reduces  by  S  words  the  size  of  the 
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EXTMEM ( P , S ) 


SIZMEW(P,S) 


CLRMEH ( P ) 

PTRMEM(P1 ,P2) 


CRUNCH 


table  pointed  to  oy  table  pointer  P* 

extends  by  S  words  the  size  of  the 
table  pointed  to  by  table  pointer  ?• 

sets  S  to  the  size  (in  words)  of  the 
table  pointed  to  by  table  pointer  P« 

destroys  the  table  pointed  to  by 
table  pointer  P> 

changes  the  table  pointer  for  the 
table  pointed  to  by  table  pointer  Pi 
to  be  P2. 

forces  a  compaction  of  the  dynamically- 
managed  memory 


2.4.3  Memory  Allocation  Algorithms 

The  memory  manager  maintains  a  "table  entry"  table  at  the 
high-address  end  of  managed  memory.  This  table  contains  a  four- 
word  entry  for  each  allocated  table  of  the  form: 


wacj  contents 

1  table  origin  (array  subscript) 

2  allocated  table  size  (in  words) 

3  requested  block  size  (in  words) 

4  address  of  table  pointer 


The  memory  manager  uses  the  redundant  information  in  the 
"table  entry"  table  for  error-checking  (verifying#  for  example, 
that  the  value  of  a  table  pointer  is  the  same  as  the  appropriate 
entry  in  the  "table  entry"  table).  Note  that  any  allocated  table 
may  in  fact  be  larger  than  required;  that  fact  would  be  indicated 
by  an  "allocated  size"  which  was  greater  than  the  "requested 
size".  The  memory  manager  assumes  that  arguments  are  passed  "by 
address"  since  it  uses  the  address  of  the  table  pointer  to  update 
its  value  should  the  table  origin  change. 
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There  are  two  basic  algorithms  used  in  the  memory  manager 
one  to  allocate  a  new  tables  the  other  to  extend  an  existing 
table-  The  allocation  algorithm  can  be  described  by  the  follow¬ 
ing: 


call  CRUNCH; 

if  (memory  insufficient)  get  more  memory; 

add  entry  at  end  of  "table  entry"  table  for  new  table; 

store  table  origin  into  table  pointer  variable; 


The  variable  MEMAVL  (memory  available)  in  the  memory  manager 
always  contains  the  number  of  words  of  memory  which  are  free  (al¬ 
located  to  some  Dlockf  but  more  tnan  requested).  Therefore  the 
"insufficient  memory"  test  is  a  simple  check  of  MEMAVL. 

The  algorithm  to  extend  a  table  is  the  following: 


if  (required  size  >  allocated  size) 

{  if  (insufficient  memory)  get  more  memory; 
call  COMPRS  to  move  tables  below  the  table 
being  extended  towards  the  origin  of 
the  dynamically-managed  memory; 
call  COMPRS  to  move  tables  above  the  table 
being  extended  towards  the  end  of  the 
dynamically-managed  memory; 

) 

store  new  required  size  in  "table  entry"  table; 


The  taole  extension  algorithm  takes  advantage  of  the  fact 
that  when  a  given  table  is  extendedi  roost  of  the  time  the  next 
call  on  the  memory  manager  will  be  to  extend  the  same  table 
again.  Note  that  neither  algorithm  changes  the  origin  of  the 
first  table  allocated  (assuming  that  it  is  not  subsequently  des¬ 
troyed  ) . 
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3  Program  Structure 

The  SPICE2  program  consists  of  15000  FORTRAN  and  COMPASS 
statements  divided  into  7  major  overlays  together  with  an  overlay 
root. 

A  composite  memory  map  for  SPICE  execution  is  shown  below. 
At  any  one  moment#  only  the  root  and  one  of  the  indicated  first- 
level  overlays  is  resident  in  memory.  The  dynamically-managed 
memory  begins  after  the  longest  overlay  (as  indicated  below). 


SPICE  overlay  root 


•  • 

•  • 

READIN  :  ERRCHK  :  SETUP 

•  • 

•  • 

DCTRAN 

•  • 

•  • 

DCOP  :  ACAN  :  OVTPVT 

•  • 

•  • 

•  • 

•  • 

DYNAMICALLY -MANAGED  MEMORY  (COtdi^ON-block  /BLANK/) 

<-IFAMWA 

<-ICORE 


The  variables  IFAMWA  and  ICORE  are  explained  in  Sections  2.4.2 
and  13.4.1. 


Each  section  of  the  program  is  discussed  separately  in  the 
following  chapters  of  this  report. 
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4  Overlay  Root 


The  SPICE  overlay  root  drives  the  rest  of  the  program#  cal¬ 
ling  the  first-level  overlays  to  process  the  circuit  description, 
analyze  the  circuit#  and  output  the  simulation  results.  The  root 
consists  of  the  main  program  together  with  routines  TITLE#  SET- 
MEM#  COMPRS#  MEMPTR,  DMPMEM#  TMPUPD#  HAGPHS#  XOR#  OUTNAM#  ALFWUM# 
FIND,  MEMORY#  DCDCHP#  DCSOL#  MOVE#  COPY#  and  ZERO.  The  main  loop 
is  described  in  the  following  figure: 


initialize; 

while  (.not.  end-of -input ) 

{  call  READIN;  call  ERRCHX;  call  SETUP; 
repeat 

{  TEMP  «  next  analysis  temperature;  ■  call  TMPUPD; 
if  (dc  transfer  curves  requested) 

{  MODE  «  1;  MODEDC  «  3; 

call  DCTRAN;  call  OVTPVT?  > 

if  (dc  operating  point  or  ac  analysis  requested) 

{  MODE  »  1;  MODEDC  «  1; 

call  DCTRAN;  call  DCOP;  } 
if  (ac  analysis  requested) 

{  MODE  -  3;  call  ACAN;  call  OVTPVT;  } 
if  (transient  analysis  requested) 

{  MODE  »  1;  MODEDC  «  2; 
call  DCTRAN;  call  DCOP; 

MODE  »  2; 

call  DCTRAN;  call  OVTPVT;  > 

) 

until  (no  more  analysis  temperatures); 

> 

stop 


The  READIN  overlay  reads  input  until  either  a  line  contain¬ 
ing  '.END'  is  read  or  end-of-file  is  found*  As  the  overlay  exe¬ 
cutes#  it  builds  the  linXed-lists  which  SPICE  uses  to  store  the 
circuit  description#  and  it  sets  certain  flags  in  COMI-ION  to  indi¬ 
cate  what  analyses  were  requested.  In  the  ERRCHK  overlay,  SPICE 
expands  the  suocircuit  calls#  prints  a  circuit  summary#  and 
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checks  for  cornmon  input  errors.  The  SETUP  overlay  constructs  the 
pointers  which  are  used  to  manipulate  the  sparse  equation  coeffi¬ 
cient  matrix  in  both  the  DCTRAN  and  ACAN  analysis  overlays.  This 
overlay  also  generates  the  machine  code  which  is  executed  to 
solve  the  circuit  equations. 

The  main  analysis  loop  consists  of  the  "repeat*  loop  in  the 
above  figure.  For  each  analysis  temperaturei  it  checks  succes¬ 
sively  for  a  dc  transfer  curve  analysis,  an  dc  operating-point 
analysis,  an  ac  (small-signal)  analysis,  and  a  transient  (large- 
signal)  analysis.  The  dc  operating-point  analysis  is  always  per¬ 
formed  if  an  ac  small-signal  analysis  is  requested  because  SPICE 
requires  the  small-signal  equivalent  models  for  the  nonlinear 
elements  in  order  to  perform  the  ac  analysis. 

The  DCTRAN  overlay  performs  the  dc  transfer  curve,  dc 
operating-point.  transient  initial-condition,  and  transient  ana¬ 
lyses.  Exactly  which  of  these  analyses  it  performs  is  determined 
by  the  values  of  the  COMMON  variables  MODE  and  MODEDC.  The  ACAN 
overlay  performs  the  ac  analysis.  Finally.  the  OVTPVT  overlay 
generates  the  tabular  output  listings  and  line-printer  plots  from 
the  results  of  the  multi-point  analyses. 

The  overlay  root  contains  several  utility  subroutines  which 
are  used  throughout  the  rest  of  the  program.  A  brief  description 
of  those  routines  follows  (A  and  B  are  array  names): 
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Dame(arauments) 

MOVE(A,I,B,J,N) 

COPY4 (A,B,N) 


COPY 6 (A,B,N) 
COPY16 (A,B ,N ) 
ZER04 (A,N ) 
ZEROS (A,N ) 
ZEROie (A.N) 


description 

moves  N  characters  from  B  to  A,  starting 
at  the  Jth  and  Ith  character  positions 
respectively;  eight  characters/word 
is  assumed 

moves  N  Integers  from  A  to  B;  if  A  and  3 
overlap  in  memory  -  for  example#  if 
B  is  really  A(10}>  to  move  data  within 
a  single  array  -  C0PY4  copies  either 
A(l)-to-A(N)  or  A(N)-to-A(l)  to  prevent 
the  transfer  from  destroying  the  data 

C0PY4()#  but  for  real  variables 

C0PY4 (  )  ,  out  for  complex  variaoles 

zeroes  A(l)  -  A{N)  integers 

ZER04 ( ) I  but  for  real  variaoles 

ZER04 ( } ,  but  for  complex  variables 
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5  Readin  Overlay 

The  READIN  overlay  perforins  the  initial  input-processing  for 
5PICE>  This  overlay  consists  of  the  routines  READIN*  KEYSRC* 
EXTNAM,  RUNCON,  OUTDEF,  CARD,  GETLIN.  and  NXTCHR.  It  builds  the 
linked-lists  which  describe  the  circuit  to  be  analyzed,  and  it 
sets  certain  flags  in  C0[<1M0N  to  indicate  which  analyses  have  been 
requested.  A  flowchart  for  the  overlay  follows: 


if  (end-of-file  found)  exit; 
if  (end-of-file  found)  exit; 


read  title; 
initialize; 
call  CARD; 
repeat 

<  switch  (type  of  input  line) 
{  element  description: 

“.ENDS": 
".MODEL": 
“.END": 
“.SU3CKT": 
other  “.“  line: 


process  in  READIN; 
process  in  READIN; 
process  in  READIN; 
exit; 

process  in  READIN; 
call  RUNCON; 


} 

call  CARD; 

>  until  (end-of-file  found); 
exit ; 


5.1  Scanning 

5.1.1  Readin 

The  READIN  subroutine  processes  element  and  device-model  de¬ 
finitions.  The  routine  calls  CARD  to  scan  the  physical  input 
lines.  Subroutine  FIND  is  called  to  preset  storage  for  each  in¬ 
put  element.  Analysis  and  output  directives  are  processed  oy 
RUNCON. 
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5«1.2  Card 

Subroutine  CARD  is  called  by  READIN  to  scan  the  next  logical 
input  line.  Since  physical  lines  may  be  logically  continued  {by 
putting  a  in  column  one  of  the  next  physical  line).  CARD  as¬ 
sumes  on  entry  that  the  next  line  has  already  been  read  from  in¬ 
put.  and  always  leaves  the  next  input  line  (unscanned  except  for 
the  initial  '  +  '  check)  in  the  array  AFIELD  (in  COMi'lON-block 
LINE) . 

Each  logical  input  line  is  parsed  into  a  set  of  "numbers" 
and  "names".  A  field  is  assumed  to  be  a  "name"  if  it  begins  with 
either  an  alphabetic  character  or  a  period  (".")  followed  by  an 
alphabetic  character.  A  field  is  a  "number"  if  it  begins  with 
either  a  digit  (0-9)  or  a  period  (".")  followed  ny  a  digit. 
Fields  are  delimited  by  any  character  in  the  set  {"  ".  ".". 

"(".  ")">;  in  the  present  version,  all  delimiters  are  equivalent 
except  with  processing  output  variable  definitions  (see 
subroutine  OUTDEF).  The  results  of  the  scan  are  stored  in  the 
tables  IFIELD.  ICODE.  IDELIM,  and  ICOLUM,  all  of  which  track  with 
each  field  of  the  line  (the  first  word  of  each  of  the  tables 
refers  to  the  first  field,  the  second  word  for  the  second  field, 
etc.).  The  entry  in  ICODE  is  +1,  0,  or  -1  to  indicate  "name", 
"number",  and  end-of-line.  The  IFIELD  entries  are  either  a 
floating-point  number  or  a  name  stored  in  8H  format  (left- 
justified.  blank-filled  to  3  characters).  The  IDELIM  table  con¬ 
tains  the  field  delimiter  character.  Finally,  the  ICOLUM  tacle 
contains  the  column  number  at  which  the  field  oegan  in  the  input 
line. 
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5.1.3  Runcon 

The  RUNCON  subroutine  processes  run  control  lines  -  for  ex¬ 
ample/  a  ".TRAN"  card  indicating  that  a  transient  analysis  is  to 
be  done.  It  sets  the  pertinent  flags  in  COMI’ION  to  indicate  which 
analyses  SPICE  is  to  perform  (see  Section  13.4  for  a  description 
of  the  COI^J'ION  flags).  Some  error-checlcing  is  performed  -  for  ex¬ 
ample/  that  required  analysis  parameters  have  not  been  omitted. 

5.1.4  Find 

The  FIND  subroutine  is  used  to  locate  a  particular  list  ele¬ 
ment  or  to  allocate  space  if  the  element  cannot  be  found.  It  has 
four  parameters/  described  in  the  following: 


Parameter  description 


ANAME 

ID 

LOC 

IFORCE 


name  of  list  element 

number  of  list  to  search 

set  by  FIND  to  location  of  element 

1  *>  element  must  not  already  exist 


If  NSBCKT  is  nonzero/  a  subcircuit  definition  is  being  read. 
In  that  event/  FIND  searches  the  current  suocircuit  definition 
list  (pointed  to  by  the  subcircuit  whose  definition  is  pointed  to 
by  NODPLC ( IS3CKT+NS3CKT ) )  instead  of  the  nominal  element  list 
(pointed  to  by  LOCATEdD)). 


5.2  Storage 
5.2.1  Elements 

All  input  data  to  SPICE2  is  maintained  in  a  series  of  linked 
lists  stored  in  table  IELi4NT.  The  LOCATE  array  (in  COMHON-block 
CIRDAT)  contains  pointers  (subscripts  to  the  NODPLC  array)  to  the 
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beginning  of  each  list.  By  convention,  the  integer-  and  real¬ 
valued  data  are  accessed  using  separate  suoscripts  for  the 
(equivalenced }  arrays  NODPLC  and  VALUE*  Each  list  entry  contains 
a  "next-pointer"  which  is  the  NODPLC-subscript  of  the  next  list 
element  (zero  if  none).  Each  entry  also  contains  a  VAI,UE- 
subscript  (usually  called  LOCV),  which  points  to  the  first  word 
of  real-valued  storage  for  the  list  element. 

A  detailed  description  of  every  linked-list  element  struc¬ 
ture  is  contained  in  Section  13*2. 

5*2*2  Analyses  and  Outputs 

SPICE2  performs  three  principal  analyses:  dc,  ac,  and  tran¬ 
sient*  For  each  of  these  analyses,  a  separate  COMMON-block  is 
used  to  store  pertinent  information*  For  example,  COMHON-block 
/TRAN/  contains  variables  TSTEP,  T3TOP,  and  TSTART,  which  are  set 
to  the  print-step,  final  simulation  time,  and  starting  simulation 
time,  respectively.  A  detailed  description  of  the  COl-U'lON- 
variables  is  contained  in  Section  13*4* 

Analysis  output  variables  are  stored  in  table  LOUTPT  as 
analysis  proceeds*  After  each  sweep  point  (of  some  source  (dc), 
frequence  (ac),  or  time  (transient))  is  solved,  tne  LOUTPT  table 
is  extended  and  the  values  of  the  output  variables  at  that  point 
stored  away*  The  variable  ISEU  (see  Section  13*2*24)  contains 
the  offset  for  the  particular  output  variable  within  each  exten¬ 
sion  of  the  LOUTPT  table*  Since  the  results  of  a  given  analysis 
are  printed  before  another  analysis  is  begun,  only  one  output 
table  is  necessary*  The  output  variaole  definitions  themselves 
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are  also  stored  in  linked-list  form  in  table  lELMWT  (with  ID=41 
through  45). 

5.3  Subcircuit  Structures 

During  READIN#  each  subcircuit  definition  is  stored  as  a 
linked  list  of  elements  and  device  models.  The  "subckt  info" 
field  in  each  element  of  a  subcircuit  definition  list  contains 
the  element  ID.  Subcircuits  are  completely  local  entities  in 
SPICE2;  elements  or  device  models  introduced  within  a  subcircuit 
definition  are  not  known  outside  that  definition.  As  a  result, 
if  the  line 

01  1  2  3  MODA 

is  read  as  part  of  a  subcircuit  definition.  SPICE  does  not  in 
general  know  where  device  model  MODA  is  defined  (it  may  be  de¬ 
fined  later  within  this  subcircuit  definition,  or  it  may  be  de¬ 
fined  in  some  enclosing  subcircuit  definition,  or  it  may  be  de¬ 
fined  in  the  nominal  circuit  definition).  To  solve  this  problem 

\ 

correctly.  an  "unsatisfied  name"  table.  lUNSAT.  was  added;  this 
table  contains  all  the  names  (such  as  "MODA").  During  READIN, 
the  various  element  definitions  contain  pointers  into  the  lUNSAT 
table.  Proper  resolution  of  name  references  is  made  in 
subroutine  LNKREF  (in  the  ERRCUK  overlay)  and  is  discussed  in 
Section  6.4. 
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5 >4  Forward  References 

For  efficiency  reasons,  the  input  file  is  read  only  once  by 
SPICE*  As  a  result,  a  reference  to  some  element  may  be  read  be¬ 
fore  the  element  definition  itself*  An  example  of  this  situation 
is  shown  below: 

•  •  • 

*DC  VIN  0  5  *5 

*  *  *  (other  input) 

VIN  4  0  PULSE (0  5) 

•  •  • 

When  the  ".DC"  line  is  read.  FIND  is  called  with  the  "force- 
action"  parameter  IFORCE  set  to  0*  The  effect  of  this  call  is  to 
add  an  "undefined  but  referenced"  element  with  name  "VIN"  to  the 
circuit.  For  such  a  case,  subroutine  FIND  sets  the  (LOC  +  2)th 
word  of  the  lis'-  entry  to  "*UN"  (for  "undefined").  (Properly 
speaking.  the  lUNSAT  table  should  also  be  used  for  this  purpose: 
at  present,  elements  within  subcircuits  are  not  treated  in  the 
same  manner  as  elements  in  the  nominal  circuit  description. ) 
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6  Errchk  Overlay 

The  ERRCHK  overlay  finishes  the  processing  of  input  data  and 
perforins  miscellaneous  error  checking*  The  overlay  consists  of 
the  routines  ERRCHK,  SHLSRT,  PUTNOD,  GETNOD,  SUBCKT,  FWDNAM, 
NEWNOD,  ADDELT,  CPYTAB,  LNKREF,  SUBNAM,  ELPRNT,  MODCHK,  and 
TOPCHX*  A  flowchart  describing  the  actions  of  this  overlay  fol¬ 
lows: 


check  forward  references; 

generate  compact  node  set; 

call  SUBCKT  to  expand  suocircuit  calls; 

call  LNKREF  to  link-up  external  names; 

set  source  function  defaults  and  limits; 

call  ELPRNT  to  print  circuit  element  summary; 

call  MODCHK  to  print  and  preprocess  device  models; 

call  TOPCHK  to  check  topology  and  print  node  table; 

invert  resistance  values; 

change  K  to  M  for  mutual  inductance; 

generate  breakpoint  table; 

check  analysis  limits; 

sequence  output  variables; 


6*1  Forward  Reference  Check 

The  forward-reference  problem  is  explained  in  Section  5*4* 
In  ERRCHK,  all  that  need  be  done  is  to  sequence  through  all  the 
element  lists  checking  whether  the  (LOC  +  2)th  word  of  any  entry 
has  the  value  ".UN".  Since  all  element  lists  store  either  a  node 
number  or  other  small  integer  in  the  (LOC  +  2)th  position,  the 
presence  of  that  value  indicates  that  the  element  was  never  actu¬ 
ally  read  in,  but  only  referenced. 
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6 >2  Compact  Element  Node  Set 

After  the  READIN  overlay  has  completed/  the  circuit  element 
list  entries  contain  the  input  node  numbers*  Since  the  analysis 
portion  of  SPICE2  requires  a  compact/  consecutive  set  of  node 
numbers/  ERRCHK  renumbers  all  element  nodes  and  stores  the  new 
node  numbers  in  place  of  the  user  node  numbers*  A  translation 
table/  JUNODE/  is  constructed  and  used  when  SPICE  has  to  print 
the  original  (user)  node  number;  if  NODPLC (LOC+2 )  contains  the 
renumbered  node  number/  NODPLC (JUNODZ+NODPLC (LOC+2 ) )  contains  the 
input  node  number* 

6*3  Subcircuit  Expansion 

At  the  end  of  execution  of  the  READIN  overlay/  all  subcir¬ 
cuit  definitions  have  been  read  in  and  stored  on  the  ID«20  linked 
list/  and  all  subcircuit  calls  have  been  read  in  and  stored  on 
the  ID*19  list*  Furthermore/  before  subcircuit  expansion  is  per¬ 
formed/  the  initial  compact  node  set  has  been  constructed/  so 
that  all  user  node  numbers  (and  all  internal  node  numbers)  are 
known*  This  last  point  is  important  since  additional  (noncon¬ 
flicting)  "user"  node  numbers  may  have  to  be  generated  during 
subcircuit  expansion* 

All  elements  and  device  models/  contained  within  any  subcir¬ 
cuit  definition  are  treated  as  purely  local  to  that  definition  by 
SPICE2*  Consider  the  following  input: 
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XI  1  2  3  SUBl 
.SUBCKT  SJ31  17  15  23 
01  17  3  15  2i40D 
•  MODEL  QMOD37  NPN  (  ) 

02  23  17  6  OMOD37 
R123  3  6  IK 
XI  17  23  67  SU32 
•ENOS  SUBl 


•MODEL  OHOD37  PNP ( ) 
05  1  17  45  OMOD37 

•  •  • 


Mote  that  OMOD »  the  device  model  for  suocircuit  element  Ql » 
is  not  defined  within  the  sudcircuit^  SPICE2  will  search  enclos¬ 
ing  subcircuit  definitions  (if  any)/  and  then  the  nominal  circuit 
description/  for  that  device  model  definition^  Note  that 
although  OMOD37  is  defined  ooth  within  the  subcircuit  definition 
and  outside  that  definition/  the  transistor  02  within  the  suDCir- 
cuit  will  use  the  device  model  definition  contained  within  the 
suocircuit/  and  transistor  05  will  not^  (For  transistor  05/ 
SPICE2  will  try  to  find  a  device  model  0i'iOD37  which  is  defined  in 
an  outer/  enclosing  subcircuit  definition  (if  any)  or  within  the 
nominal  circuit  description^ ) 

The  actual  expansion  algorithm  used  in  subroutine  SU3CKT  is 
as  follows: 
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LOCX  *  LOCAT£(19);  (pointer  to  first  subcict  call) 
while  (LOCX  is  nonzero) 

{  call  FWDNAM  to  determine  subcircuit; 
check  for  recursive  definition; 

LOC  *  first  subcircuit  element; 
while  (LOC  is  nonzero) 

<  call  FIND  to  add  element  space; 
set  "subckt  info"  appropriately; 
call  ADOELT  to  add  element  data; 

LOC  *  pointer  to  next  suockt  element; 

> 

LOCX  »  pointer  to  next  X-element; 


Subcircuits  are  expanded  "top-down"  since  any  "X-element"  within 
a  subcircuit  definition  being  expanded  will  be  added  to  the  end 
of  the  subcircuit  call  list.  The  test  for  recursion  checks  that 
the  X-element  which  caused  a  particular  subcircuit  to  be  expanded 
was  not  itself  added  as  a  result  of  a  previous  expansion  of  that 
subcircuit.  The  test  is  simplified  by  the  fact  that  "subckt  in¬ 
fo"  for  the  "X-elemert"  causing  the  expansion  points  to  the  sub¬ 
circuit  definition  of  which  it  was  a  part  (if  any). 

The  call  to  subroutine  FIND  is  somewhat  modified  since  no 
naming  conflict  is  possible  at  this  point  in  execution.  There- 
forei  a  unique  "element  name"  is  constructed  from  the  number  of 
elements  of  the  type  being  added. 

The  FNDNAM  subroutine  determines  what  element  is  meant  by 
some  name  in  the  lUNSAT  table  (see  Section  5.3).  The  routine 
first  searches  the  subcircuit  expansion  list  (which  is  ootained 
from  the  "subckt  info"  word  of  each  list  entry)  in  reverse  order. 
If  no  match  is  found,  FNDNAM  then  searches  the  nominal  circuit 
description  list. 
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6-4  lU.'JSAT  Resolution 

After  suDcircuit  expansion.  "unsatisfied"  name  references 
must  oe  resolved  before  the  complete  circuit  description  is 
available  (see  Section  5.3).  This  resolution  is  performed  by 
subroutine  LNKREF.  which  sequences  through  the  linked  lists,  cal¬ 
ling  FNDNAM  to  resolve  the  reference (s)  for  each  circuit  element. 

6.5  Circuit  Summary 

Before  printing  the  circuit  summary,  the  ERRCHK  suoroutine 
first  makes  a  pass  through  the  independent  sources,  setting  de¬ 
fault  function  values.  Then  subroutine  ELPR.'JT  is  called.  which 
prints  out  a  summary  of  the  circuit  elements. 

To  print  out  the  device  model  summary,  suoroutine  MODCHK  is 
called.  '  It  first  assigns  default  model  parameters  from  the  de¬ 
fault  parameter  list  OEFVAL.  The  device  models  are  then  printed. 
After  printing  the  parameter  values.  MODCHK  performs  one-time 
model  parameter  computations,  such  as  replacing  resistances  by 
conouctances  (details  regarding  the  pre-processing  are  contained 
in  Sections  13.2.18  -  13.2.21).  Finally.  MODCHK  generates  addi¬ 
tional  node  nuraoers  for  those  nodes  adoed  to  accommodate  nonzero 
device  model  resistors. 

6.6  Node  Table  /  Topology  check 

Subroutine  TOPCHK  is  called  by  ERRCHK  to  construct  the  cir¬ 
cuit  node  table  and  to  check  that  three  topological  constraints 
are  satisfied  oy  the  circuit  to  be  analyzed:  no  voltage  source 
and/or  inductor  loop,  no  current  source  and/or  capacitor  cutset, 
and  no  node  with  less  than  two  branches  connected  exists. 
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The  node  table  is  constructed  using  taoles  lUR/  ITABLEf  ana 
ITABID.  The  entries  in  the  lUR  table  are  offsets  within  the  IT- 
ABLE  and  ITABID  tables*  The  entries  in  the  ITABLE  table  are 
pointers  to  circuit  elements;  entries  in  ITABID  are  the  element 
ID  for  the  corresponding  pointer  in  ITABLE*  After  the  tables 
have  been  constructed#  the  elements  connected  to  node  number  N 
are  pointed  to  by 

NOD?LC(ITABLE+NODPLC(IUR+N) ) 
through 

MODPLC ( ITABLE+NODPLC ( I UR+N+1 ) -1 ) 

6*7  Breakpoint  Table 

In  transient  analysis.  SPICE2  always  uses  a  program- 
calculated  timestepf  regardless  of  the  user-specified  print  in¬ 
terval*  However,  the  independent  source  waveforms  frequently 
have  sharp  transitions  which  could  cause  an  unnecessary  reduction 
in  the  timestep  in  order  to  find  the  exact  transition  time*  To 
overcome  this  problem.  ERRCHK  generates  a  "nreaxpoint"  table, 
LSBKPT,  which  contains  a  sorted  list  of  all  the  transition  points 
of  the  independent  sources*  During  transient  analysis,  whenever 
the  next  timepoint  is  sufficiently  close  to  one  of  the  break¬ 
points,  the  timestep  is  adjusted  so  that  the  program  lands  exact¬ 
ly  on  the  breakpoint* 
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7  Setup  Overlay 

The  SETUP  overlay  generates  the  pointer  arrays  necessary  .  to 
manipulate  the  sparse  equation  matrix  used  in  all  of  the  ana* 
lyses*  The  overlay  consists  of  routines  SETUP*  HATPTR*  RESERV* 
REORDR,  MATLOC,  INDEX,  CODGEN ,  and  MINS.  A  flowchart  of  this 
overlay  follows: 

NSTOP  =  number  of  equations; 

call  MATPTR  to  establish  matrix  structure; 

call  REORDR  to  reorder  equations  for  sparsity; 

call  MATLiOC  to  store  matrix  locations; 

if  (machine  code  requested) 

call  CODGEN  to  perform  code  generation; 

The  subroutines  MATPTR,  REORDR,  and  MATLOC  respectively  es¬ 
tablish  the  equation  matrix  structure,  reorder  the  equations  to 
maximize  sparsity,  and  store  the  matrix  pointers  for  each  circuit 
element.  The  CODGEN  subroutine  generates  machine  instructions  to 
solve  the  set  of  equations. 

7.1  Matrix  Structure 

The  equation  matrix  is  very  sparse  -  typically,  over  85%  of 
its  entries  are  zero.  The  matrix  structure  is  stored  internally 
in  two  different  forms  during  the  setup  process.  Initially,  the 
nonzero  matrix  entries  are  recorded  using  a  set  of  linJced  lists, 
one  for  each  row  of  the  matrix,  together  with  some  auxiliary 
tables.  The  i'th  entry  of  the  ISR  table  is  a  pointer  to  the  be¬ 
ginning  of  a  lin)ced  list,  each  of  whose  entries  records  the 
numoer  of  some  nonzero  column  in  the  i'th  row.  The  linked  list 
entries  themselves  are  stored  at  the  end  of  the  NUMOFF  table. 
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The  auxiliary  tables  used  are  NDIAGi  NUMOFF,  and  NMOFFC.  The 
NDIAG  table  records  whether  or  not  the  diagonal  matrix  elements 
are  nonzero,  the  first  N  words  of  the  NUMOFF  table  record  the 
number  of  nonzero  columns  in  each  row  of  the  matrix,  and  the 
NMOFFC  table  contains  the  number  of  nonzero  rows  u.n  each  column 
of  the  matrix. 

The  linked-list  structure  used  for. matrix  setup  has  the  ad¬ 
vantage  that  only  a  constant  (small)  overhead  is  required  to 
modify  the  matrix  structure.  However,  the  linked  list  method  is 
not  desirable  for  analysis  (given  that  SPICE2  does  not  use  pivot¬ 
ing)  since  such  a  structure  increases  both  memory  and  cpu  re¬ 
quirements..  Because  of  these  considerations,  after  the  final  ma¬ 
trix  structure  has  been  estaolished,  SPICE  generates  a  compact 
matrix  pointer  system. 

The  compact  pointer  system  divides  the  equation  matrix  into 
three  parts:  the  matrix  diagonal,  an  upper-,  and  a  lower- 
triangular  section-  The  diagonal  is  stored  separately  since  it 
is  not  sparse.  For  the  upper  triangle,  the  tables  lUR  and  lUC 
are  used.  Consider  the  i'th  row  of  the  equation  matrix;  let 
II  =  NODPLC (lUR  +  i  )  and  let  12  =  NODPLC ( I UR  +  i  +  1 ) -1 .  Then  the 
numbers  of  the  nonzero  columns  in  the  upper-triangular  part  of 
the  i'th  row  are  stored  in  NODPLC ( lUC+Il )  through  NODPLC (IUC+I2 ) . 
(If  I1>I2  then  there  are  no  such  columns-)  Tne  actual  matrix  en¬ 
tries  are  stored  in  taole  LVN  "in  parallel"  to  the  lUC  table-  In 
other  words,  having  found  the  offset  within  the  lUC  table  for 
some  (i,^)  matrix  element,  than  the  actual  location  within  the 
compressed  matrix  storage  of  the  element  value  is  immediately 
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known.  The  lower-triangula r  part  of  the  matrix  is  stored  using  a 
dual  set  of  matrix  pointers^  in  taoles  ILC  (the  ILR  index)  and 
ILR. 


7.2  Row-Swap 

SPICE2  uses  the  Modified  Nodal  formulation  [3] ,  [4]  to  con¬ 

struct  the  set  of  equations  describing  the  input  circuit. 
Although  that  set  of  equations  is  nonsingular  for  any  circuit 
satisfying  the  topological  constraints  mentioned  in  Section  6.6« 
problems  may  still  arise  when  no  provision  for  pivoting  is  incor¬ 
porated  in  the  program.  These  problems  arise  because  the  subset 
of  equations  added  to  solve  for  currents  flowing  in  "voltage- 
defined"  branches  may  be  singular  [1»  pp.  67-69]. 

Basically/  the  solution  to  this  singularity  problem  is  to 
swap  the  matrix  row  which  solves  for  a  branch  current  with  the 
corresponding  row  which  solves  for  the  voltage  at  the  positive 
node  of  that  branch.  However/  the  order  in  which  rows  are 
swapped  is  not  arbitrary  if  there  are  trees  of  voltage-defined 
elements.  The  algorithm  used  in  SPICE  to  decide  on  the  order  of 
row-swaps  is  the  following: 

Tables  already  constructed: 

NODEVS:  #  of  voltage-defined  elements  connected  to  node 

NDIAG:  nonzero  if  conductance  connected  to  node 

ISEj:  vector  of  pointers  to  all  voltage-defined  elements 

let  NUMVDE  *  f  of  voltage-defined  elements 
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for  (I  -  1  to  HUMVDE) 

{  FOUND  »  . . false. ; 

for  (J  »  I  to  NUMVDE) 

{  (examine  element  pointed  to  by  ISEQ(J)) 
for  (NODE  =  each  of  the  2  element  nodes) 

{  if  (NODE  *  ground)  continue; 

if  (NODEVS(NODE)  >  1)  continue; 
if  (NDIAG(NODS)  «  0)  go  to  swap; 

FOUND  «  .true.;  save  J; 

> 

} 

if  (FOUND  «  .false.  ) 

<  abort  -  voltage-defined  element  loop  exists  > 

swap: 

exchange  rows  for  element  indicated  oy  J ; 
delete  element  from  ISEQ  table; 

> 


7.3  Reorder 

As  previously  noted,  the  typical  equation  matrix  is  very 
sparse.  Howver,  the  order  in  which  unJcnowns  are  solved  can  have 
a  considerable  effect  on  sparsity  due  to  the  generation  and  pro¬ 
pagation  of  "fill-in"  terras  (zero  matrix  entries  which  become 
nonzero  during  the  LU  decomposition).  Various  algorithms  to  de¬ 
cide  on  .  the  order  of  variable  elimination  have  been  proposed 
(Barry  [5],  Markowitz  (6],  Nakhla  (7)).  SPICE2  uses  the  Mar¬ 
kowitz  algorithm;  the  number  of  matrix  terms  per  row  and  column, 
which  the  method  requires,  is  obtained  from  the  NUMOFF  and  NMOFFC 
tables,  respectively. 


7.4  Element  Matrix  Pointers 

After  the  final  matrix  structure  has  been  established, 
subroutine  MATLOC  is  called.  To  minimize  element  load  time,  the 
offsets  within  the  LVN  table  to  each  matrix  term  which  a  given 
circuit  element  increments  are  stored  with  each  circuit  element. 
In  this  way,  the  overhead  of  mapping  an  (i,^)  matrix  location  to 
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its  proper  offset  within  the  LVN  table  is  incurred  only  once. 
For  example,  a  resistor  between  nodes  W1  and  il2  causes  an  addi¬ 
tion  (or  suotraction)  to  matrix  terms  (i']l,i'jl)«  (NlfN2},  (N2,i']l), 
and  (N2,I']2).  The  locations  of  each  of  these  terms  is  stored  with 
each  resistor. 

7.5  Machine  Code  Generation  and  Execution 

Suoroutines  CODGEiJ  and  HIWS  are  used  to  generate  the  loop¬ 
less  machine  instructions  which  are  executed  to  perform  the  LU 
decomposition  and  forward/backward  suostitution  steps  which  solve 
the  linearized  circuit  equations.  The  code  generation  process  is 
simplified  oy  noting  that  two  forms  of  assignment  statements  can 
be  used  to  accomplish  the  entire  n^trix  solution;  these  forms 
are; 

A  *  A/B  and  A  =  A  -  B*C 

for  different  values  of  A,  B,  and  C.  The  technique  used  in 
SPICE2  is  to  perform  a  mock  equation  solution,  generating  the  ap¬ 
propriate  machine  instructions  in  place  of  performing  the  actual 
arithmetic.  Subroutine  CODEXC  is  used  as  an  interface  to  the 
generated  machine  code,  since  standard  FORTRAi-^  contains  no  provi¬ 
sion  for  calling  suoroutines  which  have  not  been  previously  com¬ 
piled  . 

Timing  comparisons  indicate  that  the  code  generated  for  the 
CDC  6400  computer  at  the  University  of  California  at  Berkeley  ex¬ 
ecutes  between  3  and  4  times  faster  than  the  equivalent  PORTRAIT 
subroutine  (which  must  contain  additional  code  to  again  map  (i,^) 
matrix  entries  to  their  location  within  the  LVN  taole  as  the 
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equations  are  solved).  However,  the  code  can  occupy  a  conslder- 
aole  a.Tiount  or  memory;  for  the  UA741,  approximately  2000  woras. 
Whether  or  not  to  use  the  generated  raachine  code  depends  on 
memory  availability  and  computer  resource  charging  algorithms. 
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8  Analysis 

This  section  briefly  describes  the  element  loadi  equation 
solution!  and  junction  limiting  algorithms  in  SPIC£2« 

6.1  Element  Load 

The  construction  of  the  circuit  equation  coefficient  matrix 
can  ne  performed  incrementally i  one  element  at  a  time  (seer  for 
exampler  (81).  A  concise  description  of  where  each  element  type 
adds  to  the  matrix  can  be  found  in  Section  13.6. 

8.2  Matrix  Equation  Solution 

The  solution  technique  for  the  system  of  linearized  circuit 
equations  is  to  first  perform  an  LU  factorization  of  the  coeffi¬ 
cient  matrix#  followed  by  forward  and  backward  substitution  steps 
to  solve  for  a  particular  right-hand  side  forcing  vector.  This 
technique  has  two  main  advantages  over  simple  Gaussian  elimina¬ 
tion.  First  of  all#  multiple  right-hand  sides  can  be  solved  with 
much  less  effort#  since  the  L  and  U  factors  need  to  be  computed 
only  once.  Secondly#  SPICE2  uses  the  adjoint  analysis  technique 
19]  to  efficiently  perform  small-signal  noise  and  distortion#  and 
dc  sensitivity  analyses*  The  equation  matrix  of  the  adjoint  cir¬ 
cuit  is  the  transpose  of  the  matrix  for  the  original  circuit; 
conveniently#  the  L  and  U  factors  of  the  transposed  matrix  are 
the  transpose  of  the  original  L  and  U  factors.  Thus#  the  adjoint 
analyses  may  be  performed  for  just  the  cost  of  a  forward  and 
backward  substitution  step  using  the  transpose  of  the  already- 
factored  equation  matrix- 
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8.2*1  LU  Decotnpositon 

As  noted*  SPICE2  factors  the  coefficient  matrix  into  lower- 
and  upper-triangular  matrices  L  and  U.  The  actual  algorithm  used 
for  the  factorization  is  so  arranged  that  the  L  and  U  factors  oc¬ 
cupy  the  same  memory  locations  as  the  original  matrix*  The  algo¬ 
rithm  is  given  below: 


Transform  "Ax  =  b"  into  "LUx  *  b" 

NSTOP  =  number  of  circuit  equations;  note  that  the  ground 
node  is  included  in  the  matrix  but  not  solved* 

for  (I  =  2  to  NSTOP) 

{  for  (J  *  I+l  to  NSTOP) 

{  A(J,I)  =  A{J,I)/A(I,I); 
for  (K  =  I+l  to  NSTOP) 

{  A(J,K)  «  A(J,K)  -  A(J,I)*A(I,K) ;  > 

> 

> 

Result:  A(I,J)  «  U{I*J)  for  J  >  I 

«  L{I,J)  for  J  <  I 


8*2*2  Forward/Backward  Substitution 

The  algorithm  used  in  SPICE2  to  perform  the  forward  and 
backward  substitution  steps  follows: 


Forward  substitution: 

-1 

evaluate:  y  »  Ux  *  L  b 

for  (I  =  2  to  NSTOP) 

{  for  (J  «  I+l  to  NSTOP) 

{  B(J)  »  B(J)  -  L(J*I)*B(I);  } 

> 
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Backward  substitution: 

-1  -1  -1 

evaluate:  x*U  y»U  L  b 

for  (I  «  NSTOP  to  2  oy  -1) 

{  for  (J  «  I+l  to  NSTOP) 

{  B(I)  «  B(I)  -  U(I,J);  } 

B(I)  «  B(I)/U(I,I); 

} 

8 >3  Devices 

SPICE2  contains  built-in  models  for  the  p-n  junction  diode# 
the  bipolar  junction  transistor  (BJT)#  the  junction  field-effect 
transistor  (JFET)#  and  the  metal-oxide-seraiconductor  field-effect 
transistor  (MOSFET).  The  model  equations  for  the  diode#  BJT#  and 
JFET  are  described  in  [1]. 

8.3.1  MOSFET  Model 

The  MOSFET  model  used  in  SPICE2  is  derived  from  the 
Frohman-Grove  model  [10] ,  with  the  addition  of  some  subthreshold 
[11]  and  short-channel  [12]  effects.  The  nonlinear  oxide  capaci¬ 
tance  formulation  of  J.  Meyer  [13]  is  included.  The  actual  equa¬ 
tions  used  are  given  in  Section  13.5. 

8.3.2  Junction  Limiting 

The  Newton-Raphson  algorithm  used  to  aid  convergence  can  run 
into  numerical  problems  due  to  the  exponential  nonlinearities  in 
the  semiconductor  device  models.  To  solve  this  difficulty#  addi¬ 
tional  constraints  are  placed  on  the  per-iteration  change  in  dev¬ 
ice  voltages.  Subroutine  PNJLIM  contains  (and  most  concisely 
describes)  the  limiting  algorithm  used  for  p-n  junction  voltages; 
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subroutine  FETLIH  is  used  to  limit  per-iteration  changes  in 
voltages . 


FET 


-32- 


SPICE2  REPORT 


DCTRAiJ  OVERLAY 


Ii'JTRQLjJCTIOW 


9  Dctran  Overlay 


The  DCTRAi'J  overlay  is  the  largest  one  in  SPICE2;  it  per¬ 
forms  the  dc  transfer  curve r  dc  operating  point*  initial  tran¬ 
sient  operating  point*  and  transient  analyses.  The  overlay  con¬ 
sists  of  routines  DCTRAN*  COMCOF*  TRUNC*  TERR*  SORUPD*  ITERS* 
CODEXC*  LOAD,  NLCSRC  *  UPDATE*  EVPOLY  *  EVTER.-l,  NXTPWR*  1NTGR8, 
PNJLIM,  DIODE*  BJT*  FETLIM*  JFET*  MOSFET*  MOSE^N *  and  MOSCAP. 


9.1  DC  Operating  Point 

The  dc  operating  point  computation  is  indicated  when  vari- 
aoles  MODE  and  MODEDC  have  values  1  and  1,  respectively.  A 
flowchart  for  this  analysis  followsJ 


initialize; 

TIME  *  0; 

call  SORUPD  to  set  sources  to  time-zero  values; 
liMITF  =  2; 
call  ITERS; 
if  (converged) 

{  INITf  =  4; 

call  DIODE;  call  BJT;  call  JFET;  call  MOSFET; 
print  operating-point  solution; 

} 


The  actual  Newton-Raphson  iteration  is  controlled  oy 
suoroutine  ITERS.  A  flowchart  for  that  suoroutine  follows: 


IGOOF  =  ITERNO  =  NDRFLO  =  NONCON  =  0; 

DONE  =  .false.; 
while  (not  DONE) 

{  call  LOAD; 

if  ( (NOSOLV  is  nonzero)  and 

(analysis  *  initial  transient))  exit; 

ITERNO  =  ITERNO  +  1; 
switch  (INITF)  of 

{  "1“:  if  (NONCON  =  0)  exit;  goto  solve; 

"2":  INITF  =  3;  goto  solve; 

"3'*:  if  (NONCON  =  0)  INITF  =  1;  goto  solve; 
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“4", -5", "6":  INITf  «  1; 
} 

solve: 


} 


if  (ITERNO  >  iteration  limit)  exit; 
call  DCDCMP;  call  DCSOL; 
if  (IGOOF  >  0) 

{  IGOOF  »  0;  NDRFLO  «  NDRFLO  +1;  } 

if  (NONCON  «  0) 

<  check  node  voltages  for  convergence; 
if  (voltages  converged)  DONE  »  .true.; 


> 


At  the  conclusion  of  the  dc  operating-point  analysis  (as  a 
result  of  the  calls  to  the  device  modeling  routines  with  INITF 
set  to  4)/  the  LXO  table  contains  all  the  small-signal  device 
parameters f  which  are  subsequently  printed  out  by  subroutine  DCOP 
(in  the  DCOP  overlay)* 


9-2  Transient  Initial-Conditions 

The  transient  initial-condition  computation  establishes  a 
consistent  set  of  circuit  conditions  prior  to  the  start  of  the 
transient  analysis.  (By  "consistent"  is  meant  conditions  which 
satisfy  Kirchboff's  laws.)  The  DCTRAN  overlay  performs  this 
analysis  when  MODE  and  MODEDC  have  values  of  1  and  2,  respective¬ 
ly.  A  flowchart  for  the  initial  transient  solution  algorithm 
follows : 


initialize; 

TIME  =  0; 

call  SORUPD  to  set  sources  to  time-zero  values; 
INITF  =2; 
call  ITERS; 
if  (converged) 

{  print  solution;  } 
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9.3  DC  Transfer  Curves 

The  dc  transfer  curve  computation  is  simply  a  repetitive  dc 
operating  point  computation  performed  for  a  range  of  values  for 
one  independent  voltage  or  current  source  in  the  circuit.  The 
analysis  is  performed  when  MODE  and  MODEDC  have  values  of  1  and 
3*  respectively.  A  flowchart  for  this  analysis  follows: 


initialize; 

TIME  =  0; 

call  SORUPD  to  set  sources  to  time-zero  values; 


INITF  =  2; 

for  (each  source  value) 

{  call  ITERS; 

if  (not  converged) 
store  outputs; 
INITF  =  6; 


> 


stop  analysis; 


(the  OVTPVT  overlay  is  called  from  the  overlay  root 
to  output  any  requested  plots) 


The  value  of  b  for  INITF  causes  a  linear  extrapolation  from 
the  previous  solutions  to  obtain  the  initial  guess  used  for  the 
subsequent  iteration;  that  is*  the  initial  guess  for  sweep  point 
n+1  is  the  linear  extrapolation  of  the  solutions  at  sweep  points 
n-1  and  n. 


9-4  Transient  analysis 

The  transient  analysis  is  performed  when  MODE  has  a  value  of 
2.  A  flowchart  of  the  analysis  steps  follows: 


initialize;  TIME  =  0;  DELTA  =  TSTEP;  INITF  =  5; 
savout:  store  outputs  in  LOUTPT  table; 
if  (TIME  >  TSTOP)  exit; 

adjust  DELTA  for  breakpoint  table  values; 
newtim:  TIME  ®  TIME  +  DELTA; 
call  SORUPD; 

if  (INITF  not  equal  to  5)  INITF  *  6; 
call  ITERS; 
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if  (converged)  goto  tsterr; 

TIME  «  TIME  -  DELTA;  DELTA  «  DELTA/8;  goto  tstdel; 
tsterr:  call  TRUi'lC; 

if  (error  acceptable)  goto  savout; 

TIi*lE  «  TIME  -  DELTA; 

DELTA  «  DELNEW  (computed  in  TRUNC); 
tstdel:  if  (DELTA  <  DELMIH)  Stop  analysis; 
goto  newtim; 

9«4«1  Numerical  Integration  Methods 

SPICE2  uses  either  Trapezoidal  (14]  or  variable-order  Gear 
[15]  to  perform  the  numerical  integration  necessary  for  the  tran¬ 
sient  analysis-  The  choice  is  based  on  the  value  of  the  COMMON- 
variable  METHOD- 

9-4-2  Automatic  Timestep 

There  are  two  methods  available  to  the  user  to  control  the 
computation  of  the  appropriate  transient  analysis  timestep: 
iteration  count  and  truncation  error  estimation-  The  iteration 
count  method  uses  the  number  of  Newton-Raphson  iterations  re¬ 
quired  to  converge  at  a  timepoint  as  a  measure  of  the  rate  of 
change  of  the  circuit-  If  the  number  of  iterations  is  less  than 
ITL3 >  the  timestep  is  doubled;  if  greater  than  ITL4 >  the 
timestep  is  divided  oy  8-  This  method  has  the  advantage  of 
minimal  computation  overhead;  however,  it  suffers  in  that  it 
does  not  take  into  account  the  true  rate-of-change  of  circuit 
variables  - 

The  second  method  of  timestep  determination  is  based  on  the 
estimation  of  local  truncation  error-  The  appropriate  function 
derivative  is  approximated  by  a  divided  difference  -  hence  the 
"estimation- "  Although  requiring  somewhat  more  cpu  time  than 
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iteration  count#  the  truncation  error  estimation  method  allows  a 
meaningful  error-hound  to  he  set  on  the  computed  output  values. 
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10  Acan  Overlay 


The  ACAN  overlay  performs  the  small-signal •  noise#  and  dis¬ 
tortion  The  overlay  contains  routines  ACAN#  ACDCMP#  ACSOL# 
ACLOAD#  NOISE#  ACASOL#  DINIT#  and  DISTO.  All  analyses  performed 
use  the  linearized  models  for  all  the  nonlinear  devices*  A 
flowchart  of  this  overlay  follows: 


initialize; 

for  (FREQ  =  FSTART  to  FSTOP ) 

{  call  ACLOAD; 
call  ACDCMP; 
call  ACSOL; 

save  small-signal  outputs; 

if  (noise  analysis  requested)  call  NOISE; 

if  (distortion  analysis  requested)  call  DISTO; 

} 

(the  OVTPVT  overlay  is  called  from  the  overlay  root 
to  output  any  requested  plots) 


10*1  Small-Signal  Analysis 

The  small-signal  analysis  is  performed  at  each  frequency 
point  by  first  calling  subroutine  ACLOAD*  This  subroutine  zeroes 
and  then  loads  the  complex-valued  equation  matrix  using  the 
results  of  the  dc  operating  point  analysis  stored  in  the  LXO 
table*  Then  subroutines  ACDCMP  and  ACSOL  are  called,  which  per¬ 
form  an  LU  decomposition  and  f orward/bacJcward  substitution  on  the 
set  of  circuit  equations  to  obtain  its  solution. 


10*2  Noise  and  Distortion  Analyses 

The  noise  analysis  used  in  SPICE2  is  described  by  references 
116]  and  [17]*  The  distortion  analysis  algorithm  is  described  in 
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11  Dcop  Overlay 

The  DCOP  overlay  performs  three  functions:  it  prints  out 
the  device  operating  point  information*  computes  the  small-signal 
transfer  function*  and  computes  dc  sensitivities  of  output  vari- 
ables  with  respect  to  circuit  parameters*  The  overlay  contains 
routines  DCOP*  SSTF*  SENCAL*  and  ASOL* 

11.1  Small-Signal  Transfer  Function 

The  small-signal  transfer  function  analysis  is  performed  by 
subroutine  SSTF;  it  evaluates  the  gain*  input  resistance*  and 
output  resistance  of  a  specified  two-port  using  the  linearized 
device  models  computed  as  a  part  of  the  dc  operating  point 
analysis. 

11.2  DC  Sensitivities 

The  dc  sensitivities  of  output  variables  with  respect  to 
circuit  parameters  are  computed  oy  subroutine  SENCAL.  The 
analysis  is  performed  using  the  adjoint  technique  t9] *  thus  re¬ 
quiring  only  one  equation  solution  per  output  variable  (rather 
than  per  circuit  parameter). 
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12  OVTPVT  overlay 

The  OVTPVT  overlay  generates  the  taoular  listings  and  line- 
printer  plots  displaying  the  results  of  the  multi-point  analyses* 
It  contains  routines  OVTPVT,  NTRPL8 ,  SETPRN,  SETPLT,  PLOT,  SCALE, 
and  FOURAN*  The  overlay  is  called  after  each  multi-point 
analysis,  to  avoid  the  necessity  of  storing  unnecessary  output 
values  in  memory. 

Subroutine  FOURAN  is  used  to  calculate  the  Fourier  coeffi¬ 
cients  of  transient  analysis  output  waveforms.  The  coefficients 
are  computed  recursively  using  a  "truncated  infinite  series"  al¬ 
gorithm  [1,  pp.  254-259]. 
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USERFS  GUIDE  FOR  SPICE 


SPICE  IS  A  GENERAL-PURPOSE  CIRCUIT  SIMULATION  PROGRAM  E OR  NONLINEAR  DC, 
NONLINEAR  transient,  and  LINEAR  AC  ANALYSES#  CIRCUITS  MAY  CONTAIN  RESISTORS, 
CAPACI’’0RS,  INDUCTORS,  MUTUAL  INDUCTORS,  INDEPENDENT  VOLTAGE  AND  CURRENT 
SOURCES,  POUR  types  OF  DEPENDENT  SOURCES,  TRANSMISSION  LINES,  AND  THE  =OUP  MOST 
COMMON  SEMICONDUCTOR  DEVICES:  DIODES,  BJTS,  JFFts,  AND  MOSFETS. 

S=ICE  HAS  SUILT-IN  models  for  the  semiconductor  devices,  AND  THE  USER  NEED 
SPECIFY  ONLY  THE  PERTINENT  MODEL  PARAMETER  VALUES*  THE  MODEL  FOR  THE  BJT  IS 
BASED  ON  the  integral  CHARGE  MODEL  OF  GUMMSL  AND  PCON ;  HOWEVER,  I®  THE  GUMMEL- 
PCON  PARAMETERS  ARE  NOT  SPECIFIED,  THE  MODEL  REDUCES  TO  THE  SIMPLER  EBERS-MCLL 
MODEL.  IN  EI-hER  case,  CHARGE  STORAGE  EFFECTS,  OHMIC  RESISTANCES,  AND  A 
C UPRENT-DEPENDENT  output  CONDUCTANCE  MAY  BE  INCLUDED,  THE  DIODE  MODEL  CAN  BE 
USED  FOR  EI’HER  JUNCTION  DIODES  OR  SCHOTTKY  BARRIER  DIODES.  THE  J^ET  MODEL  IS 
BASED  ON  •'HE  MODEL  OF  SHICHMAN  AND  HODGES.  THE  MODEL  FOR  THE  MOSFET  IS 

based  on  the  frohman-grove  model;  however,  channel-length  modulation,  sub- 
ThRESHOLD  conduction,  AMD  SOME  SHOP T- CHANN EL  EFFECTS  ARE  INCLUDED. 


PROGRAvi  limitations 


SPICE  USES  DYNAMt;;  MEMORY  MANAGEMENT  TO  STORE  ELEMENTS,  MODELS,  AND  OUTPUT 
VALUES.  ”HUS ,  the  only  LIMITATION  IMPOSED  BY  THE  PROGRAM  ON  THE  SIZE  OR  THE 
complexity  nr  twE  CIRCUIT  TO  BE  SIMULATED  I-S  THAT  ALL  NECESSARY  DATA  FIT  IN 
MFynOY,  FQR  EXAMPLE,  A  lOO-PQINT  TRANSIENT  ANALYSIS  0=  THE  UA*’A1  OPERATIONAL 
amplifier  REQUIRES  A PORQX I  MAT  ELY  BCOOO  (OCTAL)  WORDS  OF  MEMORY  ON  THE  CDC  6*00 
AT  THE  university  of  CALIFOPNIA  AT  BERKELEY. 


transient  analysts 


■^ME  transient  analysts  oortidn  of  spice  computes  the  transient  output 

VARIABLES  as  A  ^UNCTION  OF  TIME  OVEP.  A  U  SEP -S  PEC  I  ^  I  ED  TIME  INTERVAL.  THE 
initial  CONt^ITIONS  ARE  AUTOMATICALLY  DETERMINED  BY  A  DC  ANALYSIS.  ALL  SOURCES 
WHICH  ARE  NOT  TIME  DEPENDENT  (FOP  EXAMPLE,  POWER  SUPPLIES)  ARE  SET  TO  THEIR  DC 
value.  for  LARGE-SIGNAL  SINUSOIDAL  SIMULATIONS,  A  FOURIER  ANALYSIS  OF  THE 
OUTPUT  WAVEFORM  CAN  EE  SPECIFIED  TO  OBTAIN  THE  FREQUENCY  DOMAIN  FOURIER  COEF- 
=  ICIENTSo  tme  transient  TIME  INTERVAL  AND  THE  FOURIER  ANALYSIS  OPTIONS  ARE 
SPECI'^IED  ON  THE  .TRAN  AND  .FOURIER  CONTROL  CARDS. 

-  ANALYSIS  AT  different  TEMPERATURES 

ALL  INOU"^  DATA  FOR  SPICE  IS  ASSUMED  TO  HAVE  BEEN  MEASURED  AT  27  DEG  C 
(300  DEG  K).  THE  SIMULATION  ALSO  ASSUMES  A  NOMINAL  TERMPERATURE  OF  27  DEG  C. 
THE  CIRCUIT  CAN  BE  SIMULATED  AT  OTHER  TEMOERATURES  BY  USING  A  .TEMP  CONTROL 
C  ARD. 

•'EWRERATURE  APPEARS  EXPLICITLY  IN  THE  EXPONENTIAL  TERMS  0=^  THE  BJT  AND 
DIODE  MODEL  EQUATIONS.  IN  ADDITION,  SATURATION  CURRENTS  HAVE  A  BUILT-IN 
TEMRERATU=E  dependence.  the  temperature  dependence  of  the  SATURATION  CURRENT  ' 
IN  THE  BJT  MODELS  IS  DETERMINED  BY: 

ISC'EMD)  =  IOW(TEMP  *PT  )  *EXP(-Q*EG/(  <*TEMP  )  ) 

WHERE  K  IS  B0LT2MANS  CONSTANT,  0  IS  THE  ELECTRONIC  CHARGE,  10  IS  A  CONSTANT,  EG 
IS  the  ENERGY  GAP  WHICH  IS  A  MODEL  PARAMETER,  AND  PT  IS  THE  SATURATION  CURRENT 
TEMPERATURE  EXPONENT  (ALSO  A  MODEL  PARAMETER,  AND  USUALLY  EQUAL  TO  3).  THE 
temperature  dependence  of  the  SATURATION  CURRENT  IN  THE  JUNCTION  DIODE  MODEL  IS 

determined  by: 

IS(TEMP)  =  IOF(TEMP*( PT/N) )*EX0(-Q*EG/(K*N*TEMP) ) 

WHERE  N  IS  THE  EMISSION  COEFFICIENT,  WHICH  IS  A  MODEL  PARAMETER,  AND  THE  OTHER 
SYMSHlS  have  the  same  MEANING  AS  ABOVE.  NOTE  THAT  FOR  SCHCTTXY  BARRIER  DIODES, 
THE  VALUE  OF  THE  SATURATION  CURRENT  TEMPERATURE  EXPONENT,  PT ,  IS  USUALLY  2. 

temperature  APPEARS  EXPLICITLY  IN  THE  VALUE  OF  JUNCTION  POTENTIAL,  PHI, 

‘=’0R  ALL  the  device  MODELS.  THE  t^mpepaTURE  DEPENDENCE  IS  DETERMINED  BY: 

PHI(TEMP)  =  K*-EMP/Q*LOG(NABND/'NI  (  TEMP) +2  ) 

W-iERE  K  IS  POLTZMANS  CONSTANT,  Q  IS  THE  ELECTRONIC  CHARGE,  NA  IS  THE  ACCEPTOR 
impurity  density,  ND  is  the  donor  impurity  DENSITY,  NI  IS  THE  INTRINSIC  CONCEN- 
•RA’ION,  AND  EG  IS  THE  ENERGY  GAP. 

temperature  APPEARS  EXPLICITLY  IN  THE  VALUE  OF  SURFACE  MOBILITY,  UO,  FQR 

THE  MOSFET  model.  the  temperature  dependence  IS  determined  by: 

UO(TEMP)  =  UC(TNCM)/(TEMO/TNOM) ♦(! ,5) 

temperature  APPEARS  EXPLICITLY  fqr  RESISTOR  VALUES  ACCORDING  TO  THE 
FOLLOWING  expression: 

value  (  TEMP)  =  value  (TNOM)  *(  1  +  TCl  *(  TEm  P-TNOM  >  •*■  TC2*  (  T  E  MP-TNOM  )  *2  )  )  , 

V 

WHERE  t£mP  is  the  CIRCUIT  TEMPERATURE,  TNOM  IS  the  NOMINAL  TEMPERATURE,  AND 
TCI  AND  TC2  ARE  the  first-  AND  SECOND-ORDER  TEMPERATURE  COEFFICIENTS. 


TYB£S  OF  analysis 


- DC  ANALYSIS 

THE  DC  ANALYSIS  »09TICN  OF  SPICE  DETERMINES  THE  DC  OPERATING  OOINT  of  THE 
CIRCUIT  with  inductors  SHORTED  AND  CAPACITORS  OPENED.  A  DC  ANALYSIS  IS  AUTO¬ 
MATICALLY  OERFORMED  RPIOP  tQ  A  TRANSIENT  ANALYSIS  TO  DETERMINE  TWE  TRANSIENT 
initial  CONDITIONSf  AND  PRIOR  TO  AN  AC  SMALL-SIGNAL  ANALYSIS  TO  DETERMINE  THE 
LINEARIZED,  small-signal  mqDELS  FOR  NONLINEAR  DEVICES.  IF  REQUESTED,  THE  DC 
SMALL-SIGNAL  VALUE  OF  A  TRANSFER  FUNCTION  (RATIO  OF  OUTPUT  VARIABLE  TO  INPUT 
SOURCE),  INPUT  resistance,  AND  OUTPUT  RESISTANCE  WILL  ALSO  9E  COMPUTED  AS  A 
BART  0=’  THE  DC  solution.  ThE  DC  ANALYSIS  CAN  ALSO  BE  USED  TO  GENERATE  DC 
tranS'^ER  curves:  a  specified  INDEPENDENT  VOLTAGE  OR  CURRENT  SOURCE  IS  STEPPED 

OVER  A  user-specified  RANGE  AND  THE  DC  OUTbuT  VARIABLES  ARE  STORED  FOR  EACH 

sequential  source  value.  IR  requested,  SPICE  ALSO  WILL  determine  the  DC  small- 

signal  SENSITIVITIES  0=  SPECIFIED  OUTPUT  VARIABLES  WITH  RESPECT  TO  CIRCUIT 
PARAMETERS.  THE  DC  ANALYSIS  CPTionc  App  SPECIFIED  ON  THE  .DC,  .TF,  .OP,  AND 
.SENS  CON-ROL  CAROS. 


' - AC  SMALL-SIGNAL  ANALYSIS 

THE  AC  small-signal  PORTION  OF  SPICE  COMPUTES  THE  AC  OUTPUT  VARIABLES  AS  A 
PUNCTION  OF  B-PPQUPMCY.  THE  PROGRAM  FIRST  COMPUTES  THE  DC  OPERATING  POINT  OF 
the  CIRCUIT  AND  determines  LINEARIZED,  SMALL-SIGNAL  mqdELS  FOR  ALL  OF  THE  NON¬ 
LINEAR  DEVICES  IN  THE  CIRCUIT.  THE  RESULTANT  LINEAR  CIRCUIT  IS  THEN  ANALYZED 
OVER  A  USER-SPECIFIED  RANGE  OF  FREQUENCIES.  thE  DESIRED  OUTPUT  OF  AN  AC  SMALL- 
SIGNAL  ANALYSIS  IS  USUALLY  A  TRANSFER  FUNCTION  (VOLTAGE  GAIN,  TR A NS  I MPEOA NCE , 
ETC).  IF  THE  CIRCUIT  HAS  ONLY  ONE  AC  INPUT,  IT  IS  CONVENIENT  TO  SET  THAT  INPUT 
TO  UNITV  and  ZERO  PHASE,  SO  THAT  OUTPUT  VARIABLES  HAVE  THE  SAME  VALUE  AS  THE 
TRANSFER  function  OF  THE  OUTPUT  VARIABLE  WITH  RESPECT  TO  THE  INPUT. 

the  GENERA-'ION  of  white  noise  by  RESISTORS  AND  SEM I CONO  UCTCP  DEVICES  CAN 
ALSO  BE  SIMULATED  WITH  THE  AC  SMALL-SIGNAL  PORTION  OF  SPICE.  EQUIVALENT  NOISE 

SOURCE  values  ape  determined  automatically  from  the  small-signal  operating 
POINT  0^  the  circuit,  and  the  contribution  QR  each  noise  source  is  ADDED  AT  A 
given  summits  point,  the  total  output  noise  level  AND  THE  EQUIVALENT  INPUT 
NOISE  LEVEL  ARE  DETERMINED  A~  EACH  FREQUENCY  OQINT.  THE  OUTPUT  AND  INPUT  NOISE 
LEVELS  ARE  normalized  WITH  RESPECT  TO  THE  SCUAPE  ROOT  OR  THE  NOISE  BANDWIDTH 
ANO  have  the  units  VOLTS/PT  HZ  OR  AMPS/RT  HZ.  THE  OUTPUT  NOISE  AND  EQUIVALENT 
INPUT  NOISE  CAN  BE  PRINTED  OR  PLOTTED  IN  Tm=  SAME  =ASH10N  AS  O'OHER  OUTPUT 
VARIABLES.  NO  ADDITIONAL  INPUT  DATA  IS  NECESSARY  FOR  THIS  ANALYSIS. 

flicker  noise  SOURCES  CAN  BE  SIMULATED  IN  THE  NOISE  ANALYSIS  BY  INCLUDING 
VALUES  FOR  THE  PARAMETERS  KF  ANO  AF  ON  THE  APPROPRIATE  DEVICE  MODEL  CAROS. 

the  otstop-tion  characteristics  of  A  CIRCUIT  IN  The  small-signal  moos  can 
BE  simulated  as  a  BART  OF  THE  AC  SMALL-SIGNAL  ANALYSIS.  THE  ANALYSIS  IS 
pERFORmso  assuming  that  CNF  OR  TWO  signal  frequencies  are  IMPOSED  AT  the  INPUT, 

the  TpCQuENCY  range  and  the  NOISE  AND  distortion  ANALYSIS  PARAMCjERS  ARE 
SPECIFIED  0>J  the  ,4C,  .noise,  AND  eCIS^CRTION  CONTROL  CARDS. 


C  0MVERG5NCE 


BOTH  DC  AND  TPANSIENT  SOLUTIONS  ARE  OBTAINED  BY  AN  ITERATIVE  PROCESS  WHICH 
IS  TERMINATED  WHEN  BOTH  OF  yhE  FOLLOWING  CONDITIONS  HOLD: 

1)  THE  NONLINEAR  BRANCH  CURRENTS  CONVERGE  TO  WITHIN  A  TOLERANCE  OF 
0.1  PERCENT  OR  1  oi COAMP  (l.OE-12  AMP),  WHICHEVER  IS  LARGER 

2)  THP  node  voltages  CONVERGE  TO  WITHIN  A  TOLERANCE  OF  0» 1  PERCENT 
OR  1  MICROVOLT  (l,0E-6  VOLT),  WHICHEVER  IS  LARGER 

although  the  algorithm  used  in  SPICE  mas  been  found  to  be  very  reliable,  in 
SOME  CASES  IT  WILL  FAIL  TO  CONVERGE  TO  A  SOLUTION.  WHEN  THIS  FAILURE  OCCURS, 
-HE  PROGRAM  WILL  PRINT  THE  NODE  VOLTAGES  AT  THE  LAST  ITERATION  AND  TERMINATE 
yhE  JOB,  IN  SUCH  CASES,  THE  NODE  VOLTAGES  THAT  ARE  PRINTED  ARE  NOT  NECESSARILY 
CORRECT  OR  EVEN  CLOSE  TO  THE  CORRECT  SOLUTION. 

failure  to  converge  in  the  DC  analysis  is  usually  DUE  TO  AN  ERROR  IN 
specifying  CICCUIT  CCNNECTICNS,  element  values,  or  MODEL  PARAMETER  VALUES. 
REGENERATIVE  SWITCHING  CIRCUITS  OR  CIRCUITS  WITH  POSITIVE  FEEDBACK  PROBABLY 
WILL  NOT  CONVERGE  IN  THE  DC  ANALYSIS  UNLESS  THE  POFF*  OPTION  IS  USED  FOR  SOME 
OF  THE  DEVICES  IN  THE  FEEDBACK  PATH. 


INPUT  FORMAT 


THE  INPUT  PORMAT  FOR  SPICE  IS  OF  THE  FPEE  FORMAT  TYPE.  FIELDS  ON  A  CARO 
APE  separated  by  ONE  OR  MOPE  BLANKS,  A  COMMA,  AN  EQUAL  (=)  SIGN,  OR  A  LEFT  OR 
RIGHT  PARENTHESIS)  EXTRA  SPACES  ARE  IGNORED.  A  CARD  MAY  BE  CONTINUED  BY 
PUNCHING  A  •*.  (PLUS)  IN  COLUMN  1  OF  THE  FOLLOWING  CARD;  SPICE  CONTINUES  READING 
BEGINNING  WITH  COLUMN  2, 

A  NAME  field  MUST  BEGIN  WITH  A  LETTER  (A  THROUGH  Z)  AND  CANNOT  CONTAIN 
ANY  delimiters.  only  THE  FIRST  EIGHT  CHARACTERS  OF  THE  NAME  ARE  USED. 

A  NUMOFR  RIELD  MAY  BE  AN  INTEGER  FIELD  (12,  -AA )  ,  A  FLOATING  POINT  R I  EL D 
(3.14159),  EITHER  AN  INTEGER  OR  FLOATING  POINT  NUMBER  FOLLOWED  BY  AN  INTEGER 
EXPONENT  (lE-14,  2.6EE3),  OR  EITHER  AN  INTEGER  OR  A  FLOATING  POINT  NUMBER 

FOLLOWED  BY  ONE  OF  THE  FOLLOWING  SCALE  FACTORS; 


G=lEa  MEGS1E6  K=1E3  MILs25.AE-4  M=lE-3  U=lE-6  N=lE-9  P=1E-12 

LETTERS  IMMEDIATELY  FOLLOWING  A  NUMBER  THAT  ARE  NOT  SCALE  FACTORS  ARE  IGNORED, 
AND  LETTERS  IMMEDIATELY  FOLLOWING  A  SCALE  FACTOR  ARE  IGNORED.  HENCE,  10,  lOV, 
lOVOLTS,  AND  lOHZ  ALL  REPRESENT  -HE  SAME  NUMBER,  AND  M,  MA,  MSEC,  AND  MMHOS  ALL 
REPRESENT  the  same  SCALE  FACTOR.  NO-E  THAT  1000,  1000.  0,  lOOOHZ,  1 E3 ,  1.0E3, 

IKHZ,  AND  IK  ALL  REPRESENT  THE  SAME  NUMBER. 


CIRCUIT  O^SCPIBTION 


the  CIOC'JIT  to  be  analyzed  is  described  to  SOICE  BY  A  SET  0^  ELEMENT 
CAPOS,  WHICH  DE'=’INE  THE  CIPCUIT  TOPOLOGY  AND  ELEMENT  VALUES,  AND  A  SE^  OP 
CONTPOL  CAPOS,  WHICH  DEFINE  THE  MODEL  PARAMETERS  AND  THE  RUN  CONTROLS*  THE 
PIPST  CA»0  IN  THE  INPUT  DECK  MUST  BE  A  TITLE  CAPO,  AND  THE  LAST  CARD  MUST  BE  A 
•  END  CARD.  THE  ORDER  O’^  THE  REMAINING  CAPOS  IS  ARBITRARY  (EXCEPT,  CP  COURSE, 
that  continuation  CAPOS  MUST  immediately  POLLOW  the  CARO  BEING  CONT  INUEO > • 

EACH  ELEMENT  IN  THE  CI®CUIT  IS  SPECIFIED  BY  AN  ELEMENT  CARO  THAT  CONTAINS 
THE  element  name,  THE  CIRCUIT  NODES  TO  WHICH  THE  ELEMENT  IS  CONNECTED,  AND  yh£ 
values  op  the  PARAMETERS  that  OEYERMINE  THE  ELECTRICAL  CHARACTERISTICS  Op  THE 
element.  the  ®IRST  letter  CP  THE  ELEMENT  NAME  3PECIP1ES  the  ELEMENT  TYOE. 
the  FORMAT  FOP  the  spice  ELEMENT  TYPES  IS  GIVEN  IN  WHAT  FOLLOWS*  THE  STRINGS 
wxxxxxxxs,  lYYYYYYY*,  AND  *ZZ2Z22Zw  DENOTE  ARBITRARY  ALPHANUMERIC  STRINGS*  FOR 
example,  A  PESISTOR  NAME  MUST  BEGIN  WITH  THE  LETTER  R  AND  CAN  CONTAIN  FROM  ONE 
TO  eight  CHARACTERS*  HENCE,  R,  R1 ,  RSE,  ROUT,  AND  R3AC22Y  ARE  VALID  RESISTOR 
NAMES. 


DATA  =IELDS  THAT  APE  ENCLOSED  IN  SQUARE  BRACKETS  (  [  ]  )  ARE  OPTIONAL*  ALL 

INDirA-'ED  PUNCTUATION  (PARENTHESES,  EQUAL  SIGNS,  ETC.)  ARE  REQUIRED.  WITH 
RESPECT  TO  BRANCH  VOLTAGES  AND  (DJPRENTS,  SPICE  UNIFORMLY  USES  THE  ASSOCIATED 
o£=EPENCE  CONVENTION  INDICATED  IN  THE  DRAWING  BELOW; 


I 


> 

>> 

>>>>>>>>>>> 

>> 

> 


NODES  MUST  BE  NONNEGATIVE  INTEGERS  BUT  NEED  NOT  BE  NUMBERED  SEQUENTIALLY, 
the  datum  (GROUND)  NODE  MUST  BE  NUMBERED  ZERO.  THE  CIRCUIT  CANNOT  CONTAIN  A 
_00P  0=  VOLTAGE  SOURCES  AND/OR  INDUCTORS  AND  CANNOT  CONTAIN  A  CUTSET  OF  CURRENT 
SOURCES  ANO/CR  capacitors.  EACH  NODE  IN  THE  CIRCUIT  MUST  havE  A  DC  RATH  TO 
GROUND,  EVE»Y  NODE  MUST  HAVE  AT  LEAST  TWO  CONNECTIONS  EXCEPT  FOR  TRANSMISSION 
LINE  NODES  (TO  PERMIT  UN- ERMI NA -^ED  TRANSMISSION  LINES). 


PLEMEN”  CAPOS 


RESISTORS 


GENERAL  POPM  OXXXXXXX  N1  N2  VALUE  [T C= T Cl [ , T C 2  ]  ] 

EXAMdlE  RCl  12  l'^  IK  TC  =  0.001 , 0.  015 

N1  AND  N2  APE  the  two  ELEMENT  NODES.  VALUE  IS  THE  RESISTANCE  (IN  OHMS) 
AND  MAY  BE  =>0SI'IVE  OR  NEGATIVE  BUT  NOT  ZERO.  TCI  AND  TC2  ARE  THE  (OPTIONAL) 
thmpePA~UPE  COEpFIC  IENTS;  IF  NOT  SPECIFIED,  ZERO  IS  ASSUMED  FOR  BOTH,  THE 
VALUE  OP  the  resistor  AS  A  FUNCTION  OF  TEMPERATURE  13  GIVEN  BY: 

VALUE (TEMP)  =  VALUE(TNOM)*(  1 +TC 1 *( TEMP-TNOM )  +  TC 2* ( "E MP-TNOM ) *  2 )  ) 


CAPACITORS  AND  INDUCTOPS 


GENERAL  POPM  CXXXXXXX  N+  N-  ©0  [PI  ?2  •••]  [IC=1NC0ND] 

LYYYYYYY  N+  N-  PO  I  PI  P2  •••1  lIC  =  INCOND] 

examples  CBYO  13  0  lUF 

COSC  IT  23  lOU  5U  IC=3V 
LL INK  A2  69  lUH 

LSHUNT  23  51  lOU  22U  IC=15,7A 

N+  AND  N-  ARE  THE  POSITIVE  AND  NEGATIVE  ELEMENT  NODES,  RESPECTIVELY. 

°0,  (PI,  02,  ETC.)  ARE  THE  COEFFICIENTS  OF  A  POLYNOMIAL  DESCRIBING  THE  ELEMENT 
value.  pop  CAOACT-^ORS,  the  CAPACITANCE  (IN  FARADS)  IS  EXPRESSED  AS  A  FUNCTION 
TF  the  voltage  ACROSS  T hE  ELEMENT,  FOR  INDUCTORS,  THE  INDUCTANCE  (IN  HENRIES) 
IS  EXPRESSED  AS  A  FUNCTION  OP  THE  CURRENT  THROUGH  THE  ELEMENT,  TO  ILLUSTRATE, 
“HE  SECOND  example  above  DESCRIBES  A  CAPACITOR  WITH  A  VALUE  DEFINED  BY 

C  =  lOE-6  +  5E-6FV 

WHERE  V  IS  the  VOLTAGE  ACROSS  THE  CAPACITOR, 

FOR  THE  CAPACITOR,  THE  (OPTIONAL)  INITIAL  CONDITION  IS  THE  INITIAL  (TIME- 
ZE=0)  value,  of  capacitor  voltage  (IN  VOLTS).  FOR  THE  INDUCTOR,  THE  (OPTIONAL) 
INITIAL  condition  IS  THE  INITIAL  (TIME-ZERO)  VALUE  OF  INDUCTOR  CURRENT  (IN 
AMPc)  that  FLOWS  PROM  N+,  THROUGH  THE  INDUCTOP,  TO  N-,  NOTE  THAT  THE  INITIAL 
CONDI” IONS  (!=  ANY)  APPLY  *ONLY»  IF  THE  UIC  OPTION  IS  SPECIpIED  ON  THE  .TRAN 
CARD, 


C0UOL£r  (MU'^UAL)  INOUCTCRS 


GENERAL  ^CRM  ^XXXXXXX  LYYYYYYY  L  ZZZZZZZ  VALUE 

EXAMdlES  <A3  LAA  LBB  0.999 

KXPRMP  LI  L2  0*87 


LYYYYYYY  AND  LZZZZZZZ  ARE  THE  NAMES  OF  THE  Y  v»0  COUPLED  INDUCTORS,  AND 
VALUE  IS  THE  coefficient  0*^  COUPLING,  <,  WHICH  MUST  EE  GREATER  THAN  0  AND  LESS 
than  DR  EQUAL  tq  i,  THE  COUPLED  INDUCTORS  CAN  BE  ENTERED  IN  EITHER  ORDER  AND 
MUST  BE  CONSTANT  (ONLY  PO  CAN  BE  SPECIFIED  ON  THE  INDUCTOR  CARDS).  IN  TERMS  OF 
THE  DOT  CONVENTION,  THIS  COUPLING  PUTS  A  DOT  BY  thE  FIRST  NODE  0^  EACH  OF  THE 
CCUFLED  INDUCTORS. 


TRANSMISSION  LINES  (LOSSLESS) 


GENERAL  FOPM  TXXXXXXX  N1  N?  NO  NA  ZOaVALUE  ITD«VAL'JE]  [FaFREO  [NL=NRmlEN]] 

+  [  IC  =  V1  ,  11 , V2 , 12 ] 


example 


Y1  1  0  2  0  20*50  TDaiONS 


N1  AND  N2  ARE  Tf-E  NODES  AT  OORT  1;  NO  AND  NA  ARE  THE  NODES  AT  PORT  2, 

20  IS  the  CHAOACTERISTIC  impedance.  the  length  of  the  line  MAY  BE  EXPRESSED  IN 
either  of  two  =0BMS.  THE  TRANSMISSION  DELAY,  TO,  MAY  BE  SPECIFIED  DIRECTLY  (AS 
TDaiONS,  =00  EXAMPLE).  alternatively,  A  FREQUENCY  F  MAY  BE  GIVEN,  TOGETHER 
WITH  NL ,  THE  normalized  ELECTRICAL  LENGTH  OF  THE  TRANSMISSION  LINE  WITH  RESPECT 
TO  the  wavelength  IN  "^HE  LINE  AT  ”hE  FREQUENCY  F,  tf  i,  FREQUENCY  IS  SPECIFIED 
BUT  NL  IS  omitted,  0.25  IS  ASSUMED  (THAT  IS,  THE  FREQUENCY  IS  ASSUMED  TO  BE  THE 
QUABTER-WAVE  FREQUENCY).  NOTE  THAT  ALTHOUGH  BOTH  FORMS  =0R  EXPRESSING  THE  LINE 
length  APE  INDICATED  AS  OPTIONAL,  ONE  0=  THE  TWO  *MUST*  EE  SPECIFIED. 

note  that  this  ELEMEN"  MODELS  ONLY  ONE  PROOAGATinG  MODE.  IF  ALL  FOuP 
NODES  APE  DISTINCT  IN  THE  ACTUAL  CIRCUIT,  THEN  TWO  MODES  MAY  BE  EXCITED.  TO 
SIMULATE  SUCH  A  SITUATION,  TwC  TR AN SM I S S I  0 N -L I NE  ELEMENTS  ARE  REQUIRED.  (SEE 
the  EXAMPLE  IN  APPENDIX  A  FC®  FURTHER  CLARIFICATION.) 

THC  (OPTIONAL)  INITIAL  CONDITION  S PE C 1  =  I  CAT  I  ON  CONSISTS  OF  THE  VOLTAGE 
AND  CURRENT  AT  EACH  OF  THE  TRANSMISSION  LINE  OQRTS.  NOTE  THAT  THE  INITIAL 
conditions  (1=  ANY)  APPLY  *ONLY*  IF  THE  UIC  OPTION  IS  SPECIFIED  ON  THE  .TRAN 
C  ARD. 


****  LINEAP  dependent  sources 


S=ICE  ALLOWS  CIRCUITS  TO  CONTAIN  LINEAR  DEPENDENT  SOURCES  CHARACTERIZED  BY 
ANY  OF  the  four  EQUATIONS 


I=GFV  VsE*V  I=F*I  V=H*I 

WHERE  G,  E,  F,  AND  H  ARE  CONSTANTS  REPRESENTING  T R ANS CONDUCT ANCE,  VOLTAGE  GAIN, 
CURRENT  GAIN,  AND  TR  ANS  RE  S  I  S'^  AN  CE  ,  RESPECTIVELY.  NOTE:  A  MORE  COMPLETE  DES¬ 

CRIPTION  OF  dependent  sources  as  implemented  in  SPICE  IS  GIVEN  IN  APPENDIX  B. 

***w  LINEAR  voltage-controlled  CURRENT  SOURCES 

GENERAL  FORM  GXXXXXXX  N+  N-  NC+  NC-  VALUE 

example  GI  2050  O.IMMHC 

N+  AND  N-  ARE  the  POSITIVE  AND  NEGATIVE  NODES,  RESPECTIVELY,  CURRENT  FLOW 
IS  FROM  THE  POSITIVE  NODE,  THROUGH  THE  SOURCE,  TO  THE  NEGATIVE  NODE.  NC+  AND 

NC-  ARE  the  POSITIVE  AND  NEGATIVE  CONTROLLING  NODES,  RESPECTIVELY,  VALUE  IS 

THE  TR  ANSCONDUC  TANCE  (IN  MHOS). 

****  LINEAR  voltage-controlled  VOLTAGE  SOURCES 
GENERAL  RORM  EXXXXXXX  N+  N-  NC+  NC-  VALUE 

example  El  2  3  14  1  2.0 

N+  Is  THE  POSITIVE  NODE,  AND  N-  IS  THE  NEGATIVE  NODE,  NC+  AND  NC-  ARE  THE 
POSITIVE  AND  NEGATIVE  CONTROLLING  NODES,  RESPECTIVELY,  VALUE  IS  THE  VOLTAGE 
GAIN. 

♦***  LINEAR  CURRENT-CONTROLLED  CURRENT  SOURCES 
GENERAL  form  FXXXXXXX  N+  N-  VNAM  VALUE 

example  FI  13  5  VSENS  5 

N+  AND  N-  ARE  THE  COSITIVE  AND  NEGATIVE  NODES,  RESPECTIVELY,  CURRENT  FLOW 
IS  FROM  Tuc  POSITIVE  NODE,  THROUGH  THE  SOURCE,  TO  THE  NEGATIVE  NODE.  VNAM  IS 
the  name  of  a  voltage  SOURCE  THROUGH  WHICH  thE  CONTROLLING  CURRENT  FLOWS.  THE 
DIRECTION  OF  OQSITIVE  CONTROLLING  CURRENT  FLOW  IS  FROM  THE  POSITIVE  NODE, 
through  the  source,  to  “HE  NEGATIVE  NODE  OF  VNAM,  VALUE  IS  THE  CURRENT  GAIN, 

4**4  LINEAR  current-controlled  VOLTAGE  SOURCES 

general  form  HXXXXXXX  N+  N-  VNAM  VALUE 

example  HX  5  17  V2  0  •  5K 

N+  AND  N-  ARE  THE  POSITIVE  AND  NEGATIVE  NODES,  RESPECTIVELY.  VNAM  IS  THE 
NAME  OF  A  voltage  SOURCE  THROUGH  WHICH  THE  CONTROLLING  CURRENT  FLOWS.  THE 
DIOECTION  OF  POSITIVE  CONTROLLING  CURRENT  FLOW  IS  FROM  THE  POSITIVE  NODE, 
through  the  SOURCE,  TO  THE  NEGATIVE  NODE  OF  VNAM,  VALUE  IS  THE  tra  nSP.E  SI  STANCE 
(  I  N  OHMS)  , 


C)  U  «/) 


****  IMOSPEN'OENT  sources 


GENEPAL  =0PM  VXXXXXXX  N-*-  N-  hoc]  DC/TRAN  VALUE]  [AC  [ACMAG  [ACPHASE]]] 

lYYYYYYY  N-t-  N-  [[DC]  OC/TSAN  VALUE]  [AC  [ACMAG  [  A  C  =  '-' A  SE  ]  ]  ] 

examples  VCC  i  0  0  DC  6 

VIN  13  2  0.001  AC  1  S:n(0  1  IMEG) 

ISPC  23  21  AC  0.333  A5.0  SFFM(0  1  lOK  5  IK) 

VMEAS  12  9 


N+  AND  N-  APE  the  positive  AND  NEGATIVE  NODE  S »  RESPECTIVELY.  NOTE  THAT 
voltage  SOURCES  NEED  NOT  BF  GROUNDED.  CURRENT  IS  ASSUMED  TO  FLOW  FQOM  THE  POS¬ 
IT  t  vh  NOOEt  THROUGH  THE  SOURCE,  TO  THE  NEGATIVE  NODE. 


OC/TRAN  IS  THE  DC  AND  TRANSIENT  ANALYSIS  VALUE  OF  THE  SOURCE.  IF  THE 
SOURCE  value  is  zero  both  pnp  DC  AND  TRANSIENT  ANALYSES,  THIS  VALUE  MAY  BE 
omitted.  TR  the  source  value  IS  TIME-INVARIANT  ( £. G. ,  A  POWER  SUP°LY),  THEN 
•^HE  value  may  OPTIONALLY  BE  PRECEDED  BY  THE  LETTERS  DC. 


ACMAG  IS  THE  AC  magnitude  AND 
this  value  in  THE  AC  ANALYSIS. 

,  A  value  of  unity  is  ASSUMED. 
SUMED.  IP  THE  SOURCE  IS  NOT  AN 
AC  values  APE  omitted. 


ACPHASE  IS  THE  AC  PHASE.  THE  SOURCE  IS  SET 
IF  ACMAG  IS  OMITTED  FOLLOWING  THE  KEYWORD 
IF  ACPHASE  IS  OMITTED,  A  VALUE  OF  ZERO  IS 
AC  small- SIGNAL  INPUT,  THE  KEYWORD  AC  AND  THE 


ANY  INDEPENDENT  SOURCE  CAN  BE  ASSIGNED  A  T IM E- DEPENDENT  VALUE  FOR 
transient  analysis.  if  a  source  is  ASSIGNED  A  TI mE-DEPENOENT  value,  the  TIME- 
ZERO  VALUE  IS  USED  PQR  DC  ANALYSIS.  THERE  ARE  FIVE  INDEPENDENT  SOURCE  FUNC¬ 
TIONS:  PULSE,  exponential,  SINUSOIDAL,  PIECE-WISE  LINEAR,  AND  SINGLE-FREQUENCY 

PM,  IP  parameters  other  than  source  values  are  OMITTED  OR  SET  TO  ZERO,  THE 
DERAUL"  VALUES  SHOWN  WILL  BE  ASSUMED,  (TSTEP  IS  THE  PRINTING  INCREMENT  AND 
TSTO°  IS  the  FINAL  time  (SEE  "^HE  .TRAN  CARD  FOR  EXPLANATION)), 

1,  PULSE  PULSE  CVl  V2  TD  TR  TF  Pw  P=R) 

Examole  VIN  3  0  PULSE(-1  1  2NS  2NS  2NS  SONS  lOONS) 

parameters  and  default  values  units 


V  1 

INITIAL  value 

volts  or  AMPS 

V  2 

PULSED  VALUE 

— 

VOLTS  OR  AMPS 

delay  time 

0.0 

SECONDS 

T  P 

PISE  time 

TSTEP 

SECONDS 

-F 

^ALL  *^1  ME 

TSTEP 

seconds 

PULSE  width 

TSTOP 

SECONDS 

DSP 

PE  PI  DO 

tstop 

SECONDS 

SINGL 


PULSE  SP  SPECIPTED  IS  DESCPIBED  BY  THE  FOLLOWING  TABLE 


T  IMiE  VALUE 


0 

VI 

^D 

VI 

TO+tr 

V2 

TDttq+dw 

V2 

TD-^tR^.d  w+TF 

VI 

TSTOP 

VI 

IN  TER  MED  I  ate 

POINTS  ARE  DET 

2#  SINUSOIDAL 

SIN(VO  VA  FREQ 

E  XAMPLE 

VIN  3  0  SI  N(  0 

epmined  by  linear  inteppolat  ion# 

TD  THETA ) 

1  I OOMEG  INS  lElO ) 


=>ARAMETEPS  and  default  VALUES 


UN  ITS 


VO  0=‘F5ET 

VA  A^^PLI^UOE 

='FEQ  «='REOUENCY 

•^D  DELAY 

THETA  DAMPING  FACTOR 


-  VOLTS  CR  AMPS 

-  VOLTS  OR  AMPS 

1/TSTOP  H2 

0#0  SECONDS 

0#0  I/SECONDS 


THE  SHAPE  OF  THE  WAVEFORM  IS  DESCRIBED  BY  THE  FOLLOWING  TABLE: 
TIME  VALUE 


0  TO  TD 
TD  TO  TSTO° 


VO 

VO  +  VAAEXP{ - (T IME-TD ) *THETA)  *SINE( TWOOI^FREQ^C  TI ME-TD)  ) 


3#  exponential  EXP(V1  V2  TDl  TAUl  TD 2  TAU2) 

EXAMOLE  VIN  3  0  EXP(-4  -1  2NS  30NS  60NS  40NS) 

PARAME*^ERS  AND  DEFAULT  VALUES  UNITS 


V  ! 

INITIAL  VALUE 

— 

volts  CR  AMPS 

V  2 

PULSE 

D  value 

— 

VOLTS  OR  AMPS 

^D1 

PISE 

DELAY  TIME 

0  #0 

SECONDS 

T  AUl 

R  ISE 

TIME  CONSTANT 

TSTEP 

SECONDS 

T  02 

FALL 

DELAY  T IME 

TD1-«-TSTEP  SECONDS 

T  AU2 

^  ALL 

TIME  constant 

TSTEP 

SECONDS 

THE  SHAPE  OF  THE  WAVEFORM  IS  DESCRIBED  BY  THE  FOLLOWING  TABLE: 
-^TME  VALUE 

VI 

VI  +  ( V2- VI  )* ( 1 -EXP (- (T IME-TDl  ) /T  AUX  >  ) 
Vl-f(V2-Vl)*<l-EXP(-<TI  ME- TDl  )  /TAUl  )  ) 

V1-V2  (1  -EXP  (  -  <T  IME-TC2  >/T  AU2  )  ) 


0  TO  ^D1 
TDl  TO  TD2 
TD2  *0  TSTOP 


A.  OIECS-WISE  LINEAR 

owlCTI  VX  its  V2  -^3  V3  TA  V4 

example  VCLOCK  7  5  PWL(0  -7  IONS  IIN  S  -3  XTNS  -3  1 8N  S  -7  SONS  -7) 

PARAMETERS  and  DEFAULT  VALUES 

EACH  PAIR  OF  VALUES  (Tl,  VI)  SPECIFIES  THAT  THE  VALUE  OF  THE  SOURCE  IS  VI 
(IN  VOLTS  OR  AMPS)  AT  TIMEsTI.  THE  VALUE  CF  THE  SOURCE  AT  INTERMEDIATE  VALUES 
OF  TIME  IS  DETERMINED  EY  USING  LINEAR  INTERPOLATION  ON  THE  INPUT  VALUES* 


5*  S I NGLE-PREOUENCY  FM 

SFFm(V0  VA  FC  mdI  FS  ) 

example  VI  12  0  SFRM( 0  IM  20K  S  IK) 

oarampters  and  default  values  units 

VO  offset  —  VOLTS  OP  AMPS 

VA  AMPLITUDE  VOLTS  OR  AMPS 

FC  CARRIER  frequency  1/TSTOP  HZ 

MDI  MODULATION  INDEX  -  — 

PS  SIGNAL  FREQUENCY  1/TSTOP  HZ 

THE  SHAPE  OF  the  WAVEFORM  IS  DESCRIBED  BY  THE  FOLLOWING  EQUATION; 

VALUE  S  VO  VA*SI  NE(  (TWOP  I*FC*T  IME  )  +  MO  I *S  I NE  ( T WOP  I  *FS *7  I  ME )  ) 


****  S=viicCNDUCTOR  DEVICES 


TMC  elements  that  have  been  DESCRIB'D  TC  THIS  POINT  TYPICALLY  REQUIRE  ONLY 
A  FEW  DARAMETER  VALUES  TO  SPECIFY  COMPLETELY  THE  ELECTRICAL  CHARACTERISTICS  OF 
the  element.  however,  the  models  for  the  FOUR  SEMICONDUCTOR  DEVICES  THAT  ARE 
INCLUDED  IN  THE  SPICE  PROGRAM  REQUIRE  MANY  PARAMETER  VALUES.  MOREOVER,  MANY 
DEVICES  IN  A  CIRCUIT  often  aRE  DEFINED  BY  THE  SAME  SET  OF  DEVICE  MODEL  PARA¬ 
METERS.  FOR  these  REASONS,  A  SET  OF  DEVICE  MODEL  PARAMETERS  IS  DEFINED  ON  A 
SEPARATE  .MODEL  CARO  AND  ASSIGNED  A  UNIQUE  MODEL  NAME,  THE  DEVICE  ELEMENT 
CARDS  IN  SPICE  THEN  REFERENCE  THE  MODEL  NAME.  THIS  SCHEME  ALLEVIATES  THE  NEED 

TO  Specify  all  of  the  model  parameters  cn  each  device  element  card. 

EACH  DEVICE  element  CARD  CONTAINS  THE  DEVICE  NAME,  THE  NODES  TO  WHICH  THE 
DEVICE  IS  CONNECTED,  AND  THE  DEVICE  MODEL  NAME.  IN  ADDITION,  TWO  OPTIONAL 
=ARAMETERS  MAY  BE  SPECI'IED  FOR  EACH  DEVICE;  AN  AREA  FACTOR,  AND  AN  INITIAL 
CONDITION. 

THE  AREA  factor  DETERMINES  THE  NUMBER  OF  EQUIVALENT  PARALLEL  DEVICES  OF  A 
SPECI'IED  MODEL.  THE  AFFECTED  PARAMETERS  ARE  MARKED  WITH  AN  ASTERISK  UNDER  THE 
HEADING  FAREA*  IN  THE  MODEL  DESCRIOTICNS  BELOW. 

TWO  DIFFERENT  FORMS  O'  INITIAL  CONDITIONS  MAY  BE  SPECIFIED  FOR  DEVICES. 

THE  'IRST  FORM  IS  INCLUDED  TO  IMPROVE  THE  DC  CONVERGENCE  FOR  CIRCUITS  THAT 
CONTAIN  MORE  THAN  ONE  S'^ABLE  STATE.  IF  A  DEVICE  IS  SPECIFIED  OFF,  THE  DC 
operating  point  is  DETERMINED  WITH  THE  TERMINAL  VOLTAGES  FOR  THAT  DEVICE  SET  TO 
ZERO.  after  convergence  IS  OBTAINED,  THE  PROGRAM  CONTINUES  TO  ITERATE  TO 

Obtain  the  exact  value  for  the  terminal  voltages,  if  a  circuit  has  more  than 

ONE  DC  STABLE  STATE,  THE  OFF  OPTION  CAN  BE  USED  TO  FORCE  THE  SOLUTION  TO 
CORRESPOND  TO  A  DESIRED  STATE.  IF  A  DEVICE  IS  SPECIFIED  OFF  WHEN  IN  REALITY 
the  device  is  conducting,  THE  PROGRAM  WILL  STILL  OBTAIN  THE  CORRECT  SOLUTION 
(ASSUMING  THE  SOLUTIONS  CONVERGE)  BUT  MORE  ITERATIONS  WILL  BE  REQUIRED  SINCE 
the  program  must  independently  converge  to  two  separate  SOLUTIONS. 

the  second  form  of  INITIAL  CONDITION  SPECIFICATION  (USING  IC*...)  IS 
PROVIDED  ~0  ALLOW  THE  USER  TO  BYOASS  THE  DC  OPERATING  POINT  CALCULATION 

normally  made  before  the  start  of  transient  analysis.  these  initial  CONDI¬ 
TIONS  ARE  USED  BY  SPICE  *ONLY*  IF  thE  UIC  OPTION  IS  GIVEN  ON  THE  .TRAN  CARD. 


****  JUNCTION  DIODES 

general  form  DXXXXXXX  NT  N-  MNAME  [APEaI  [OFF]  tIC=VDl 

example  DERIDGE  2  10  DIODE! 

DCLMP  3  7  DMOD  3.0  1C=0.2 

N+  AND  N-  ARE  THE  POSITIVE  AND  NEGATIVE  NODES,  RESPECTIVELY.  MNAME  is  THE 
model  name,  area  IS  THE  AREA  FACTOR,  AND  OFF  INDICATES  AN  (OPTIONAL)  INITIAL 
CONDITION  ON  THE  DEVICE  FOR  DC  ANALYSIS.  IF  THE  AREA  FACTOR  IS  OMITTED,  A 
VALUE  OF  1,0  IS  assumed,  THE  (OPTIONAL)  INITIAL  CONDITION  SPECIFICATION  USING 
1C  =  VD  IS  in-tended  fop  use  with  the  UIC  OPTION  ON  THE  •’'RAN  CARO,  WHEN  A  TRAN¬ 
SIENT  ANALYSIS  IS  DESIRED  STAR'^ING  FROM  OTHER  THAN  THE  QUIESCENT  OPERATING 

POINT, 


****  btdolAR  junction  transistors 


GENERAL  *^ORM  OXXXXXXX  NC  NB  NE  MNAME  (AREA]  (  OPF  ]  IICsVBE,VCEl 

example  022  10  2A  13  OMOD  IC=0«6,5*0 

NC  ,  N8,  ANO  ne  are  "^HE  collector,  ease,  and  emitter  nodes,  respectively. 

MNAME  IS  THE  MODEL  NAME,  AREA  IS  THE  AREA  FACTOR,  AND  OFF  INDICATES  AN 
(OOTICNAL)  initial  CONDITION  ON  THE  DEVICE  FOR  THE  DC  ANALYSIS.  IF  THE  AREA 
factor  is  omitted,  a  value  OF  1,0  IS  ASSUMED,  THE  (OPTIONAL)  INITIAL  CONDITION 
SPECIFICATION  USING  IC=VPE,VCE  IS  INTF^oED  FOR  USE  WITH  THE  UIC  OPTION  ON  THE 
.'•RAN  CARO,  WHEN  A  TRANSIENT  ANALYSIS  IS  DESIRED  STARTING  FROM  OTHER  THAN  THE 
DUIESCENT  0®ERATING  point. 


****  JUNCTION  field-effect  transistors 

general  =GPM  JXXXXXXX  ND  NG  ns  MNAME  [AREA]  [OFF]  [IC=VCS,VGSl 

EXAmoi_e  J1  7  2  3  JMl  OFF 

ND ,  NG,  AND  NS  ARE  THE  DRAIN,  GATE,  AND  SOURCE  NODES,  RESPECTIVELY,  MNAME 
IS  THE  model  name,  area  is  THE  AREA  FACTOR,  AND  OFF  INDICATES  AN  (OPTIONAL) 
initial  CONDITION  ON  THE  DEVICE  FOR  DC  ANALYSIS.  IF  THE  AREA  FACTOR  IS  OMIT¬ 
TED,  A  value  of  1,0  IS  ASSUMED.  THE  (OPTIONAL)  INITIAL  CONDITION  SPECIFICATION 
USING  IC*VDS,VGS  IS  INTENDED  FOR  USE  WITH  THE  UIC  OPTION  ON  THE  .TRAN  CARD, 
when  a  transient  ANALYSIS  IS  DESIRED  STARTING  FROM  OTHER  THAN  THE  QUIESCENT 
operating  POINT, 


*F*F  MOSFETS 

GENERAL  fqdm  MXXXXXXX  NC  NG  NS  NB  MNAME  [W=VALl  [L=VAL]  (A0=VAL]  [AS=VAL] 

+  [OFF]  I IC=VCS, VGS, VBS 3 

EXAMOlE  M31  2  17  6  10  mqdm  L=2VIL  W=0.5MIL 

ND,  NG,  NS,  AND  NB  ARE  THE  DRAIN,  GATE,  SOURCE,  AND  BULK  (SUBSTRATE) 

NODES,  RESPFC"! VELY,  '  MNAME  is  THE  MODEL  NAME.  W  AND  L  ARE  THE  CHANNEL  WIDTH 
AND  length,  in  centimeters;  IF  OMITTED,  BOTH  W  AND  L  ARE  ASSUMEQ  tO  BE  1,0, 

AD  AND  AS  APS  THE  AREAS  OF  THE  DRAIN  AND  SOURCE  DIFFUSIONS,  IN  CMt2;  IF  NOT 
SPECIFIED,  50'H  AREAS  ARE  ASSUMED  TO  BE  l.OE-6.  OFF  INDICATES  AN  (OPTIONAL) 
INITIAL  CONDITION  ON  THE  DEVICE  FOR  DC  ANALYSIS,  THE  (OOTIONAL)  INITIAL  CONDI¬ 
TION  S=>EC  I^IC  AT  ION  USING  I  C  =  VDS ,  VGS  ,  VBS  IS  INTCmoeo  pqR  USE  WITH  the  UIC  OPTION 
DN  THE  .TRAN  CARD,  WHEN  A  TRANSIENT  ANALYSIS  IS  DESIRED  STARTING  FROM  OTHER 
than  the  quiescent  OPERATING  POINT, 


»  lii 


^MODEL  CARD 

GENERAL  «='nRV  •VCDEL  MNAME  TY  PE  ( PN  A  WEI  =  PV  AL 1  PNAME2=PVAL2  •••  ) 

example  .model  MCDl  NPN(BFs50  IS=1E-IE  VA^SO) 

THE  .MODEL  CARD  SPECIFIES  A  SET  OF  MODEL  PARAME"^ERS  THAT  WILL  BE  USED  BY 
ONE  OR  MORE  DEVICES.  MNAM^  IS  THE  MODEL  NAMEf  AND  TYPE  IS  ONE  OF  THE  FOLLOWING 
SEVEN  TYPES: 

NPN  NPN  BJT  MODEL 

ONP  PNP  BJT  MODEL 

D  DIODE  MODEL 

NJF  N-CHANNEL  JFET  MODEL 

oj^  P-CHANNEL  JFET  MODEL 

NMOS  N-CHANNEL  MOSFET  MODEL 

PMOS  D-CHANNEL  MOSFET  MODEL 

parameter  values  are  DE'=^INED  by  appending  the  parameter  name,  as  GIVEN 

BELOW  FOR  EACH  MODEL  TYPE,  FOLLOWED  BY  AN  EQUAL  SIGN  AND  THE  PARAMETER  VALUE. 
MODEL  PARAMETERS  THAT  APE  NOT  GIVEN  A  VALUE  ARE  ASSIGNED  THE  DEFAULT  VALUES 
GIVEN  BELOW  FOR  EACH  MODEL  TYPE. 


-  DIODE  MODEL 

THE  DC  CHARACTERISTICS  OF  THE  DIODE  ARE  DETERMINED  BY  THE  PARAMETERS  IS 
AND  N.  AN  OHMIC  RESISTANCE,  RS,  IS  INCLUDED.  CHARGE  STORAGE  EFFECTS  ARE 
MODELED  BY  A  TRANSIT  TIME,  TT ,  AND  A  NONLINEAR  DEPLETION  LAYER  CAPACITANCE 
MICH  IS  determined  by  the  PARAMETERS  CJC,  PB,  AND  M.  THE  TEMPERATURE  DEPEND- 
NCE  OF  THE  SATURATION  CURRENT  !$  DEFINED  BY  THE  PARAMETERS  EG,  THE  ENERGY  GAR, 
AND  PT,  THE  SATURATION  CURRENT  TEMPERATURE  EXPONENT.  REVERSE  BREAKDOWN  IS 
MODELED  BY  AN  EXPONENTIAL  INCREASE  TN  THE  REVERSE  DIODE  CURRENT  AND  IS  DETER¬ 
MINED  BY  the  parameters  BV  AND  IPV  (BOTH  QR  WHICH  APE  POSITIVE  NUMBERS). 


AREA  NAME 


parameter 


default  typical 


F  IS 
F  RS 
N 

tt 

F  c  JO 
PB 
M 

EG 


9 

DT 

1  0 

KF 

1  1 

A=^ 

1  2 

FC 

1  3 

B  V 

1  4 

!BV 

1 

2 

•7 

4 

s 

6 

7 

Q 


SATURATION  CUPPEN-^  l.OE-14 

OHMIC  RESISTANCE  0 

EMISSION  COE^«=ICIENT  1 

transit  time  0 

ZERO-BIAS  junction  CAPACITANCE  0 

JUNCTION  POTENTIAL  1 

GRADING  COEFFICIENT  0.5 

ENERGY  GAP  l.ll 

SATURATION  CURRENT  TEMP.  EXPONENT  3.0 

FLICKER  NOISE  COEFFICIENT  0 

FLICKER  NOISE  EXPONENT  1 

CRWAPD-BI  AS  NON  IDEAL  JUNCTION 
CAPACITANCE  COEFFICIENT  0.5 

REVERSE  BOEAKDOWN  VOLTAGE  100.0 

CURRENT  AT  breakdown  VOLTAGE  l.OE-3 


I  «  OE-14 
1  0 
1.0 

0.  1  NS 
2op 
0.  6 
0.  5 

1.11  51 

0.69  SBD 
0.67  GE 
3.  0  JN 
2.0  SBD 


40.  0 


- BJT  MCC£LS  (BCTH  NIPN  AND  PNP) 

the  bipolar  junction  TRANSTSTOO  model  in  spice  is  an  adaptation  0^  THE 
INTEGRAL  CHARGE  CCNTPOL  MODEL  OF  GUMMEL'  AND  PCON;  HOwEVE^,  IF  THE  ADDITIONAL 
GUMMEL-PCON  =ARAMETEPS  are  not  specified,  the  simpler  EBEPS-MOLL  model  of 
SOTCE,  version  1,  IS  USED.  THE  DC  MODEL  IS  DEFINED  BY  THE  PARAMETERS  C2f 

IK,  AND  NE  ,  WHICH  DETERMINE  THE  FORWARD  CURRENT  GAIN  CHARACTERISTICS,  BR,  C4, 
TKP,  AND  NC,  WHICH  DETERMINE  THE  REVERSE  CURRENT  GAIN  CHARACTERISTICS,  VA  AND 
VB,  WHICH  determine  the  output  CONDUCTANCE  FOR  FORWARD  AND  REVERSE  REGIONS,  AND 
the  SATucaTION  current,  is.  THREE  OHMIC  RESISTANCES,  RB,  RC ,  AND  RE,  ARE 
INCLUDED.  RASE  CHARGE  STORAGE  IS  MODELED  BY  FORWARD  AND  REVERSE  TRANSIT  TIMES, 
TF  AND  TP,  AND  NONLINEAR  DEPLETION  LAYER  CAPACITANCES  WHICH  ARE  DETERMINED  BY 
CUE,  PE,  AND  ME  POR  THE  B- E  JUNCTION,  AND  CJC,  PC,  AND  MC  FOP  •^HE  B-C  JUNCTION. 
A  CONSTANT  COLLECTOO-SUBSTP ate  CAPACITANCE,  CCS,  IS  ALSO  INCLUDED.  THE 
temdeRATURE  DEPENDENCE  OF  THF  SATURATION  CURRENT  IS  DETERMINED  BY  THE  ENERGY 
GAP,  EG,  AND  ^HE  SATURATION  CURRENT  TEMPERATURE  EXPONENT,  PT • 

AREA  NAME  PARAMETER  DE^^AULT  TYPICAL 


1 

BF 

2 

BR 

3 

a 

IS 

4 

3k 

RB 

5 

3k 

R  C 

6 

A 

RE 

7 

V  A 

S 

VB 

<5 

I  K 

1  0 

C2 

1  1 

NE 

1  2 

I  KP 

1  3 

CA 

1  A 

NC 

1  5 

TF 

I  A 

TP 

1  7 

3k 

CCS 

1  9 

t 

CUE 

1  9 

DP 

20 

ME 

2  1 

F 

CJC 

22 

PC 

23 

MC 

2  A 

EG 

2  5 

PT 

2  6 

KF 

2T 

AF 

2  8 

IDEAL  FOPWAPO  CUPRENT  GAIN 
IDEAL  REVERSE  CURRENT  GAIN 
SATURATION  CURRENT 
BASE  OHMIC  RESISTANCE 
COLLECTOR  OHMIC  RESISTANCE 
EMITTER  OHMIC  RESISTANCE 
^CRWAPD  EARLY  VOLTAGE 
REVERSE  EARLY  VOLTAGE 
FORWARD  HIGH-CURPENT  KNEE  CURRENT 
FORWARD  LOW-CURRENT  NONIDEAL 

BASE  CURRENT  COEFFICIENT 

ncnideal  low-curpent  base-emitter 

EMISSION  COEFFICIENT 
REVERSE  HIGH-CURRENT  KNEE  CURRENT 
REVERSE  LOW-CURRENT  NCNIDEAL 

BASE  CURRENT  COEPFiciENT 
NONIDEAL  LOW-CURRENT  BASE-COLLECTOR 

EMISSION  COEFFICIENT 
FOPWAPO  TRANSIT  TIME 
REVERSE  TRANSIT  time 
COLLECTCR-SUeSTP ATE  CAPACITANCE 
ZERC-®IAS  B-E  JUNCTION  CADACITANCE 
B-E  JUNCTION  POTEN^^TAL 
9-E  JUNCTION  GRACING  COEFFICIENT 
ZERO  BIAS  B-C  JUNCTION  CAPACITANCE 
B-C  JUNCTION  potential 
B-C  GRADING  COEFFICIENT 
ENERGY  GAP 


100 

1 

1  .  OF-14 
0 
0 
0 

INFINITE 
I  NF  I  N  I  TE 
INFIN  ITE 

0 

2.  0 

INFIN  ITE 
0 

2.  0 
0 
0 
0 
0 

1  .  0 

0.5 

0 

1.0 

0.5 

1.11 


SATURATION  current  TEMP.  EXPONENT  3.0 

flicker  noise  coefficient  0 


flicker  noise  exponent 


1 


=^0RWAR0-3I  AS  ncnideal  JUNCTION 
capacitance  COEFFIC I ENT 


100 

0.1 

1 . OE-16 
100 
1  0 

1 

200 

200 

lOMA 

1000 

2.  0 
1  OOMA 

1.0 

2b  0 

O.INS 
1  ONS 
2PF 
2PF 
0.  T 

0.33 
IPF 
0.  5 
0.33 

1.11  31 

0.67  GE 

6.6S-16  NPN 
6.3E-13  PNP 
1.0  NPN 
1.5  ONP 


0.5 


jpp-r  mCOELS  (BOTH  N  AND  p  CHANNEL) 


THE  JF5T  MODEL  IS  DERIVED  FPpM  thE  FET  MODEL  OF  SHICHMAN  AND  HDDGES.  THE 
DC  CHARACTEPI  STI C5  AOF  DEFINED  0Y  THE  PARAMETERS  VTD  AND  BETAt  WHICH  DETERMINE 
the  variation  of  DRAIN  CURRENT  WITH  GATE  VOLTAGE»  LAMBDA,  WHICH  DETERMINES  THE 
OUTPUT  CONDUCTANCE,  AND  IS,  THE  SATURATION  CURRENT  OF  THE  TWO  GATE  JUNCTIONS. 
TWO  OHMIC  RESISTANCES,  RD  AND  RS,  ARE  INCLUDED#  CHARGE  STORAGE  IS  MODELED  BY 
NONLINEAR  DEPLETION  LAYER  CAPACITANCES  FOR  BOTH  GATE  JUNCTIONS  WHICH  VARY  AS 
the  -1/2  POWER  OF  JUNCTION  VOLTAGE  AND  ARE  DEFINED  BY  THE  PARAMETERS  CGS,  CGO, 
AND  P8# 


AREA  NAME 


PARAMETER 


DEFAULT  TYPICAL 


1 

2  * 

3 

4  « 

5  * 

6  * 

7  * 

6 

o  A 

1  0 
1  1 
I  2 


VTO 

BE^A 

L  AV°D 

PD 

RS 

CGS 

CGD 

PB 

I  S 

KF 

AF 

FC 


threshold  VQL'^AGE  -2#0 

TPANSCONDUCTANCE  PARAMETER  l#0E-4 

CHANNEL  length  MODULATION  PARAMETER  0 

DRAIN  OHMIC  RESISTANCE  0 

SOURCE  OHMIC  RESISTANCE  0 

ZERO-EIAS  G-S  JUNCTION  CAPACITANCE  0 

ZERO-EIAS  G-D  junction  CAPACITANCE  0 

gate  JUNCTION  POTENTIAL  1 

gate  JUNCTION  SATURATION  CURRENT  X#OE-14 

FLICKER  NOISE  COEFFICIENT  0 

FLICKER  NOISE  EXPONENT  1 

FCRWARD-eiAS  NONIDEAL  JUNCTION 
CADACI'^ANCE  COEFFICIENT 


-2#  0 
1 #05-3 
1 • OE-4 
100 
100 
5PF 
IPF 
0#  6 


1 .OE-IA 


0# 


-  MQSrET  MODELS  { 50TH  N  AND  P  CHANNEL) 


THE  MOS-ET  MODEL  IS  DEPIVED  PROM  THE  F B OHM AN -G PO VE  MODEL.  THE  DC  CHARAC- 
TS®ISTir.S  OP  the  MOSPET  ARE  DEFINED  SY  THE  PARAMETERS  VTO  t  SETA,  LAMBDA,  PHI, 
AND  gamma,  VTO  IS  POSITIVE  (NEGATIVE)  FOR  ENHANCEMENT  MODE  AND  NEGATIVE  (POSI¬ 
TIVE)  POP  DEPLETION  MODE  N-CHANNEL  (P-CHANNEL)  DEVICES.  CHARGE  STORAGE  IS 
MODELED  BY  three  CONSTANT  CAPACITORS,  CGS,  CGD,  AND  CGB,  BY  THE  NONLINEAR  OXIDE 
capacitance  WHICH  IS  DISTRIBUTED  AMONG  THE  GATE-SOURCE,  GATE-CRAIN,  AND  GATE- 
SULK  REGIONS  USING  t hE  PORMULATION  OF  J.  E.  MEYER,  AND  BY  THE  NONLINEAR  DEPLE¬ 
TION-LAYER  CAPACITANCES  FOP  SOTH  SUBSTRATE  JUNCTIONS  WHICH  VARY  AS  THE  -1/2 
POWER  OF  JUNCTION  VOLTAGE  AND  ARE  DETERMINED  BY  THE  PARAMETERS  C3D,  CBS,  AND 
RS. 


NAME 

OARAMETE  p 

default 

"TYPICAL 

UNITS 

1 

VT  0 

ZEPO-eiAS  THPESHCLD  VOLTAGE 

0«  0 

-0#1 

V 

2 

KP 

INTRINSIC  TP.ANSCONDUCT  ANCE  PARAMETEP 

1.  OE-5 

3*  IE- 5 

A/V  ^2 

3 

GA  MMA 

SULK  threshold  PARAMETER 

0^  0 

0*37 

V  ♦  (  1  /  2  ) 

4 

OH  I 

sup^ace  potential  at  strong  inversion 

0.6 

0«65 

V 

c 

LAMBD  A 

CHANNEL- LENGTH  MOO  U  LA I  CN  ©ARAMETEP 

0«0 

0*02 

/V 

6 

DPAIN  OHMIC  RESISTANCE 

0.  0 

1.0 

OHMS 

T 

PS 

SOURCE  ohmic  resistance 

0«0 

1  *0 

OHMS 

s 

CGS 

GATE-SOUPCE  OVEPLAO  CAPACITANCE 

©ER  CM  CHANNEL  WIDTH 

0.  0 

4. 06-  13 

F/CM 

9 

CGD 

GATE-DRAIN  OVERLAP  CAPACITANCE 

PER  CM  CHANNEL  WIDTH 

0.  0 

4*  OE-13 

F/CM 

1  0 

CGB 

gate-bulk  overlap  caoacitance 

OER  CM  CH/-NNEL  length 

0.  0 

2. OE-12 

F/CM 

1  1 

CBD 

ZERO-BIAS  S-D  junction  CAPACITANCE 

PEP  CM^2  OF  JUNCTION  AREA 

0.0 

2tOE-8 

F/CM^2 

1  2 

CSS 

ZERO-BIAS  B-S  JUNCTION  CAPACITANCE 

PEP  CM^2  OF  JUNCTION  AREA 

0.0 

2m OE- 8 

F/CM*  2 

1  3 

TDX 

OXIDE  thickness 

INFINITY 

1  #06-5 

CM 

IF  nsub  specified: 

U  OE-5 

1  4 

P9 

SULK  junction  potential 

0«  6 

0  *37 

V 

1  5 

JS 

SULK  junction  reverse  saturation  current 

PEP  CM^2  OF  junction  ADEA 

1#0f-3 

!• OE- 8 

A/CM*  2 

1  e 

NSU8 

EP=^ECTIVE  SUBSTRATE  DOPING 

Om  0 

OEl  5 

/CM*3 

1  T 

NS  S 

effective  surface  sta^e  density 

0.  0 

It  OEl  0 

/CM  *2 

1  3 

NF  5 

EFFECTIVE  FAST  SURFACE  STATE  DENSITY 

Oa  0 

1 tOEl 0 

/CM  *2 

1  9 

X  J 

METALURGICAL  junction  depth 

0 

1 • OE-4 

CM 

2  0 

LD 

LATERAL  DIFFUSION  COEFFICIENT 

0«  8 

0«  3 

21 

NG  ATE 

PQLYSILICON  GATE  DOPING 

AL  GATE 

1  •0E2  0 

/CMf2 

22 

TDS 

TYOE  of  POLYSILICON:  <-l)  FOR 

OPPOSITE  (SAME)  AS  substrate 

-*“1 .0 

23 

UO 

SUPPACE  MOBILITY 

700 

600 

CM*  2/V- 

24 

UCPI^ 

critical  field  for  mobility  degradation 

1.0E4.4 

1 . OE^  4 

V/CM 

2  5 

US  xo 

critical  field  exponent  (mobility) 

o«  o 

0«1 

26 

UT^  A 

transverse  field  coefficient  (MOBILITY) 

0.  0 

0.  3 

27 

KF 

FLICKER  noise  coefficient 

0*0 

2S 

AF 

flicker  noise  exponent 

It  0 

2  Q 

^C 

forward-eias  nonideal  JUNC“^I0N 

cap  act  tan CE  COEFfic lENT 

0»  5 

****  SUBCIPCUITS 


A  SUSCIPCUIT  tha-  consists  OF  5° ICE  ELEMENTS  CAN  BE  DEFINED  AND  REFERENCED 

IN  A  fashion  similar  tc  device  models*  the  sueciRcuiT  IS  defined  in  the  input 

DECK  FY  A  GROUPING  OF  ELEMENT  CARDS;  THE  PROGRAM  THEN  AUTOMATICALLY  INSERTS 
THE  GROUP  O'"  ELEMENTS  WHEREVER  THE  SUBCIRCUIT  IS  REFERENCED.  THERE  IS  NO  LIMIT 
ON  THE  SIZE  OR  COMPLEXITY  OF  SUBCIRCUITSt  AND  SUBCIPCUITS  MAY  CONTAIN  OTHER 
SUBCIRCUITS.  AN  example  OF  SUBCIRCUIT  USAGE  IS  GIVEN  IN  APPENDIX  A. 


****  .SUBCKT  CARD 

GENERAL  FORM  .SUBCKT  SUBNAM  N1  [  N2  N3  ***] 

example  .SUBCKT  OPAMP  1234 

A  SUPCIRCUIT  definition  IS  BEGUN  WITH  A  .SUBCKT  CARD.  SUBNAM  IS  THE  SUB- 
CIRCUIT  NAME,  AND  Nl,  N2,  o.*  ARE  THE  EXTERNAL  NODES,  WHICH  CANNOT  SE  ZERO. 

THE  GROUP  OF  ELEMENT  CARDS  WHICH  IMMEDIATELY  FOLLOW  THE  .SUBCKT  CARD  DEFINE  THE 
SUBCIRCUIT,  THE  LAST  CARD  IN  A  SUBCIRCUIT  DE=IN1TI0N  IS  THE  .ENDS  CARO  (SEE 
BELOW).  CtONT  POL  CARDS  MAY  NOT  APPEAR  WITHIN  A  SUBCIRCUIT  DEFINITION;  HOWEVER, 
SUBCIRCUIT  definitions  MAY  CONTAIN  ANYTHING  ELSE,  INCLUDING  OTHER  SUBCIRCUIT 
DEF-INI-^IONS,  DEVICE  MODELS,  AND  SUBCIRCUIT  CALLS  (SEE  BELOW).  NOTE  THAT  ANY 
DEVICE  MODELS  OR  SUECIRCUIT  DERINITIONS  INCLUDED  AS  PART  OF  A  SUBCIRCUIT  DEFI¬ 
NITION  APE  STRICTLY  LOCAL  (I.E.,  SUCH  MODELS  AND  DEFINITIONS  ARE  NOT  KNOWN 
OUTSIDE  the  SUBCIRCUIT  DEFINITION).  ALSO,  ANY  ELEMENT  NODES  NOT  INCLUDED  ON 
THE  .SUBCKT  CARD  ARE  STRICTLY  LOCAL,  WITH  thE  EXCEPTION  OF  0  (GROUND)  WHICH  IS 
ALWAYS  GLOBAL. 


*•*'**  .ENDS  CARD 

GENERAL  FORM  .ENOS  [SUBNAM) 

EXAMPLE  .ENOS  OPAMP 

T^IS  CARO  MUST  BE  THE  LAST  ONE  FOR  ANY  SUBCIRCUIT  DEFINITION,  THE  SUB- 
CTRC'JI-^  NAME,  I'^  INCLUDED,  INDICATES  WHICH  SUBCIRCUIT  DEFINITION  IS  BEING 

terminated;  if  omitted,  all  SUPCIRCUITS  BEING  defined  are  terminated,  THE 

NAME  IS  NEEDED  ONLY  WHEN  NESTED  SUBCIRCUlT  DEFINITIONS  ARE  BEING  MADE. 


SUECIRCUIT  CALLS 

GENERAL  FORM  XYYYYYYY  Nl  [ N2  N3  ***]  SUBNAM 

example  XI  2  4  17  3  1  MULTI 

SUBCIRCUITS  ARE  USED  IN  S»ICE  BY  SPECIFYING  =S EUOO-ELE MENT S  BEGINNING  WITH 
■HE  letter  X,  RDLLOWED  BY  THE  CIRCUIT  NODES  *0  BE  USED  IN  EXPANDING  THE  SUB- 
rioCUIT,  FOLLOWED  BY  THE  SUBCIRCUIT  NAME,  THE  NODES  MUST  SE  IN  thE  ORDER  THAT 
they  are  defined  in  CN  the  .SUBCKT  CARD. 


CON'^^01.  CARDS 


****  TT  H-E  CARO 

examrle  =>owep  amplifier  circuit 

this  card  must  EE  THE  FIRST  CARO  IN  THE  INPUT  DECK*  ITS  CONTENTS  ARE 
PRINTED  VERBATIM  AS  THE  HEADING  FOR  EACH  SECTION  OF  OUTPUT, 

****  .END  CARD 
EXAMo^_h  .end 

this  card  must  always  be  the  last  card  in  the  input  deck.  NOTE  THAT  "^HE 
RERIOD  is  an  integral  part  CF  the  name. 

CDMHRny  card 

general  RORM  a  ANY  COMMENTS 

example  *  RFslK  GAIN  SHOULD  BE  100 

THIS  CARD  IS  PRINTED  OUT  IN  THE  INPUT  LISTING  BUT  IS  OTHERWISE  IGNORED. 

.TEMB  CARD 

general  form  .temp  ti  Its  [te  ...Ij 

example  .TEMP  -55.0  25,0  125.0 

THIS  CARD  SPECIFIES  THE  TEMPERATURES  AT  WHICH  THE  CIRCUIT  IS  TO  BE  SIMULA¬ 
TED.  Ti,  T2,  ...  ARE  THE  DIFRERENT  TEMPERATURES,  IN  DEGREES  C.  TEMPERATURES 

LESS  THAN  -223.0  DEG  C  APE  IGNOPED.  MODEL  DATA  IS  SPECIFIED  AT  TNO M  DEGREES 
(SEE  TME  .OPTION  CARD  FOR  TNOM  )  ;  IF  THE  .TEMP  CARD  IS  OMITTED,  THE  SIMULATION 
ALSO  WILL  BE  PERFORMED  AT  A  TEMPERATURE  OF  TNOM. 

****  .width  card 

GENERAL  FOPM  .width  IN=C0LNUM 

example  .width  in  =  72 

COLNUM  is  the  last  column  read  from  each  line  0=  input;  THE  se~ting  takes 
EFFC<;7  with  the  nEX”  line  READ.  the  DEFAULT  VALUE  f^S  COLNUM  IS  60. 


****  .OOTIONS  CARD 

GENERAL  RORM  . OPTIONS  OFTl  0OT2  •••  (OR  CPT=OPTVAL  •••) 

example  .options  acct  list  node 

THIS  CARD  ALLOWS  THE  USER  TO  RESET  PROGRAM  CONTROL  AND  USER  OPTIONS  FOR 
SPECIFIC  SIMULATION  PURPOSES.  ANY  COMBINATION  OF  THE  FOLLOWING  OPTIONS  MAY  BE 
INCLUDED,  IN  ANY  ORDER.  *X*  (BELOW)  REPRESENTS  SOME  POSITIVE  NUMBER. 

OPTION  EFFECT 


ACCT  CAUSES  THE  EXECUTION  TIME  FOR  THE  VARIOUS  SECTIONS  OF 

the  program  to  BE  PRINTED  AS  WELL  AS  OTHER  ACCOUNTING 
INFORMATION  TO  BE  PRINTED. 

LIST  CAUSES  THE  SUMMARY  LISTING  OF  ThE  INPUT  DATA  TO  BE 

printed. 

NOMOD  SUPPRESSES  THE  PRINTOUT  OF  THE  MODEL  PARAMETERS. 

NOPAGE  SUPPRESSES  PAGE  EJECTS 

NODE  CAUSES  THE  NODE  TABLE  TO  EE  PRINTED. 

OPTS  CAUSES  THE  OPTION  VALUES  TO  BE  PRINTED. 

GM1N=X  RESETS  THE  VALUE  OF  GMIN,  THE  MINIMUM  CONDUCTANCE 

ALLOWED  BY  THE  PROGRAM.  THE  DEFAULT  VALUE  IS  l.OE-12. 

RELTOL=X  RESETS  THE  RELATIVE  ERROR  TOLERANCE  OF  THE  PROGRAM,  THE 

default  value  is  0.001  (0.1  PERCENT). 

A8ST0L=X  RESETS  THE  ABSOLUTE  CURRENT  ERROR  TOLERANCE  OF  THE 

PROGRAM.  THE  DEFAULT  VALUE  IS  1  PICOAMP. 

VNTOL=X  RESETS  THE  ABSOLLTE  VOLTAGE  ERROR  TOLERANCE  OF  THE 

PROGRAM.  THE  default  VALUE  IS  1  MICROVOLT. 

TRTOL=X  RESETS  the  TRANSIENT  ERROR  TOLERANCE.  THE  DEFAULT  VALUE 

IS  "^.O.  THIS  parameter  IS  AN  ESTIMATE  OF  THE  FACTOR  BY 
WHICH  SPICE  OVERESTIMATES  THE  ACTUAL  TRUNCATION  ERROR. 

CHGTOL=X  RESETS  THE  CHARGE  TOLERANCE  OF  ThE  PPOGRAM,  THE  DEFAULT 

VALUE  IS  I.OE-IA. 

NUMDGTsX  RESETS  THE  NUMBER  OF  SIGNIFICANT  DIGITS  PRINTED  FOR 

OUTPUT  variable  values.  X  MUST  SATISFY  THE  RELATION 

0  <  X  <  e.  the  default  value  is  a.  note:  this  ortion  is 

INDEPENDENT  OF  THE  ERROR  TOLERANCE  USED  BY  SPICE  (I.E.,  IF 
THE  VALUES  OF  CP^TICNS  RELTQL,  ABSTCL,  Etc.  ARE  NOT  CHANGED 
THEN  ONE  may  be  PRINTING  NUMERICAL  FNOISE*  FOR  NUMDGT  >  A, 

TNOM=X  RESETS  THE  NOMINAL  TEMPERATURE.  THE  DEFAULT  VALUE  IS 

27  DEG  C  (300  DEG  X  ) . 


Itl1=  X 
! TL2=X 

I  TL"?*X 

I TL4»X 

Itl5=x 

LIMTI M 

LI  mdt  S' 

L VLC3D 

LVL”I  M: 

•method 


RESETS  THS  DC  ITERATION  LIMIT,  THE  DEFAULT  IS  100, 

RESETS  THE  DC  TOanSFER  CURVE  ITERATION  LIMIT,  tmE 
DEFAULT  IS  20, 

-RESETS  "^he  lower  transient  analysis  iteration  limit, 
the  default  value  is  4, 

resets  the  transient  analysis  timepoint  iteration  limit, 
THE  default  is  10, 

RESETS  THE  TRANSIENT  ANALYSIS  TOTAL  ITERATION  LIMIT, 

THE  OEi=AULT  IS  5000, 

X  RESETS  the  AMOUNT  OF  TIME  RESERVED  BY  SPICE  '=0R 

generating  olOTS,  -THE  OES^AULT  VALUE  IS  2  (SECONDS), 

X  RESETS  the  total  NUMBER  0^  »OINTS  THAT  CAN  BE  PRINTED 

OR  PLOTTED  IN  A  DC,  AC,  OR  TRANSIENT  ANALYSIS,  THE 
default  value  is  201, 

X  if  X  is  2  (TWO),  THEN  MACHINE  CODE  FOR  THE  MATRIX  SOLU¬ 

TION  WILL  PE  GENERA-EO,  OTHERWISE,  NO  MACHINE  CODE  IS 
GENERATED,  THE  DEFAULT  VALUE  IS  2, 

X  1*^  X  IS  1  (ONE),  THE  ITERATION  TIMESTEP  CONTROL  IS  USED, 

IR  X  IS  2  (TWO),  THE  TRUNCATION-ERROR  "'IMESTEo  IS  USED, 

THE  DEFAULT  VALUE  IS  2,  IF  METHOO=GEAR  ANO  MAX0RD>2  THEN 
LVLTIM  IS  SET  TO  2  BY  SPICE. 

NAME  SETS  TME  NUMERICAL  INTEGRATION  METHOD  USED  BY  SPICE, 

POSSIBLE  NAMES  ARE  GEAR  OR  TRAPEZOIDAL,  THE  DEFAULT  IS 

trapezoidal, 

SETS  the  maximum  order  for  the  integration  METHOD  IF 
GEAP*S  variable-order  method  is  USED,  X  MUST  BE  BETWEEN 
2  AND  6,  THE  default  VALUE  IS  2, 


M  A  XO  R  D=  X 


,0Q  CAPO 


GENEPAL  P  OFW  *  OF 

THE  INCLUSION  OF  THIS  CARD  IN  AN  INPUT  DECK  WILL  F QPCE  SPICE  TO  DETERMINE 
the  DC  OPERATING  PO  I  N"^  DF  THE  CIRCUIT  WITH  INDUCTORS  SHORTED  AND  CAPACITORS 
OPENED#  NOTE;  A  DC  ANALYSIS  IS  AUTOMATICALLY  PERFORMED  PRIOR  TO  A  TRANSIENT 
ANALYSIS  TC  DETERMINE  THE  TRANSIENT  INITIAL  CONOITIONSt  AND  PRIOR  TQ  AN  AC 
small-signal  analysis  to  determine  THE  LINEARIZED,  SMALL-SIGNAL  MODELS  FOR 
NONLINEAR  DEV  ICES# 

S^ICE  PERFORMS  A  DC  OPERATING  pOIN^  ANALYSTS  IF  NO  OTHER  ANALYSES  ARE 
REQUESTED# 


#DC  CARD 

general  P'ORM  .DC  SRCNAM  VSTART  VSTOP  VINCR 

example  iOC  VTN  0#25  5#  0  C#  25 

THIS  CARD  defines  thE  DC  TRANSFER  CURVE  SOURCE  AND  SWEEP  LIMITS#  SRCNAM 
IS  THE  NAME  OF  AN  INDEPENDENT  VOLTAGE  OR  CURRENT  SOURCE#  VSTART,  VSTOP,  AND 
VINCR  ARE  THE  STARTING,  FINAL,  AND  INCREMENTING  VALUES  RE SPEC TI VEL Y •  THE  ABOVE 
example  will  CAUSE  THE  VALUE  OF  THE  VOLTAGE  SOURCE  VIN  TO  BE  SWEPT  FROM  0#25 
VOLTS  TO  5#0  VOLTS  IN  INCREMENTS  OF  0#25  VOLTS# 


.Tcr  CARD 

GENERAL  FORM  « TF  OUTVAF  INSRC 

example  .TF  V(E,3)  VIN 

THIS  CARD  defines  THE  SMALL-SIGNAL  OUTPUT  AND  INPUT  FOR  THE  DC  SMALL- 
SIGNAL  ANALYSIS#  OUTVAR  IS  THE  SMALL-SIGNAL  OUTPUT  VARIABLE  AND  INSPC  IS  THE 
SMALL-SIGNAL  INPUT  SOURCE#  IF  thIS  CARD  IS  INCLUDED,  SPICE  WILL  COMPUTE  THE 
DC  SMALL-SIGNAL  VALUE  OF  JHE  TRANSFER  FUNCTION  ( OUT PUT / I NPUT ) ,  INPUT  RESIS¬ 
TANCE,  AND  CU"T3UT  resistance.  FOR  THE  ABOVE  EXAMPLE,  SOICE  WOULD  COMPUTE  THE 
RATIO  OF  V(^,5)  TO  VIN,  THE  SMALL-SIGNAL  INPUT  RESISTANCE  AT  VIN,  AND  THE 
SMALL-SIGNAL  OU'^pU'^  RESISTANCE  MEASURED  ACROSS  NODES  5  AND  3# 


.SENS  CARD 

GENERAL  FORM  .SENS  DVl  I OV2  •••  ] 

example  .sens  VC9)  V(4,3)  V(17) 

1=  A  .SENS  CARD  IS  INCLUDED  IN  THE  INPUT  DECK,  SPICE  WILL  DETERMINE  THE 
DC  small-signal  sensitivities  of  EACH  SPECIFIED  OUTPUT  VARIABLE  WITH  RESPECT  TO 
EVERY  CIOCUIT  parameter®  NOTES  FOP  LARGE  CIRCUITS,  LARGE  AMOUNTS  OF  OUTPUT 
CAN  BE  GENERATED# 


**«*  .NCISE  CACD 


general  •nctse  outv  tnspc  nums 

example  .noise  V(5)  VIN  10 

THIS  CAPO  CONTPCLS  THE  NOISE  ANALYSIS  OF  THE  CIPCUTT.  OUTV  IS  A  VOLTAGE 
OUTPUT  VASIAELE  which  defines  the  summing  POINT.  INSRC  IS  the  name  OF  THE 
INDEPENDENT  VOLTAGE  OR  CURRENT  SOURCE  WHICH  IS  THE  NOISE  INPUT  RE*=^ERENCE.  NLWS 
IS  the  summary  INTERVAL.  SPICE  WILL  COMPUTE  THE  EQUIVALENT  OUTPUT  NOISE  AT 
the  SRSCl^IED  OUTPUT  AS  WELL  AS  THE  EQUIVALENT  INPUT  NOISE  AT  THE  SPECIFIED 
INPUT,  IN  ADDITION,  THE  C 0  NT R 1 5UT I CNS  OF  EVERY  noise  GENERATOR  IN  THE  CIRCUIT 
will  be  PRINTED  AT  EVERY  NUMS  FREQUENCY  POINTS  (THE  SUMMARY  INTERVAL),  IF  NUMS 
IS  ZERO,  NO  SUMMARY  PRINTOUT  WILL  BE  MACE. 

the  output  noise  AND  THE  EQUIVALENT  INPUT  NCISE  MAY  ALSO  BE  PRINTED  AND/OR 
PLOTTED  {SEE  THE  DESCRIPTION  OF  thE  .PRINT  AND  .PLOT  CARDS). 


****  ,TPiN  CARD 


GENERAL  POPM 

•  TP  AN 

TSTEP  TSTOP 

[TSTART  [TMAX]]  [UIC] 

E  XAMPLES 

•  TRAN 

INS  lOONS 

•  TRAN 

INS  lOOONS  5 

0  0  NS 

•  TRAN 

IONS  lUS  UIC 

TSTEP  IS  the  PRINTING  INCREMENT,  TSTOP  IS  THE  FINAL  TIME*  AND  TSTART  IS 
THE  INITIAL  ”IME,  Ip  TSTART  IS  OMITTED,  IT  IS  ASSUMED  TO  BE  ZERO.  THE 
TRANSIENT  ANALYSIS  ALWAYS  BEGINS  AT  TIME  ZERO,  IN  THE  INTERVAL  [ZERO,  TSTART), 
THE  CIRCUIT  IS  analyzed  (TO  REACH  A  STEADY  STATE),  BUT  NO  OUTPUTS  ARE  STORED. 

IN  THE  INTERVAL  [TSTART,  TSTOP],  THE  CIRCUIT  IS  ANALYZED  AND  OUTPUTS  ARE 
STORED.  tmax  IS  the  maximum  STEPSIZE  THAT  SPICE  WILL  USE  (DEFAULT  VALUE  IS 
T  STOP/ 50. 0  )  . 

UIC  (USE  INITIAL  CONDITIONS)  IS  AN  OPTIONAL  KEYWORD  WHICH  INDICATES  THAT 
THE  USER  DOES  NOT  WANT  SPICE  TC  SOLVE  FOR  THE  QUIESCENT  OPERATING  POINT  BEFORE 
beginning  the  transient  ANALYSIS.  IF  THIS  KEYWORD  IS  SPECIFIED,  SPICE  USES  THE 
values  specified  using  IC=...  on  THE  VARIOUS  ELEMENTS  AS  THE  INITIAL  TRANSIENT 
CONDITION  AND  PROCEEDS  WITH  THE  ANALYSIS. 


;  StstW 


■OUR  CARD 


GENERAL  PORM 


FOUR  PREQ  OVl  [ DV2  0V2 


example 


.FOUR  lOOKHZ  V(5) 


****  , AC  CAOD 


general  form 

•  AC 

DEC 

ND 

FSTAPT  FSTOP 

•  AC 

DCT 

NO 

FSTART  FSTOP 

•  AC 

LIN 

NP 

PSTAPT  FSTOP 

examples 

•  AC 

DEC 

10 

1  10KH2 

•  AC 

DEC 

20 

1  lOOKHZ 

•  AC 

LIN 

100 

1  100H2 

DEC  stands  pop  decade  VARIATIONr  AND  NC  IS  THE  NUMBER  OF  POINTS  PER 

DECADE.  OCT  stands  FOF  OCTAVE  VARIATION,  AND  NO  IS  THE  NUMBER  OF  POINTS  PER 

OCTAVE.  LIN  STANDS  FOR  LINEAR  VARIATION,  AND  NP  IS  THE  NUMBER  OF  POINTS. 

PSTART  is  the  starting  frequency,  AND  PSTOP  IS  THE  FINAL  FREQUENCY,.  IF  THIS 
CARD  IS  INCLUDED  IN  thE  DECK,  SRICE  WILL  PERFORM  AN  AC  ANALYSIS  OF  THE  CIRCUIT 
OVER  the  SPECIFIED  FREQUENCY  RANGE.  NOTE  THAT  IN  ORDER  FOR  THIS  ANALYSIS  TO  BE 
MEANINGFUL,  AT  least  ONE  INDEPENDENT  SOURCE  MUST  HAVE  SEEN  SPECIFIED  WITH  AN  AC 
VALUE. 


****  .DISTO  CARO 

GENERAL  FORM  .DISTO  RLCAD  I  INTER  fSKW2  [REFPWR  tSPWg]]]] 

EXAMdlE  .DISTO  RL  2  0.95  l.OE-3  0.75 

THIS  CARO  CON-^POLS  WHETHER  SPICE  WILL  COMPUTE  ThE  DISTORTION  CHARACTERIS¬ 
TICS  0=  THE  CIRCUIT  IN  A  SMALL-SIGNAL  MODE  AS  A  PART  OF  THE  AC  SMALL-SIGNAL 
SINUSOIDAL  steady-state  ANALYSIS.  thE  ANALYSIS  IS  PERFORMED  ASSUMING  THAT 
ONE  OR  TWO  signal  fpequENCIES  ARE  IMPOSED  AT  THE  INPUT;  LET  THE  TWO  FRE¬ 
QUENCIES  BE  R1  (THE  NOMINAL  ANALYSIS  FREQUENCY)  AND  F2  (=SKW2Wp1).  THE  PROGRAM 
then  computes  the  FOLLOWING  DISTORTION  MEASURES: 


HD  2 

-  THE 

MAGNI TUDE 

OF 

THE 

FREQUENCY 

component 

2FF1 

ASSUMING 

THAT 

F2 

IS  NOT  PRESENT* 

HD3 

-  the 

^^AGNITUDE 

OF 

the 

FREQUENCY 

component 

3FF1 

ASSUMI NG 

THAT 

F  2 

I S  NOT  PFESENT* 

5  IM2 

-  THE 

M  AGN'  ITUDE 

OF 

THE 

FREQUENCY 

COMPONENT 

FI  + 

F2* 

D  I  M2 

-  TWE 

MAGNITUDE 

OF 

THE 

FREQUENCY 

COMPONENT 

FI  - 

F2* 

D  TM? 

-  THE 

M  AGN  ITUDE 

OF 

THE 

FREQUENCY 

COMPONENT 

2FF  1 

-  F2* 

RLOAD  IS  THE  NAME  OF  THE  OUTPUT  LOAD  RESISTOR  INTO  WHICH  ALL  DISTORTION 
POWER  PRODUCTS  ARE  tq  EE  COMPUTED.  INTER  IS  THE  INTERVAL  AT  WHICH  THE  SUMMARY 
PRINTOUT  OF  TmE  contributions  OF  ALL  NONLINEAR  DEVICES  TO  THE  TOTAL  DISTORTION 
IS  TO  BE  PRINTED.  IF  OMITTED  OR  SET  TO  ZERO,  NO  SUMMARY  PRINTOUT  WILL  BE  MADE. 
REFDWR  IS  THE  OFerEp£NCE  POWER  LEVEL  USED  IN  COMPUTING  THE  DISTORTION  PRODUCTS. 
IF  omitted,  a  value  OF  I  MW  (THAT  IS,  DBM)  IS  USED.  SKW2  IS  TME  RATIO  OF  F2  TO 
FI.  IF  OMITTED,  A  VALUE  OF  0.9  IS  USED  (I.E.,  F2  =  0.9*fi),  sows  js  thE  AMP¬ 
LITUDE  OF  F2,  IF  omitted,  a  value  OF  I .0  IS  ASSUMED. 

THE  DISTORTION  MEASURES  H02,  HD3,  SIM2,  DIMS,  AND  DIM3  MAY  ALSO  BE  BE 
PRINTED  AND/OR  PLO'^TED  (SEE  THE  DESCRIPTION  OF  THE  .PRINT  AND  .PLOT  CARDS). 


BRINT  CAR 


CHNSRAL  PORW 


,oRINT  PPtybE  CVl  (0V2  •••  Ovel 


E  XAMPL  SS 


sPP  INT 
.op INT 
•PRINT 
•PRINT 

•PRINT 


TRAN  V(A)  nviN) 

AC  VM(4,2)  VR(7)  VR(e,2) 
DC  V(2)  I(VSRC)  V(23,17) 
NOISE  INOISE 
DISTO  HD3  SIM2(DB) 


THIS  CAPO  DEFINES  THE  CONTENTS  OF  A  TABULAR  LISTING  OF  ONE  TO  EIGHT  OUTPUT 
variables.  pRTYPE  IS  the  -^YFE  0*^  the  ANALYSIS  toe,  AC,  TRAN,  NOISE,  OR  OISTCR- 

TioNi  'OR  WHICH  THE  Specified  outputs  are  desired^  the  «orm  for  voltage  or 

CURRENT  OUTOO-^  VARIABLES  IS  AS  FOLLOWS: 


V( N1 1  ,N2  1  ) 


I  (  VXXXXXXX ) 


specifies  the  voltage  difference  between  nodes  N1 
AND  N2^  if  N2  (ANO  the  preceding  COMHA)  IS  OMITTED, 
GROUND  (0)  IS  assumed^  FOR  THE  AC  ANALYSIS,  FIVE 
ADDITIONAL  OUTPUTS  CAN  BE  ACCESSED  BY  REPLACING  THE 
letter  V  by: 


VR  “  REAL  PART 
VI  -  IMAGINARY  BART 

vM  -  magnitude 
VP  -  phase 

VCB  -  20FLOG10 (MAGNITUDE) 


SPECIFIES  THE  CURRENT  FLOWING  IN  THE  INDEPENDENT 

voltage  source  named  VXXXXXXX,  positive  current 
FLOWS  from  the  positive  node,  through  the  source,  to 

THE  NEGATIVE  NODE,  FOP  THE  AC  ANALYSIS,  THE  CORRES¬ 
PONDING  REPLACEMENTS  FOR  THE  LETTER  I  MAY  BE  MADE  IN 
THE  SAME  WAY  AS  DESCRIBED  FOR  VOLTAGE  OUTPUYS, 


OUTPUT  variables  ^OR  ywc  nq 1 S E  ANO  DISTORTION  ANALYSES  HAVE  A  0IFP£RENT 
PORM  FOOM  THAT  OF  THE  OTHER  ANALYSES,  "^HE  GENERAL  FORM  IS 


ov  r  tx ) ) 

aherE  OV  is  any  of  ONOISc  (OUTPUT  NOISE),  INOISE  (EQUIVALENT  INPUT  NOISE), 
HOE,  H03,  SIM2,  0IM2 ,  OR  DIMS  (SEE  DESCRIPTION  OF  DISTORTION  ANALYSIS),  ANO  X 
MAY  3E  any  of  : 


R  -  peal  part 
I  -  imaginary  BART 

M  -  magnitude  (DEFAULT  IF  NOTHING  SPECIFIED) 
P  -  phase 

OS  -  20*LCG1 0(maGNI TUDE) 


THUS,  SIMP  (OR  SIM2(M))  DESCRIBES  t»-e  MAGNITUDE  OF  THE  SIM2  DISTORTION  MEASURE, 
while  hD2(R)  DESCRIBES  "^HE  REAL  PART  OF  the  HD2  DISTORTION  MEASURE, 


,dlOT  C^=5D 


GENERAL  P'OPW 


»PLD^  PLTYPE  OVl  [ (  PLCl  t PH  1 1  )  ]  [ 0V2  [  ( PL02 , PH  1 2 )  ]  •••  0V8] 


EXAMPLES 


.PLOT  DC  V(4)  V(5)  V(l) 

.PLOT  TRAN  V(17t5)  (2,5)  I(VIN)  V(17)  (1,9) 

•  PLO'^  AC  VM(5)  VM(31t24)  VDB(5)  VP(5) 

•PLOT  DISTO  HC2  HC3(P)  51^2 


THIS  CARD  defines  THE  CONTENTS  0^  ONE  PLOT  OF  FROM  ONE  TO  EIGHT  OUTPUT 
VARIABLES.  PLTY=>E  IS  'ThE  TYPE  OF  ANALYSIS  (DC,  AC,  RA  N ,  NOISE,  OR  DISTORTION) 

POP  WHICH  THE  specified  OUTPUTS  APE  DESIRED.  THE  SYNTAX  FOR  THE  OVI  IS 
IDENTICAL  TO  THAT  FOP  THE  .PRINT  CARD,  DESCRIBED  ABOVE. 

THE  OPTIONAL  PLOT  LTMr’S  (PLO,PHl)  MAY  BE  SPECIFIED  AFTER  ANY  OF  THE 
OUTOUT  VARIABLES.  ALL  OUTPUT  VARIABLES  TO  THE  LEFT  OF  A  PAIR  OF  PLOT  LIMITS 
(PL0,PHI)  WILL  BE  PLOTTED  USING  THE  SAME  LOWER  AND  UPPER  PLOT  BOUNDS.  IF  PLOT 
LIMITS  ARE  NOT  SPECIFIED,  SPICE  WILL  AUTOMATICALLY  DETERMINE  THE  MINIMUM  AND 
MAXIMUM  VALUES  0^  ALL  OUTPUT  -VARIABLES  BEING  PLOTTED  AND  SCALE  THE  PLOT  TO  FIT. 
MORE  THAN  ONE  SCALE  WILL  BE  USED  IF  THE  OUTPUT  VARIABLE  VALUES  WARRANT  (I.E., 
MIXING  OUTPUT  variables  W  It h  VALUES  WHICH  ARE  ORDE RS-OF-M AGN TTUDE  DIFFERENT 
STILL  GIVES  READABLE  PLOTS). 


X  • 


the  overlap  of  TWO  OP  MORE  "GRACES  ON  ANY  ©LOT  IS  INDICATED  BY  THE  LETTER 


appendix  a 


example  data  decks 


the  following  deck  determines  the  DC  operating  point  and  SMAL.L-SI  gnal 

TRANSPEP  PUNCTTON  OF  A  SIMPLE  D I PFE ©E NT  I AL  •  PA  I R# 

simple  dipferenttal  pair 

VCC  7  0  12 
VEE  9  0  -12 
V  I  N  1  0 

RSI  1  2  IK 

RS2  6  0  1  K 

Q1  ?  2  4  MODI 

02  5  6  4  MODI 

PCI  T  3  1  OK 
PC2  7  E  1  OK 
PE  4  8  lOK 

•  Mooc^^  modi  NPN(?p  =  50  VA  =  50  IS  =  l«05-12  PR  =  l  00  > 

•  TP  V ( 5  )  V  IN 

•  END 


'^'-'E  FOLLOWING  DECK  DETERMINES  THE  DC  TRANSFER  CURVE  AND  THE  TRANSIENT 
PULSE  RESPONSE  OF  A  SIMPLE  PTL  INVERTER.  THE  INPUT  IS  A  PULSE  FROM  0  TO  5 
VOLTS  WI^H  DELAYt  PISE,  AND  FALL  TIMES  OF  2NS  AND  A  PULSE  WIDTH  OF  SONS.  THE 
TRANSIENT  INTERVAL  IS  0  TO  lOONS,  WITH  PRINTING  TQ  BE  DONE  EV  ER  Y  NANOSECOND. 

simple  PTL  INVERTER 
VCC  405 

VIN  !  0  PULS£(0  5  2NS  2NS  2NS  SONS) 

RE  1  2  lOK 
G  1  3  2  0  0 1 

RC  3  4  IK 
.  FUCT  DC  V  (3  ) 

•  PLOT  TRAN  V(  3)  (0,5) 

. PRINT  TRAN  V ( 3 ) 

.MODEL  01  NPN(8F=20  RB=100  TR=0.1NS  CJC=2PF) 

.DC  VIN  0  5  0.1 

.TRAN  INS  lOONS 

•  END 

Tme  FOLLOWING  DECK  DETERMINES  THE  AC  SMALL-SIGNAL  RESPONSE  A  ONE- 

TRANSISTOR  amplifier  OVEP  THE  FREQUENCY  RANGE  IHZ  TO  100MEGH2. 

0  NE-TR  ANS  I  S""0  R  AMPL  I  F  I  fp 
VCC  5012 

V  EE  6  0  -1 2 

V  IN  1  0  AC  I 

PS  1  2  IK 

01  3  2  4  X33 

RC  5  3  500 
PE  4  6  IK 
C  EYPASS  4  0  1  U^D 

•  PLOT  AC  VM(  *3  )  VP(  3) 

.AC  DEC  10  1H2  lOOMEGHZ 

.MODEL  X33  NPN{3p  =  30  RE^rSO  VA  =  20) 

.  END 


THE  FOLLOWING  DECK  SI^^ULATES  A  FOUF-5IT  EINAPY  ADDER,  USING  SEVERAL  SUB¬ 
CIRCUITS  TO  DESCRIBE  VARIOUS  PIECES  OF  THE  OVERALL  CIRCUIT. 

ADDER  -  A  BIT  A LL- N A ND-G A TE  BINARY  ADDE® 

***  SUBCIRCUIT  definitions 

o  SUBCKT  NANO  12  3  4 

*  nodes:  input(2),  output,  VCC 

Q1  9  5  I  QMOD 
DICLAMO  0  I  DMOD 
02  9  5  2  QMCD 
D2CLAMP  0  2  DMOD 
RB  4  5  4K 
R1  4  6  1.6K 

03  6  9  8  OMPD 

P2  8  0  IK 
RC  4  7  130 

Q4  6  10  OMOD 

DVBEDROP  10  3  DMOD 
QE  3  6  0  OMCD 

•  ENOS  NANO 

•SUBCKT  ONEPIT  123456 

nodes:  INPUT(2),  CARP.Y-IN,  output,  CARRY-OUT,  VCC 

XI  1  2  7  6  NAND 

X2  1  7  8  6  NAND 
X3  2  T  9  6  NAND 
X4  8  9  10  6  NAND 
X5  3  10  11  6  NAND 

X6  3  11  12  A  NAND 

X7  10  11  13  6  NAND 

X8  12  134  6  NAND 
X9  11  7  5  6  NAND 

•  ENDS  ONES  IT 

•SUBCKT  TWOBIT  123456700 

*  NODES:  INPUT  -  EIT0(2)  /  BIT1(2),  OUTPUT  -  BITO  /  BITl, 

*  CAPRY-IN,  CAPRY-OUT,  VCC 

X  I  1  2  7  5  1  0  9  ONES  I  T 
X2  3  4  10  6  8  9  ONEBIT 
•ENDS  TWOBIT 

•  SUBCK-  CUCURBIT  123^56-^89  10  11 

^  nodes:  INPUT  -  PIT0<2)  /  BIT1(2)  / 

*  OUTPUT  -  BITO  /  BITl  /  BIT2 

XI  1  2  3  4  9  10  1 3  1 6  15  TWOBIT 

X2  5  6  7  8  11  12  16  14  15  TWOBIT 

•  ENDS  PCUPBIT 

DEFINE  nominal  CIRCUIT 
•MODEL  DMOD  0 

.MODEL  QMOD  NdN(BF=75  PB=ipO  CJE=1PP  CJC=3PF) 
vyrc  99  0  DC  5V 

VINIA  1  0  PUL5E(  0  3  0  IONS  IONS  IONS  EONS) 


12  13  14  15 

BIT2 (2 )  /  8IT3( 2) , 

/  3IT3,  CARRY-IN,  CARRY-OUT,  VCC 


V  INI  B 

2 

0 

PUL  3  E  (  0 

3 

0 

1  ONS 

IONS 

2  0N  S 

1  OONS) 

VIN2A 

3 

0 

PULSE ( 0 

3 

0 

1  ONS 

1  ONS 

40NS 

200NS ) 

V  IN23 

A 

0 

P'JL  SE  (  0 

3 

0 

1  ON  S 

IONS 

SONS 

400NS) 

V  I  N3  A 

5 

0 

PULSE  (  0 

3 

0 

IONS 

IONS 

1  SONS 

BOONS) 

V  IN3B 

6 

0 

PUL  SE  (  0 

3 

0 

1  ONS 

IONS 

320NS 

1  600NS) 

V  TN4  A 

7 

0 

PULSE  (  0 

3 

0 

1  ONS 

IONS 

6  4  0N  S 

22  OONS) 

V  INAB 

3 

0 

PULSE  (  0 

3 

0 

IONS 

IONS 

1 2  SONS 

6400KS  > 

X  1  1  2 

4  5  3  7 

e  9 

10  11 

12  0 

1  3  99  I 

FQURBIT 

RBITO 

9 

0 

1  K 

PBITI  10  0  IK 
=551^2  11  0  IK 
T3  12  0  IK 
^CDU^  13  0  IK 

•  CLOT  -^PAN  V(l)  V(2>  V(3)  V(A)  V(5)  V(6)  VC?)  V(B) 

•  PLOT  TPANJ  V(  9)  VCIO)  veil)  V(12)  V(13) 

.^PINT  TPAN  V(l)  V(2)  V(3)  V(4)  V(5)  V(6)  V(7)  V(8) 

.point  tpam  v(o)  v(io)  veil)  V(12)  V(13) 

•  TPAN  INS  eAOONS 

±:k±  (PQO  thhco  WITH  MONEY  (AND  MEMORY)  TO  BURN) 

•  COT  ACCT  list  node  LIMPTS=6401 

•  END 


THE  ^’OLLOWING  DECK  SIMULATES  A  TR  AN  S  M I  S  S  1  ON- L  I  N5  INVERTER*  TWC  TRANSMIS¬ 
SION-LINE  ELEMENTS  ARE  REQUIRED  SINCE  TWO  PROPAGATION  MODES  ARE  EXCITED#  IN 
^HE  CASE  OF  A  COAXIAL  LINE,  THE  FIRST  LINE  ( T1 )  MODELS  THE  INNER  CONDUCTOR  WITH 
RESPECT  TO  the  SHIELD,  AND  THE  SECOND  LINE  ( T2 )  MODELS  THE  SHIELD  WITH  RESPECT 
TO  the  OUTSIDE  WORLD# 


TRANSMISSICN-LINE  INVER  ■"ER 
VI  1  0  P'J  L  SE  (  0  1  0  0  •  1  N  ) 

PI  1  2  50 

XI  20  OA  TLINE 

R2  4  0  50 

. SUSCKT  TLINE  1234 
T1  1  2  3  4  20  =  5  0  "^D=1«5NS 
T2  2  0  A  0  20=100  TD=1NS 
•ENDS  TLINE 
,TPAN  0#1NS  2  0NS 
#PLCT  ^RAN  V(2)  VC4) 

•  END 


A  P»ENO  !  X  B 


MriNLINeAR  DEPENDENT  SOURCES 


SPTCE  ALLOWS  CIRCUITS  TO  CONTAIN  DEPENDENT  SOURCES  CHARACTERIZED  BY  ANY  OF 
THE  FOUR  EQUATIONS 

I=F(V)  V=F(V)  I=F(I)  V=R(I) 

where  the  functions  MUST  BE  POLYNOMIALS!  AND  THE  ARGUMENTS  MAY  BE  MULTI¬ 
DIMENSIONAL.  THE  POLYNOMIAL  FUNCTIONS  ARE  SPECI^’IED  BY  A  SET  OF  COEFFICIENTS 
PO,  PI,  PN,  BOTH  the  number  OF  DIMENSIONS  AND  THE  NUMBER  OF  COEFFICIENTS 

ARE  arbitrary,  the  MEANING  OF  THE  COEFFICIENTS  DEPENDS  UPON  THE  DIMENSION  OF 
THE  POLYNOMIAL,  AS  SHOWN  IN  THE  FOLLOWING  EXAMPLES: 

SUPPOSE  THAH  the  FUNCTION  IS  ONE -O I  MENS I ON AL  (THAT  IS,  A  FUNCTION  OF  ONE 
ARGUMEN”),  THEN  THE  RUNCT  ION  VALUE  FV  IS  DETERMINED  BY  THE  FOLLOWING  EXPRES¬ 
SION  IN  FA  (THE  ^UNCTION  ARGUMENT); 

cv  =  DO  +  (di*fA)  +  (P2*FA*2)  +  (P3*FA*3)  +  (P4*RA»4)  +  (P3*FAt5)  +  ••• 

SUPPOSE  NOW  THAT  THE  FUNCTION  IS  TWO-DIMENSIONAL,  WITH  ARGUMENTS  FA  AND 

THEN  THE  FUNC^ICN  VALUE  FV  IS  DETERMINED  BY  THE  FOLLOWING  EXPRESSION: 

cv  =  PO  +  (PIFRA)  +  (PRFFB)  +  (P3*FA»2)  +  (P4WFAFFB)  +  (P5WFB»2> 

+  (P6*FAA3)  +  (P7*FA»2*FB)  +  (08*FA*FB*2)  ♦  (P9*FB*3)  ••• 

CONSIDER  NOW  THE  CASE  OF  A  THREE-D I  MENS  I CNAL  POLYNOMIAL  FUNCTION  WITH 
ARGUMENTS  FA,  FB,  AND  FC.  THEN  THE  RUNCTION  VALUE  FV  IS  DETERMINED  BY  THE 
FOLLOWING  EXPRESSION: 

RV  =  PO  (PI*FA)  (P2WFB)  (P3*^C)  +  (P4WRA*2)  +  (P5*FAWFB) 

+  (P6*FA*FC)  +  (P7*FB»2)  +  (P8*FB*FC)  ♦  (P9*FC»2)  +  (OlOFFA+3) 

+  (P11*FA*2*FB  )  -f  (R12FFAA2FFC  >  +  (P13*FA*FB»2  )  +  ( PI  4»F  AWFBFFC ) 

♦  (  PI  5*FA*FC»  2  )  (P16FFBA3)  +  (  P  I  7*  FB  ‘  2  WF  C  )  +  (  P 1  8*FB  ♦FC  ♦  2  ) 

(R19FFC43)  +  (P20*FA»4>  ••• 

NO'E;  IF  THE  RCLYNCMIAL  IS  ONE- D I  MENS  ICNAL  AND  EXACTLY  ONE  COEFFICIENT  IS 
SPECIFIED,  then  SPICE  ASSUMES  IT  TO  BE  Ri  (AND  PO  =  0.0),  IN  ORDER  TO  FACILI¬ 
TATE  ~HE  INPIJ"^  OF  LINEAR  controlled  SOURCES. 

FOR  ALL  FOUR  OF  THE  DEPENDENT  SOURCES  DESCRIBED  BELOW,  THE  INITIAL  CONDI¬ 
TION  parameter  is  described  as  optional,  if  not  SPECIFIED,  SPICE  ASSUMES  0,0, 
the  INITIAL  CONDITION  FOR  DEPENDENT  SOURCES  IS  AN  INITIAL  FGUESS*  FOR  THE  VALUE 
OP  THE  controlling  VARIABLE,  THE  PROGRAM  USES  THIS  INITIAL  CONDITION  TO  OBTAIN 
the  DC  OPERATING  POINT  OF  THE  CIRCUIT,  AFTER  CONVERGENCE  HAS  BEEN  OBTAINED, 

THE  program  CONTINUES  ITERATING  TO  OBTAIN  THE  EXACT  VALUE  FOR  THE  CONTROLLING 
VARIABLE.  HENCE,  TO  REDUCE  THE  COMPUTATIONAL  EFFORT  FOR  THE  DC  OPERATING 
POINT  (OR  IF  THE  POLYNOMIAL  SPECIFIES  A  STRONG  NONLINEARITY),  A  VALUE  FAIRLY 
CLOSE  TO  "^HE  ACTUAL  CONTROLLING  VARIABLE  SHOULD  BE  SPECIFIED  FOR  THE  INITIAL 
CONDI''  ION, 


VOL’ age-con~polled  current  sources 

GENERAL  GXXXXXXX  N+  N-  [OCLY(ND)]  NCl-t-  NCI-  •••  *>0  [PI  •••]  tlC=«*.] 

EXAMOi_=s  GX  1  0  5  2  0  0«1*^WHO 

GR  17  3  XT  3  0  IM  1 • EM  IC=2V 

GMLT  22  17  0CLY(2)  3  5  1  2  0  IM  17M  3»5U  1C=2.S,  1.3 

N+  AND  N-  ARE  Tt-E  POSITIVE  AND  NEGATIVE  NODES,  RESPECTIVELY.  CURRENT  FLOW 
IS  ^ROM  THE  POSITIVE  NODE,  THROUGH  THE  SOURCE,  TO  THE  NEGATIVE  NODE.  POLY(ND) 
ONLY  HAS  "^0  3E  SPECIFIED  IF  THE  SOURCE  IS  MULTI -0 1  MENSI  ONAL  ( ONE-OI  MENS  lONAL  IS 

the  DEFAULT).  IF  SOECIPIED,  ND  IS  THE  NUMEER  OF  DIMENSIONS,  WHICH  MUST  SE 

POSITIVE.  NCl-t-,  NC1-,  ...  ARE  THE  POSITIVE  AND  NEGATIVE  CONTROLLING  NODES, 
RESPECT  IVEL''.  ONE  PAIR  OF  NODES  MUST  BE  SPECIFIED  FOR  EACH  DIMENSION.  PO  ,  PI, 
=2,  ...,  ON  ARE  THE  POLYNOMIAL  COEFFICIENTS.  THE  (OPTIONAL)  INITIAL  CONDITION 
IS  •^HE  initial  guess  at  the  VALUE(S)  of  the  CONTROLLING  VOLTAGE(S).  IF  NOT 

S»£CIFT£d,  0.0  IS  ASSUMED.  THE  POLYNOMIAL  SPECIFIES  THE  SOURCE  CURRENT  AS  A 

OUNCTION  OF  THE  CONTROLLING  VOLTAGECS).  THE  SECOND  EXAMPLE  ABOVE  DESCRIBES  A 
CURRENT  SOURCE  WITH  VALUE 

I  =  1E-2*V(17,3)  +  1. 5E-3*V(1 7,3  )  *2 

note  THAT  SINCE  Hc  SOURCE  NODES  ARE  THE  SAME  AS  THE  CONTROLLING  NODES,  THIS 
SOURCE  ACTUALLY  MODELS  A  NONLINEAR  RESISTOR. 


*F**  VOLTAGE-CONTROLLED  VOLTAGE  SOURCES 

general  =0RM  EXXXXXXX  N-f  N-  rPCLY(NC)l  NCI  +  NCI-  ...  PO  [PI  •..]  IIC=..ol 

Examples  El  3  a  21  17  lo.s  2.1  1.75 

EX  17  0  P0LY(3)  13  0  15  0  17  0  0  1  1  1  I C= 1 .5 , 2 . 0 , 1 7 . 35 

N+  AND  N-  APE  the  POSITIVE  AND  NEGATIVE  NODES,  RESPECTIVELY.  POLY(ND) 

ONLY  HAS  TO  BE  SPECIFIED  IF  THE  SOURCE  IS  MULT  I -D I  MENS  I ONAL  ( ONE- D I  MEN S  ION AL  IS 

’HE  DEFAULT),  IF  SPECIFIED,  ND  IS  THE  NUMBER  OF  DIMENSIONS,  WHICH  MUST  BE 

POSITIVE.  NC1  +  ,  NC1-,  ...  ARE  The  POSITIVE  AND  NEGATIVE  CONTROLLING  NODES, 

RESPECTIVELY,  ONE  PAIR  OF  NODES  MUST  BE  S=ECIFIE0  ctqr  EACH  DIMENSION,  OO  ,  PI, 
=2,  ...,  PN  ARE  the  polynomial  coefficients,  THE  (OPTIONAL)  INITIAL  CONDI’ICN 
IS  ’HE  INITIAL  GUESS  AT  THE  VALUE(S)  OF  Th=  CONTROLLING  VCLTAGE(S),  IF  NOT 

specified,  0.0  IS  ASSUMED.  THE  POLYNDMIAL  SPECIFIES  THE  SOURCE  VOLTAGE  AS  A 

FUNCTION  OF  THE  CONTROLLING  VCLTAGE(S),  THE  SECOND  EXAMPLE  ABOVE  DESCRIBES  A 
voltage  SOURCE  WITH  VALUE 

V  =  V(13,0)  +  V(15,0)  +  V(17,0) 

(IN  other  words,  AN  IDEAL  VCLTAGE  SUMMER). 


****  current-controlled  current  sources 


GENERAL  ^DRM  RXXXXXXX  N+  N-  rPOLY(ND)l  VNJ  [ VN 2  •••]  PO  [PI  •••]  tIC=«««) 

EXAMO|_cs  FJ  12  10  VCC  IMA  1,3M 

FXFEP  12  20  VSENS  0  1 

N+  AND  N-  ARE  THE  POSITIVE  AND  NEGATIVE  NODES,  RESPECTIVELY.  CURRENT  FLOW 
!S  FPOM  THE  POSITIVE  NODE,  THROUGH  THE  SOURCE,  TO  THE  NEGATIVE  NODE*  POLY(ND) 
ONLY  HAS  TO  BE  SPECIFIED  IF  THE  SOURCE  IS  MULT  I -D I M ENS I ON AL  (ONE-DIMENSIONAL  IS 
THE  DEFAULT).  IF  SPECIFIED,  NO  IS  THE  NUMBER  OF  DIMENSIONS,  WHICH  MUST  BE 
POSITIVE.  VNl,  VN2,  ...  ARE  THE  NAMES  OF  VOLTAGE  SOURCES  THROUGH  WHICH  THE 
CCNTPOLLING  CURRENT  FLOWS!  ONE  NAME  MUST  BE  SPECIFIED  FOR  EACH  DIMENSION.  THE 
DIRECTION  OF  POSITIVE  CONTROLLING  CURRENT  FLOW  IS  FROM  THE  POSITIVE  NODE, 
through  the  SOURCE,  TO  THE  NEGATIVE  NODE  OF  EACH  VOLTAGE  SOURCE.  P  O,  PI, 

=2,  ...,  PN  ARE  the  POLYNOMIAL  COEFFICIENTS.  THE  (OPTIONAL)  INITIAL  CONDITION 

IS  the  initial  GUESS  AT  THE  VALUE! S)  OF  THE  CONTROLLING  CURRENT! S)  (IN  AMPS), 
ic  MOT  S3=cTpIED,  0.0  IS  ASSUMED.  THE  FCLYNCMIAL  SPECIFIES  THE  SOURCE  CURRENT 
AS  A  function  0'=  'HE  CONTROLLING  CUKRENT(S).  THE  ^IRST  EXAMPLE  ABOVE  DESCRIBES 
A  CURRENT  SOURCE  WITH  VALUE 

I  =  lE-2  +  1 .3E-2*I ( VCC ) 


****  CURRENT-CONTROLLED  VOLTAGE  SOURCES 

GENERAL  form  HXXXXXXX  N+  N-  [PCLY(ND)J  VNl  [ VN2  ...1  PO  [PI  **•]  [1C=...] 

EXAMPLES  HXY  13  20  P0LY(2)  VINl  VIN2  00001  IC=0.5  1.3 

HR417VX0  01 

N+  AND  N-  ARE  the  POSITIVE  AND  NEGATIVE  NODES,  RESPECTIVELY.  POLY(ND) 

ONLY  HAS  TO  BE  SPECIFIED  IF  THE  SOURCE  IS  MUL TI -D I  MENS I  ON AL  (ONE-DIMENSIONAL  IS 
-HE  DE'^AULT),  IP  SPECIFIED,  ND  IS  thE  NUMBER  OF  DIMENSIONS,  WHICH  MUST  BE 
POSITIVE.  VNl,  VN2,  ...  ARE  THE  NAMES  CF  VOLTAGE  SOURCES  THROUGH  WHICH  THE 
CONTROLLING  CURREN-  FLOWS!  CNE  NAME  MUST  BE  SPECIFIED  FOR  EACH  DIMENSION.  THE 
DIRECTION  OF  POSITIVE  CONTROLLING  CURRENT  FLOW  IS  FROM  THE  POSITIVE  NODE, 
through  tmc  source,  TO  THE  NEGATIVE  NODE  OF  EACH  VOLTAGE  SOURCE.  PO,  PI, 

°2,  PN  ARE  the  POLYNOMIAL  COEFFICIENTS.  THE  (OPTIONAL)  INITIAL  CONDITION 

IS  the  INITIAL  GUESS  AT  THE  VALUE ( S >  OF  THE  CONTROLLING  CURRENTfS)  (IN  AMPS). 

IF  not  specified,  0.0  IS  ASSUMED.  thE  POLYNOMIAL  SPECIFIES  THE  SOURCE  VOLTAGE 
AS  A  FUNC-ICN  OF  THE  CONTROLLING  CUFRENT(S).  THE  FIRST  EXAMPLE  ABOVE  DESCRIBES 
A  voltage  SOURCE  WITH  VALUE 


V 


I  ( VI Nl  ) *  I ( V  IN2 ) 


SPICZ2  REPORT 


LINKED  LIST  SPECIFICATIONS 


LIST  ELEMENTS 


13>2  LinXed  List  Specifications 

Each  list  element  generally  contains  both  integer  and  real 
data*  Even  though  both  data  types  require  only  one  word  of 
memory  on  the  CDC  6400  computer*  separate  subscripts  are  used  to 
access  the  two  types*  All  integer  data  is  referenced  using  the 
array  NOOPLC;  all  real  (and  character)  data  is  accessed  using 
the  (equivalenced)  array  VALUE*  The  VALUE-subscript  for  the 
first  real  value  is  stored  in  the  integer  part  of  the  list  ele¬ 
ment  and  is  called  LOCV  (LOCH  for  device  models);  the  NODPLC- 
subscript  is  called  LOC* 

In  the  detailed  list-element  structure  definitions  which 
follow,  notation  is  defined  only  on  first  use* 
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RESISTOR 


13.2.1  RESISTOR 


ID  »  1 


- 

1 

suDckt  info 

LOC  + 

0 

next-pointer 

1 

LOCV 

LOCV-H 

0: 

element  name 

-*■ 

2 

nl 

l: 

g(TEMP) 

+ 

3 

n2 

-f- 

2: 

r(TNOM) 

+ 

4 

(nl ,n2) 

3: 

temp,  coefficient 

1 

5 

(n2 , nl ) 

4 : 

temp,  coefficient 

2 

+ 

6 

(nl , nl ) 

7 

(n2,n2) 

Comments  s 

1)  "subcict  info"  is  used  to  indicate  subcircuit  relationships 
(see  Sections  5.3  and  6.3  for  details).  Briefly* 

a)  if  the  element  is  part  of  the  nominal  circuit 
description,  "subckt  info"  is  zero. 

b)  if  the  element  is  contained  within  a  subcircuit 
definition,  "subcxt  info"  is  the  element  ID  (1  for 
resistors ) . 

c)  if  the  element  is  added  to  the  circuit  as  a 
result  of  subcircuit  expansion,  "suockt  info"  is  a 
pointer  to  the  "X"  element  which  caused  the  expan¬ 
sion. 

2)  "next-pointer"  points  to  (is  the  NODPLC-subscript  of)  the 
next  element  of  the  same  ID;  if  there  is  no  next  element, 
"next-pointer"  is  zero.  However,  if  this  element  is  a  part  of  a 
subcircuit  definition,  then  "next-pointer"  points  to  the  next 
element  within  the  definition,  regardless  of  ID. 

3)  "LOCV"  points  to  the  real-valued  storage  for  the  element. 

4)  "ni"  is  element  node  number  i.  During  READIN,  this  is  the 
number  read  from  input;  after  ERRCHK,  this  entry  is  replaced  by 
an  index  into  the  JUNODE  array  (the  compact  renumoered  node 
list ) . 

5)  The  notation  “(a,b)"  means  a  pointer  to  matrix  location 
(a,b):  the  a'th  row,  b'th  column  entry. 

6)  "element  name"  is  the  element  name,  left- justified,  with 
blank  fill  to  8  characters. 

7)  "g(TEMP)"  is  the  element  conductance,  adjusted  for  the  value 
of  TEMP. 
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RESISTOR 


8)  “r(TNO['l)"  is  the  input  element  resistance  {assumed  to  be 

TNOM  degrees). 


at 
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CAPACITOR 


13.2.2  CAPACITOR 


ID  »  2 


-  1 

subckt  info 

LOC+  0 

next-pointer 

+  1 

LOCV 

LOCV  + 

0: 

element  name 

+  2 

nl 

l: 

computed  element  value 

+  3 

n2 

+ 

3: 

initial  condition 

+  4 

function  code 

+ 

3: 

argument  vector 

+  5 

(nl . n2 ) 

+  6 

(n2.nl) 

+  7 

tp(function  coefficients) 

+  8 

LXi  offset 

LXi  + 

OS 

q(capacitor ) 

+  9 

exponent  vector 

+ 

l: 

i  (capacitor) 

+  10 

(nl.nl ) 

+  11 

(n2 . n2  ) 

Comments : 

1)  "function  code"  is  zero  for  "polynomial",  the  only  function 
currently  implemented. 

2)  The  notation  "tp ( something ) "  means  "a  table  pointer  to  a 
table  which  contains  'something'". 

3)  "LXi  offset"  is  the  offset  for  this  element  into  any  of  the 
LXi  tables  (LXO.  LXI.  etc.)  which  are  used  during  analysis  to 
contain  intermediate  analysis  results. 

4)  "argument  vector"  and  "exponent  vector"  are  used  by  the  func¬ 
tion  evaluation  routines. 

5)  "q (element)"  means  the  charge  stored  in  element. 

6)  "i (element)"  means  the  current  flowing  in  element. 
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13.2.3  INDUCTOR 


ID  * 

3 

- 

1 

subcxt  info 

LOC+ 

0 

next-pointer 

•f 

1 

LOCV 

LOCV  + 

0: 

element  name 

+ 

2 

nl 

+ 

1  * 

computed  element  value 

+ 

3 

n2 

+ 

2.- 

initial  condition 

+ 

4 

function  code 

3: 

argument  vector 

■+ 

5 

IBR 

+ 

6 

(nl.IBR) 

+ 

7 

(n2,IBR) 

+ 

8 

(IBR,nl) 

9 

(IBR,n2) 

+10 

tp(function  coefficients) 

+11 

LKi  offset 

LXi  + 

0: 

phi ( inductor ) 

+  12 

exponent  vector 

+ 

i; 

V  ( inductor ) 

+  13 

(IBR, IBR) 

Comments : 

1)  "I3R"  is  the  equation  number  for  the  inductor  current. 

2)  "phi (element ) "  is  the  flux  in  element. 

3)  "v(element)"  is  the  voltage  across  element. 
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13.2.4  MUTUAL  INDUCTANCE 


ID  = 

4 

- 

1 

subckt  info 

LOC  + 

0 

next-pointer 

-»■ 

1 

LOCV 

LOCV+  0: 

element  name 

-»■ 

2 

ptr (LI ) 

+  l; 

value 

-»■ 

3 

ptr (L2 ) 

4 

(L1/L2) 

-»■ 

5 

(L2,L1) 

Comments: 

1)  "ptr(Li)"  means  a  pointer  to  one  of  the  inductor  elements 
which  this  element  is  coupling.  During  READIN»  this  word  is  an 
index  into  the  lUNSAT  table;  after  ERRCHK.  it  points  directly  to 
the  inductor. 

2)  "(LifL^)"  means  the  matrix  location  (A»B)  where  A  is  the 
equation  number  for  the  current  in  Li»  and  B  is  the  equation 
number  for  the  current  in  Lj. 

3)  During  READIN#  "value"  is  Kf  the  coefficient  of  coupling; 
after  ERRCHK»  "value"  is  M,  the  mutual  inductance. 
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source:  I 


G(V) 


13.2.5  VOLTAGE-CONTROLLED  CURRENT  SOURCE 
ID  *  5 


-  1:  subcict  info 
LOC+  0:  next-pointer 

+  1:  LOCV  LOCV+  0:  element  name 

+  2 :  nl 
+  3:  n2 

+  4:  dimension  of  function 
+  5;  function  code 
+  6:  tp (controlling  nodes) 

+  7:  tp (matrix  locations) 

+  8:  tp(function  coefficients) 

+  9:  tp (argument  vector) 

+10:  tp (exponent  vector) 

+11:  tp(initial  conditions) 

+12:  LXi  offset  LXi  +  0:  i (source) 

2  values/dimension:  +  l:  controlling  v 

+  2:  di (source)/dv (control ) 
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source:  V 


E{V) 


13.2.6  VOLTAGE-CONTROLLED  VOLTAGE  SOURCE 
ID  =  6 


-  1 
LOC+  0 
+  1 
+  2 
+  3 
+  4 
+  5 
+  6 
+  7 
+  8 
+  9 
+  10 
+11 
+12 
+13 


subckt  info 
next-pointer 

LOCV  LOCV+  0:  element  name 

n+ 

n- 

dimension  of  function 

function  code 

IBR 

tp (controlling  nodes) 
tp (matrix  locations) 
tp(function  coefficients) 
tp (argument  vector) 
tp (exponent  vector) 
tp( initial  conditions) 

LXi  offset  LXi  +  0:  v(source) 

+  1:  i( source) 

2  values/dimension:  +  2:  controlling  v 

+  3:  dv(source)/dv(control) 


Comments : 


1)  "IBR"  is  the  equation  number  for  the  current  through  this 
element. 
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F(I  ) 


13. 


.7  CURRENT-CONTROLLED  CURRENT  SOURCE 
ID  =  7 


-  1 

LOC+  0 
+  1 
+  2 
+  3 
+  4 
+  5 
+  6 
+  7 
+  8 
+  9 
+  10 
+  11 
+  12 


subckt  info 
next-pointer 

LOCV  LOCV+  0:  element  name 

nl 

n2 

dimension  of  function 
function  code 

tp(ptrs  to  controlling  currents) 
tp (matrix  locations) 
tp (function  coefficients) 
tp (argument  vector) 
tp (exponent  vector) 
tp (initial  conditions) 

LXi  offset  LXi  +  0:  i(source) 

2  values/dimension:  +  l:  controlling  i 

+  2:  di ( source ) /di (control ) 
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source: 


V  * 


H(I) 


13.2.8  CURRENT-CONTROLLED  VOLTAGE  SOURCE 
ID  =  8 

-  1:  subckt  info 
LOC+  0;  next-pointer 

+  1:  LOCV  LOCV+  0:  element  name 

+  2 :  n+ 

+  3 :  n- 

4:  dimension  of  function 
+  5:  function  code 
+  6:  IBR 

+  7:  tp(ptrs  to  controlling  currents) 

+  8:  tp (matrix  locations) 

+  9:  tp( function  coefficients) 

+10:  tp (argument  vector) 

+11:  tp (exponent  vector) 

+12:  initial  conditions 

+13:  LXi  offset  LXi  +  0:  v(source) 

+  1:  i(source) 

2  values/dimension:  +  2:  controlling  i 

+  3:  dv ( source )/di (control ) 
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13.2.9  INDEPENDENT  VOLTAGE  SOURCE 


ID  »  9 


- 

1 

subckt  info 

LOC  + 

0 

next-pointer 

1 

LOCV 

LOCV+ 

0: 

element  name 

2 

n+ 

l: 

dc/transient  value 

+ 

3 

n- 

+ 

2: 

ac  value:  magnitude 

+ 

4 

function  code 

+ 

3: 

ac  value:  phase 

-H 

5 

tp( function  coefficients) 

+ 

6 

IBR 

+ 

7 

(nl ,IBR) 

+ 

8 

{n2,IBR) 

+ 

9 

(IBR.nl ) 

+  10; 

(IBR,n2) 

Comments : 

1)  "function  code"  indicates  which  built-in  function  to  use  for 
this  source,  from  the  following  possibilities: 


coqe 


built-in 


0 

1 

2 

3 

4 

5 


<no  function  specified> 

pulse 

sine 

exponential 

piece-wise  linear  (PWL) 
single-frequency  fm  (SFFM) 


2)  "IBR"  is  the  equation  number  for  the  current  flowing  in  this 
element • 
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13-2.10  INDEPENDENT  CURRENT  SOURCE 


ID  =  10 


- 

1: 

suocKt  info 

LOC  + 

0: 

next-pointer 

+ 

l: 

LOCV 

LOCV  + 

0: 

element  name 

+ 

2: 

nl 

+ 

l: 

dc/transient  value 

+ 

3: 

n2 

+ 

2: 

ac  value:  magnitude 

+ 

4 : 

function  code 

+ 

3: 

ac  value:  phase 

+ 

5: 

tp (function  coefficients) 

Comments : 


1)  "function  code"  indicates  which  built-in  function  to  use  for 
this  source*  from  the  following  possibilities: 


function  code 


0 

1 

2 

3 

4 

5 


<no  function  specified> 

pulse 

sine 

exponentia 1 

piece-wise  linear  (PNL) 
single-frequency  fm  (SFFM) 
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13.2.11  DIODE 


ID 


11 


-  1 

subcxt  into 

+  0 

next-pointer 

+  1 

LOCV 

LOCV  + 

0 : 

element  name 

+  2 

np 

+ 

1 : 

area  factor 

+  3 

nn 

+ 

2: 

IC:  vd 

+  4 

np' 

+  .5 

mp 

+  6 

off 

+  7 

(np/ np' ) 

+  8 

(nn/ np ' ) 

+  9 

(np' / np ) 

+  10 

(np ' / nn ) 

+  11 

LXi  offset 

LXi  + 

0: 

v(diode ) 

+  12 

LDO  offset 

l: 

i ( diode ) 

+  13 

(np / np ) 

+ 

2: 

geq 

+  14 

(nn/nn) 

3: 

g(diode  capa 

+  15 

(np ' 1 np' ) 

4  : 

i (diode  capa 

Comments : 


1)  If  the  diode  has  no  extrinsic  resistance  (RS  >  0 

corresponding  device  model)/  then  np  =  np '  ■  the  p-doped 
the  diode/  and  nn  *  the  n-doped  side.  If  RS  is  nonzero/  then  it 
is  modeled  as  a  resistance  between  np  and  np'/  with  np'  being 
true  p-doped  side  of  the  diode- 


in  the 
side  of 


the 


2)  "mp"  is  a  pointer  to  the  device  model  for  this  element* 

3)  "off"  is  zero  unless  the  element  was  specified  as  "off"  in 
the  circuit  description  (in  which  case  its  value  is  1). 

4)  "IC:"  means  "initial  condition  specification". 

5) .  "LDO  offset"  is  the  offset  for  this  element  into  the  distor¬ 
tion  analysis  working  storage  taole  LDO. 
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13.2.12  BJT 
ID  =  12 


LOC+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 


+  10 
+11 
+12 
+13 
+  14 
+  15 
+16 
+  17 
+18 
+19 
+20 
+21 


subcict  info 

next-pointer 

LOCV 

no 

nb 

ne 

nc' 

nb' 

ne ' 

mp 

off 

( nc »  nc ' ) 

(nb/ nb' ) 

( ne I ne ' ) 

( nc ' / nc ) 

( nc ' / nb' ) 

( nc ' I ne ' ) 

( nb' / nb ) 

(nb' / nc' ) 
(nb' / ne' ) 
(ne' /ne) 

(ne' / nc' ) 
(ne' /nb' ) 


LOCV+  0: 
+  l: 
+  2: 


element  name 
area  factor 
I C :  vbe 
3:  IC:  vce 


+22 

LXi 

offset 

LXi  + 

0 

Vbe 

+  23 

LDO 

offset 

1 

vbc 

+24 

(nci 

nc ) 

+ 

2 

ic 

+  25 

(nb/ 

nb ) 

3 

ib 

+  26 

(nei 

ne ) 

+ 

4 

gpi 

+  27 

(nc' 

/ nc'  ) 

+ 

5 

gmu 

+  23 

( no' 

/nb' ) 

6 

gmo 

+29 

( ne ' 

'  / ne '  ) 

+ 

7 

go 

+ 

8 

g(cbe) 

+ 

9 

i ( coe ) 

+  1 

0 

q(CDC) 

+11 

i (coc) 

+1 

2 

C ( CCS ) 

+13. 

i (CCS ) 

Comments : 


1)  The  element  nodes  nc,  nb,  and  ne  are  the  collector#  base/  ana 
emitter  nodes/  respectively.  If  any  of  the  extrinsic  resistances 
(RB/  RC/  or  RE)  are  nonzero/  then  the  corresponding  resistor  is 
included  between  nodes  nx  and  nx'  (for  the  appropriate  'x'). 
Otherwise/  nx'  is  the  same  node  as  nx. 
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13.2.13  JFET 


ID 

»  13 

-  l: 

suhckt  info 

LOC+  0: 

next-pointer 

+  1 

LOCV 

LOCV+ 

0: 

element  name 

+  2 

nd 

1: 

area  factor 

+  3 

ng 

+ 

2: 

IC_:  vds 

+  4 

ns 

3: 

IC":  vgs 

+  5 

nd' 

+  6 

ns' 

+  7 

mp 

+  8 

off 

+  9 

(nd.nd' ) 

+  10 

(ng/ nd' ) 

+  11 

(ng / ns ' ) 

+  12 

(ns/ ns ' ) 

+  13 

(nd' / nd ) 

+  14 

(nd' / ng) 

+  15 

(nd ' / ns ' ) 

+  16 

(ns ' / ng ) 

+  17 

( ns ' f  ns ) 

+  18 

(ns' «nd' ) 

+19 

LXi  offset 

LXi 

+ 

0 

vgs 

+  20 

(nd / nd ) 

1 

vgd 

+21 

(ng/ ng) 

+ 

2 

ig  (gate) 

+  22 

(ns / ns ) 

3 

id  (drain) 

+  23 

(nd' /nd' ) 

4 

i (gate-to-drain ) 

+24 

(ns ' / ns ' ) 

-1- 

5 

gm 

+ 

6 

gds 

-K 

7 

ggs 

+ 

8 

ggd 

9 

q (cgs ) 

+  10 

i  (cgs ) 

+11 

q(cgd) 

+  12 

i  (cgd ) 

Comment 

5 : 

1)  The 

element  nodes  nd,  ng» 

and  ns  are 

the  drain/  gate/ 

source , 

respectively. 

and 
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13.2.14  MOSFET 


ID  =  14 

-  1 

subckt  info 

LOC+  0 

next-pointer 

+  1 

LOCV 

+  2 

nd 

+  3 

ng 

+  4 

ns 

+  5 

nb 

+  6 

nd' 

+  7 

ns ' 

+  8 

mp 

+  9 

off 

+  10 

(nd»  nd ' ) 

+  11 

(ng»  nb) 

+  12 

(ng»nd' ) 

+  13 

(ng, ns ' ) 

+  14 

(ns , ns ' ) 

+  15 

(nb, ng ) 

+  16 

(nb, nd ' ) 

+  17 

(nb, ns ' } 

+18 

(nd ' , nd ) 

+  19 

(nd ' , ng ) 

+  20 

(nd' , nb) 

+  21 

(nd ' , ns ' ) 

+  22 

(ns ' , ng ) 

+  23 

(ns ' , ns ) 

+24 

(ns ' , nb) 

+  25 

(ns ' , nd ' ) 

+  26 

LXi  offset 

+  27 

(nd, nd ) 

+  28 

(ng,ng) 

+  29 

(ns , ns ) 

+  30 

(nb, nb) 

+  31 

(nd' ,nd' ) 

+  32 

(ns ' , ns ' ) 

LOCV+  0 
+  1 
+  2 
+  3 
+  4 
+  5 
+  6 
+  7 
+  8 
+  9 
+  10 


element  name 
channel  length 
channel  width 
drain  diffusion  area 
source  diffusion  area 
IC:  vds 
IC:  vgs 
IC:  VOS 

devmod 

von 

vdsat 


i  +  0 

vod 

+  1 

vbs 

+  2 

vgs 

+  3 

vgd 

+  4 

id 

+  5 

ibs 

+  6 

ibd 

+  7 

gm 

+  8 

gds 

+  9 

gmos 

+  10 

gbd 

+  11 

gbs 

+  12 

q (cbd ) 

+  13 

X (cbd) 

+  14 

g (cbs ) 

+  15 

i (cbs ) 

+  16 

q (cgs ) 

+  17 

i  (cgs ) 

+  18 

q  (cgd) 

+19 

i  (cgd ) 

+  20 

q (cgb) 

.  +21 

i (cgb) 
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Conunents : 

1)  The  element  nodes  nd»  ng,  ns»  and  nb  are  respectively 
draini  gate«  source*  and  bulk  (substrate)* 

2)  "devmod"  is  the  device  mode:  +1  (-1)  for  normal  (inverse) 

3)  "von"  is  the  adjusted  threshold  voltage  for  the  device. 

4)  "vdsat"  is  the  saturation  voltage  for  the  device. 


the 
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13.2.15  TRANSMISSION  LINE 


ID  *  17 

-  1 

subckt  info 

LOC+  0 

next-pointer 

+  1 

LOCV 

LOCV  + 

0: 

+  2 

nl 

l: 

+  3 

n2 

+ 

2: 

+  4 

n3 

3: 

+  5 

n4 

4: 

+  6 

nil 

5: 

+  7 

ni2 

6: 

+  8 

IBRl 

7: 

+  9 

IBR2 

+ 

8: 

+  10 

(nl #  nl  ) 

+11 

(nl / nil ) 

+  12 

(n2,IBRl ) 

+  13 

(n3,n3) 

+  14 

(n4 ,IBR2 ) 

+15 

( nil  I nl ) 

+  16 

(nil <  nil ) 

+  17 

(nil , IBRl ) 

+18 

(ni2 , ni2 ) 

+  19 

(ni2,IBR2) 

+20 

(IBRl ,n2 ) 

+  21 

(IBRl,n3) 

+  22 

(IBRl ,n4 ) 

+23 

( IBRl , nil ) 

+  24 

(IBRl ,IBR2 ) 

+  25 

(IBR2,nl) 

+26 

(IBR2,n2 ) 

+27 

(IBR2 ,n4 ) 

+28 

(IBR2 ,ni2 ) 

+  29 

(IBR2 ,IBR1 ) 

+  30 

LTD  offset 

LTD  + 

0: 

+  31 

(n3 » ni2 ) 

+ 

1 : 

+  32 

(ni2 ,n3 ) 

Comments : 


1]  The  model  for  this  element  is  described 
SPICE  circuit  elements  (although  the  program 
tional  dependencies): 

R1  nl  nil  ZO 

VI  nil  n2  delay (v (n3 » n4 ) *TD ) +delay (i (v2 ) 
R2  n3  ni2  ZO 

V2  ni2  n4  de lay ( v (nl . n2 ) ,TD ) +delay (i ( vl ) 


TRANSMISSION  LINE 


element  name 
20  ^ 
td 

excitation:  IBRl 
excitation:  IBR2 
IC:  v(port  1) 

IC:  i(port  1) 

IC:  v(port  2) 

IC:  i(port  2) 


past  value  (LOCV  +  3) 
past  value  (LOCV  +  4) 


below  in  terms  of 
does  not  allow  func- 


,TD )  *Z0 
,TD ) *Z0 
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13.2.16  SUBCIRCUIT  CALL  ('X'  element) 

ID  «  19 

-  1:  subckt  info 
LOC+  0:  next-pointer 

+  1:  LOCV  LOCV+  0:  element  name 

+  2:  tp('X'  element  nodes) 

+  3:  ptr  to  suJDcircuit  description 

Comments : 

1)  Subcircuit  structures  and  manipulations  are  sketched  as  a 
part  of  the  description  of  suucircuit  definitions  in  Section 
13.2.17 . 
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13.2.17  SUBCIRCUIT  DEFINITION 

ID  =  20 

-  l: 

LOC+  0: 

+  l: 

+  2: 

+  3: 

Comments : 

1)  Consider  the  following  SPICE  input  (not  contained  within  any 
enclosing  subcircuit  definition): 

XI  1  2  3  SUBl 
.SUBCKT  SUBl  17  13  4 
XI  17  6  4  SUB2 
R1  17  13  IK 
Q1  1  4  17  A 
.MODEL  A  NPN() 

•ENDS  SUBl 

At  the  end  of  READIN  (oefore  ERRCHK)»  the  internal  data  struc¬ 
tures  representing  this  input  would  look  something  like  the  fol¬ 
lowing: 


subckt  info 
next-pointer 

LOCV  LOCV+  0:  subcircuit  name 

tp( subckt  definition  nodes) 
ptr  to  definition  element  list 


LOCATE (19 ) 


—  > 


0 


0 


(LOCV) 


(nodlst ) 


( subckt ) 


— >  "XI" 


—  > 


l: 

2: 

3: 


> 


"SUBl" 


(in  lUNSAT  table) 
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LOCATE (20)  — > 


( LOCV ) 
(nodlst ) 
ptr ( def ) 


— >  “SUB1‘ 


—  > 


17 
13 
.  .4 


(ID)  19 

(ID)  1 

nxt-pntr 

--> 

nxt-pntr 

<"X1"> 

<"R1"> 

( def n ) 

(def n ) 

—  > 


2)  After  subcircuit  expansion  (controlled  by  subroutine  SUBCKT 
in  overlay  ERRCHX)»  copies  of  the  elements  on  the  subcircuit  ele¬ 
ment  list  (LOC  +  3)  are  added  (with  appropriately  modified  node 
numbers)  to  the  element  lists.  For  these  added  elements/  "subckt 
info"  is  a  pointer  to  the  "X"  element  which  caused  the  given  ele¬ 
ment  to  be  added. 
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13.2.18  DIODE  MODEL 
ID  =  21 

-  1:  subcxt  info 
LOC+  0:  next-pointer 

+  l:  LOCV  LOCV+  0 

+  2:  <unused>  +  1 

+  2 
+  3 
+  4 
+  5 
+  6 
+  7 
+  8 
+  9 
+  10 
+  11 
+12 
+  13 
+  14 
+  15 
+16 
+  17 
+  18 

Comments : 

1)  The  model  parameter  keywords  are  defined 
(Section  13.1). 

2)  The  notation  "f(argsy"  refers  to  some  va 
tion  of  the  variable (s)  "args")  which 
subroutine  MODCHK. 


model  name 
IS 
RS 
N 

TT 

CJO 

PB 

M 

EG 

PT 

KF 

AF 

FC 

BV 

IBV 

f  (FC) 

f  (FC) 

f  (FC) 

vcrit 


in  the  User's  Guide 


lue  (which  is  a  func- 
is  evaluated  once  in 


3)  One-time  only  preprocessing  of  model  parameters  is  done  in 
subroutine  MODCHK  in  the  ERRCHK  overlay.  For  the  diode  model » 
the  preprocessing  is  as  follows: 


location  replaced  ^y. 


LOCV+  2 
+  12 
+  13 
+  15 
+  16 
+17 
+  18 


1/RS  (0  if  RS  =  0) 

FC*PB 

computed  start  of  reverse  exponential 
PB*(1-(1-FC)''(1-M)  )/(l-M) 

(1-FC)“ (1+M) 

1-FC* (1+M) 

N*VT*ln(N*VT/(sqrt(2)*IS) ) 
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4)  SPICE  does  not  generate  an  additional  node  for  the  diode 
resistance  if  RS  is  zero;  insteadr  the  np'  node  is  set  equal  to 
the  np  node  (see  Section  13.2.11).  This  processing  of  np»  to¬ 
gether  with  setting  1/RS  to  0  when  RS  is  0»  allows  SPICE  to  treat 
RS  consistently  during  analysis  regardless  of  the  value  of  tne 
resistance. 
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13.2.19  BJT  MODEL 


ID  ■  22 

-  1:  sudckt  info 

LOC+  0:  next-pointer 

+  1:  LOCV  LOCV+  0 

model  name 

+  2-  model  type  +  1 

BF 

+  2 

BR 

+  3 

IS 

+  4 

RB 

+  5 

RC 

+  6 

RE 

+  7 

VA 

+  8 

VB 

+  9 

IK 

+10 

C2 

+11 

NE 

+  12 

IKR 

+13 

C4 

+  14 

NC 

+  15 

TF 

+16 

TR 

+  17 

CCS 

+  18 

CJE 

+  19 

PE 

+  20 

ME 

+21 

CJC 

+  22 

PC 

+  23 

MC 

+24 

EG 

+  25 

PT 

+  26 

KF 

+  27 

AF 

+  28 

FC 

+  29 

f  (FC.PE) 

+  30 

f (FC,PE) 

+  31 

f (FCfPE) 

+  32 

f (FC,PC) 

+  33 

f (FC,PC) 

+  34 

f (FC,PC) 

+35 

f (FC,PC) 

+  36 

vcrit 

Comments: 

1)  "model  type"  is  +1  for  "npn"  and  -1  for  "pnp". 
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2)  The  preprocessing  performed  in  subroutine  MODCHK  is  as 
lows: 


jggation 

replaced  cx 

LOCV+  4 

1/RB  (0  if  RB  »  0) 

+  5 

1/RC  (0  if  RC  «  0) 

+  6 

1/RE  (0  if  RE  »  0) 

+  7 

1/VA  (0  if  VA  «  0) 

+  8 

1/VB  (0  if  VS  »  0) 

+  9 

1/IK  (0  if  IK  «  0) 

+  12 

1/IKR  (0  if  IKR  »  0) 

+  28 

FC*PE 

+  29 

PE*  (l-d-FC)”' (1-ME)  )/(l-ME) 

+30 

(1-FC)~ (1+ME) 

+  31 

1-FC* (1+ME) 

+  32 

FC*PC 

+  33 

PC*  (l-d-FC)”'  d-MC)  )/d-MC) 

+  34 

d-FC)”'  d+MC) 

+  35 

1-FC* (l+MC) 

+  36 

VT*ln(VT/(sqrt (2) *IS) ) 

fol- 


3)  See  Comment  4  in  Section  13>2>18  regarding  the  processing  of 
device  model  resistors. 
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13.2.20  JFET  MODEL 


ID  = 

23 

- 

l: 

suockt  info 

LOC+ 

O: 

next-pointer 

-f 

l: 

LOCV 

LOCV+  0 

model  name 

+ 

2: 

model  type 

+  1 

VTO 

+  2 

BETA 

+  3 

LAMBDA 

+  4 

RD 

+  5 

RS 

+  6 

CGS 

+  7 

CGD 

+  8 

PB 

+  9 

IS 

+  10 

KF 

+  11 

AF 

+  12 

FC 

+13 

f  (FC) 

+14 

f  (FC) 

+  15 

f  (FC) 

+  16 

vcrit 

Comments : 


1)  "model  type"  is  +1  for  "njf"  and  -1  for  "pjf". 

2)  The  preprocessing  performed  in  the  MODCHK  subroutine  is 
described  in  what  follows; 


location 


LOCV+  4 
+  5 
+  12 
+  13 
+14 
+15 
+16 


1/RD  (0  if  RD  *  0) 

1/RS  (0  if  RS  =  0) 

FC*PB 

PB*  (l-d-FC)'*  (1-M)  )/(l-M) 
(1-FC)“ (1+M) 

1-FC* (1+M) 

VT*ln(VT/(sqrt(2)*IS) ) 


Note:  for  the  JF£T<  M  is  not  a  model  parameter  and  is  fixed  at 

0.5. 


3)  See  Comment  4  in  Section  13.2.18  regarding  the  processing  of 
device  model  resistors. 
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13.2.21  MOSFET  MODEL 
ID  =  24 


-  1:  suDCKt  info 


LOC-t-  0:  next-pointer 

+  1:  LOCV  LOCV+  0 

model  name 

+  2:  model  type  +  1 

VTO 

+  2 

KP 

+  3 

GAiMMA 

+  4 

PHI 

+  5 

LAi'lBDA 

+  6 

RD 

+  7 

RS 

+  8 

CGS 

+  9 

CGD 

+  10 

CGB 

+11 

CBD 

+  12 

CBS 

+13 

TOX 

+  14 

PB 

+  15 

JS 

+  16 

NSUB 

+  17 

NSS 

+  18 

NFS 

+  19 

XJ 

+20 

LD 

+  21 

NGATE 

+  22 

TPS 

+23 

UO 

+  24 

UCRIT 

+  25 

UEXP 

+  26 

UTRA 

+27 

KF 

+28 

AF 

+29 

FC 

+  30 

f  (FC) 

+  31 

f  (FC) 

+  32 

f  (FC) 

+  33 

VINIT 

+34 

vbi 

+  35: 

xd 

Comments : 

1)  "model  type"  is  +1  for  "nmos"  and  -1  for  "pmos". 

2)  "VINIT"  is  the  initial  guess  used  for  the  MOSFET  diode 
tages. 
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2)  The  preprocessing  performed  oy  subroutine  MODCHK  is 
described  oelow: 


replaced 

LOCV+  2 

KP*COX 

+  3 

sqrt ( 2*EPSSIL*CHARGE*NSUB )/COX 

+  4 

2*VT*ln(NSU3/NI ) 

+  6 

1/RD  (0  if  RD  «  0) 

+  7 

1/RS  (0  if  RS  «  0) 

+  13 

EPSOX/TOX  (COX;  0  if  TOX  «  0) 

+29 

FC*PB 

+  30 

PB^d-Cl-FO'd-M)  )/d-M) 

+31 

d-FC)“  (l+M) 

+  32 

1-FC*  d+M) 

+  33 

VT*ln(VT/(sqrt(2)*JS) ) 

+  34 

VTO-GAi'U-lA*  sqrt  ( PHI ) 

+  35 

sqrt(EPSSIL~2/(CHARGE*NSU3)  ) 

Note:  for  the  MOSFETi  M  is  not  a  model  parameter  and  is  fixed  at 
0.5.  Also,  some  of  the  above  pre-processing  is  a  function  of 
which  model  parameters  were  specified  by  the  user.  See  the  list¬ 
ing  of  suoroutine  MODCHK  for  the  exact  details* 

3)  See  Comment  4  of  Section  13.2.18  regarding  the  processing  of 
device  model  resistors. 
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•PRINT  DC 


13.2.22  .PRINT  DC 


ID  =  31 


LOC+  0 
+  1 
+  2 
+  3 
+  4 
+  5 
+  6 
+  7 
+  8 
+  9 
+  10 
+  11 
+  12 
+  13 
+  14 
+  15 
+  16 
+  17 
+  18 
+  19 


subckt  info 

next-pointer 

LOCV 

analysis  type 
nuniber  of  variables 
ptr (variable  1} 
typ (variable  1) 
ptr(variable  2) 
typ(variable  2) 
ptr(variable  3) 
typ (variable  3) 
ptr (variable  4) 
typ ( variaole  4) 
ptr(variaDle  5) 
typ(variable  5) 
ptr(variaDle  6) 
typ(variable  6) 
ptr (variable  7) 
typ(variable  7) 
ptr(variable  8) 
typ (variable  8) 


LOCV+  0:  pseudo-name 


Comments : 

1)  "pseudo-name"  is  meaningless  (simply  a  unique  "name"  so  that 
subroutine  FIND  may  be  used). 

2)  “analysis  type"  is  actually  redundant  information; 


typ^ 


1  dc  transfer  curve 

2  transient 

3  ac 

4  noise 

5  distortion 

3)  The  number  of  output  variables  must  be  oetween  1  and  3. 

4)  ptr(var)  means  a  pointer  to  the  appropriate  output  descrip¬ 
tion  on  the  appropriate  output  variable  list. 


5)  typ(var)  describes  the  form  in  which  the  variable  is  to  be 
output.  Assuming  that  the  variable  is  "V(l)"  then  the  following 
table  gives  the  interpretation  of  “typ": 
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PRINT  DC 


tVP(var) 


actual  iPtgyprgtatiQn  sii  "van" 


1 

2 

3 

4 

5 

6 


V(l) 

VM  ( 1 ) 
VR(1) 
VI  (1) 
VP(1) 
VDB(1 ) 


6)  The  linked-list  entries  for  ID  »  32#  33»  34i  and  35  are  ex¬ 
actly  the  same  as  for  ID  >  31  except  that  the  analysis  types  are 
transient,  ac,  noise,  and  distortion,  respectively. 
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PLOT  DC 


13.2.23  .PLOT  DC 
ID  =  36 

-  1:  suDckt  info 
LOC+  0:  next-pointer 
+  l:  LOCV 

+  2:  analysis  type 
+  3:  number  of  variables 
+  4s  ptr(variable  1) 

+  5:  typ(variable  1) 

+  6:  ptr(variable  2) 

+  7:  typ (variable  2) 

+  8:  ptr(variable  3) 

+  9:  typ{variable  3) 

+10:  ptr(variable  4) 

+11:  typ{variable  4) 

+12:  ptr(variable  5) 

+13;  typ(variable  5) 

+14;  ptr(variable  6) 

+15:  typ(variable  6) 

+16:  ptr(variable  7) 

+17:  typ(variable  7) 

+18:  ptr(variable  8) 

+19:  typ(variable  3) 

Comments : 


LOCV+  0:  pseudo-name 

+  1:  plo(variable  1) 
+  2:  phi(variaole  I) 
+  3:  plo(variable  2) 
+  4:  phi(variaole  2) 
+  5s  plo(variable  3) 
+  6:  phi(variaDle  3) 
+  7:  plo(variaole  4) 
+  8:  phi (variable  4) 
+  9:  plo(variable  5) 
+10:  phi(variaDle  5) 
+11:  plo(variable  6) 
+12:  phi(variable  6) 
+13:  plo(variable  7) 
+14;  phi (variable  7) 
+15:  plo(variable  8) 
+16:  phi(variable  3) 


1)  "plo(var)"  and  “phi(var)“  are  the  lower  and  upper  plot  limits 
for  the  indicated  output  variable.  A  value  of  0.0  means  that  no 
plot  limit  was  specified. 


2)  The  linked-list  entries  for  ID  =  37  ,  33  ,  2'J ,  and  40  are  ex¬ 
actly  the  same  as  for  ID  =  36  except  that  the  analysis  types  are 
transient,  ac,  noise,  and  distortion,  respectively. 
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13.2.24  DC  ANALYSIS  OUTPUT  VARIABLE 


ID  *  41 


-  1 
LOC+  0 
+  1 
+  2 
+  3 
+  4 
+  5 
+  6 


subckt  info 
next-pointer 


LOCV 

nl  1 

ptr (source) 

n2  I 

— 

ISEU  1 

type  0  1 

1 

<unused> 

LOCV-i-  0:  variable  name 

I  — 

1  — 

I 

12-8 


Conunents : 


1)  “variable  name"  is  an  artificially-constructed  unique  iden¬ 
tifier  for  the  output  variable.  It  has  one  of  3  forms: 


output  variable 


"variable  name 


v (nl »  n2 ) 
i ( vxxxxx) 
noise/distortion 


(nl*2“12)*2“18  +  (n2*2''i2)*2~6 
vxxxxx  (the  name  of  the  source) 
desired  output  -  for  example#  "hd3" 


2)  "ISEO"  is  the  offset  within  the  memory  block  for  the  outputs 
of  a  given  sweep  point  at  which  the  value  for  this  output  vari¬ 
able  is  stored. 


3)  "type"  is  the  type  of  output  variable  (and  determines  the 
meaning  of  (LOC  2)  and  (LOC  +  3)  as  detailed  above): 


type 


meaning 


0 

1 

2 

3 

4 

5 

6 

7 

8 


voltage: 
current: 
noise : 
noise: 
distortion : 
distortion: 
distortion : 
distortion : 
distortion : 


v(nl » n2 ) 

i ( source ) 

onoise 

inoise 

HD2 

HD3 

DIM2 

SI  M2 

DIM3 


4)  The  linked-list  entries  for  ID  *  42,  43,  44,  and  45  are  ex¬ 
actly  the  same  as  for  ID  »  41  except  that  the  output  variables 
are  for  transient,  ac,  noise,  and  distortion  analyses,  respec¬ 
tively. 
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13*3  Table  Specifications 

13.3.1  lELMNT 

All  linked'list  elements  are  stored  in  this  table.  It  is 
always  the  first  table  allocated;  its  origin  is  never  changed  oy , 
the  memory  manager.  SPICE2  takes  advantage  of  that  fact  and  uses 
absolute  (rather  than  relative)  subscripting  for  the  linked  lists 
(that  is»  NODPLC ( variable )  rather  than  NODPLC (lELMNT+of f set)  ) . 

13.3.2  ISBCKT 

This  table  is  used  as  a  stack  of  pointers  to  suocircuit  de¬ 
finitions.  Each  entry  is  one  word,  and  consists  of  the  NODPLC- 
subscript  of  the  subcircuit  definition  currently  being  read  (list 
I0>20).  The  taole  is  required  since  subcircuit  definitions  may 
be  nested. 

13.3.3  lUNSAT 

This  taole  contains  all  names  which  cannot  be  resolved  until 
after  subcircuit  expansion  has  been  performed  (in  the  ERRCHK 
overlay).  Each  entry  is  one  word. 

13.3.4  ITEMPS 

This  table  contains  the  temperatures  at  which  analysis  has 
been  requested  to  be  performed.  The  first  entry  is  always  TNOM; 
if  the  table  contains  more  than  one  entry,  the  first  analysis 
temperature  used  is  the  second  entry  in  the  taole. 
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13.3.5  IFOUR 

This  taole  contains  pointers  to  the  output  variables  (list 
ID=42)  for  which  a  Fourier  analysis  has  been  requested.  , 

13.3.6  ISENS 

This  table  contains  pointers  to  the  output  variables  (list 
ID  =  41)  for  which  a  dc  sensitivity  analysis  has  been  requested. 

13.3.7  IFIELD,  ICODE,  IDELIM,  and  ICOLUM 

These  tables  are  filled  by  subroutine  CARD.  For  each  field 
of  information  scanned  from  the  current  logical  input  line,  one 
entry  is  made  in  each  of  the  tables,  as  follows: 


IFIELD: 


ICODE: 


IDELIM: 


ICOLUM: 


VALUE (IFIELD+n)  contains  either: 

1)  a  name,  stored  in  8H  format, 

2)  a  (floating-point)  number 

NODPLC (ICODE+n )  contains 

+1  =>  corresponding  IFIELD  is  name 
0  *>  corresponding  IFIELD  is  number 
-1  =>  previous  field  was  last  one  of  line 

VALUE ( IDELIM+n )  contains  the  character 
which  delimited  the  field,  stored  in 
IH  format 

NODPLC ( iCOLUM+n )  contains  the  column 
position  in  the  current  line  at  which 
the  field  started 


13.3.8  JUNODE 

This  table  contains  the  input  user  node  numbers.  SPICE2 
uses  a  compact  set  of  element  node  numbers,  numoered  sequentially 
from  1  (for  the  ground  node).  In  routine  ERRCHK,  an  ordered  list 
of  user  node  numbers  is  constructed  in  table  JUNODE,  and  all  ele¬ 
ment  nodes  are  replaced  by  offsets  into  the  JUNODE  taole.  The 
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overhead  is  put  on  obtaining  the  user  node  numbers  since  they  are 
printed  perhaps  three  times  during  a  simulation  runi  while  the 
compact  element  node  numoers  are  needed  thousands  of  times* 

13.3.9  LSBKPT 

This  table  contains  the  transient  analysis  orea]cpoints c  in 
order  of  increasing  time.  Each  entry  is  one  word  (see  Section 
6.7).  The  breakpoint  values  are  obtained  from  the  independent 
source  waveforms  and  time-delayed  elements  (transmission  lines). 

13.3.10  lORDER,  lUR,  ITABLE,  and  ITABID 

These  tables  are  constructed  by  TOPCHKi  and  destroyed  after 
topological  checking  has  been  performed.  The  lORDER  table  is 
used  to  check  that  a  conductive  path  to  ground  from  every  node  in 
the  circuit  exists.  The  other  tables  are  used  in  the  generation 
of  the  node  table,  with  the  following  descriptions: 

lUR:  NODPLC (lUR+n )  is  the  offset  from  the  origin  of 

the  ITABLE  and  ITABID  tables  to  where  the 
pointers  to  elements  connected  to  node  "n** 
are  stored.  The  last  pointer  for  node  "n" 
is  stored  at  the  location  in  ITABLE  indicated 
by  NODPLC (IUR+n+l)-l. 

ITABLE:  NODPLC ( ITABLE+m )  contains  a  pointer  into  the 

lELi-lI'JT  table  to  a  circuit  element. 

ITABID:  NODPLC ( ITASID+m )  contains  the  element  ID  for 

the  element  pointed  to  oy  the  corresponding 
entry  in  the  ITABLE  taole. 
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13.3.11  ISR 

This  table,  together  with  part  of  the  NUMOFF  table,  is  used 
to  record  the  nonzero  circuit  equation  coefficient  matrix  ele¬ 
ments.  The  ISR  table  contains  WSTOP+1  entries.  The  i'th  entry 
is  a  list  header  which  points  to  a  linked-list  stored  in  NUMOFF 
containing  the  numbers  of  the  nonzero  columns  of  the  i'th  row  of 
the  matrix.  The  actual  list  manipulation  is  performed  in 
subroutine  RESERV. 

13.3.12  ISEQ 

This  table  is  used  to  keep  track  of  the  "voltage-defined" 
elements  (that  is,  those  elements  whose  currents  are  unxnowns  in 
the  Modified  Nodal  equations  -  (controlled)  voltage  sources,  in¬ 
ductors,  .  and  the  dependent  sources  within  transmission-line  ele¬ 
ments).  Each  entry  points  to  the  (LOC  +  0)'th  position  of  one  of 
those  elements. 

13.3.13  ISEOl 

This  table  tracks  the  ISEQ  taole.  It  is  used  to  keep  track 
of  where  element  nodes  are  stored,  since  the  nodes  needed  for 
transmission-line  elements  are  stored  in  (relative)  locations 
different  from  those  of  the  other  elements. 

13.3.14  NEQN 

This  table  also  tracks  the  ISEQ  table.  It  records  the 
number  of  the  equation  added  to  the  Modified  Nodal  matrix  to 
evaluate  the  element  current  (see  variable  IBR  in  suoroutine 
t-IATLOC )  . 
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13.3.15  NODEVS 

This  table  records  the  number  of  "voltage-defined"  elements 
connected  to  each  circuit  node.  It  is  used  in  subroutine  REOROR 
to  help  decide  the  order  in  which  to  swap  matrix  rows  (to  elim¬ 
inate  some  singularity  problems). 

13.3.16  NDIAG 

In  the  SETUP  overlay*  this  table  records  whether  the  diago¬ 
nal  elements  of  the  matrix  are  nonzero  (due  to  the  loading  of 
some  element,  or  due  to  propagation  of  fill-in  from  the  LU  fac¬ 
torization  process).  During  dc  analysis,  this  taole  is  used  to 
hold  the  right-hand  side  (of  the  system  of  circuit  equations) 
while  the  reordering  required  oy  the  row-swap  process  is  per¬ 
formed. 

13.3.17  NMOFFC 

The  D ' th  entry  of  this  table  contains  the  number  of  nonzero 
off-diagonal  terms  in  the  j'th  column  of  the  equation  matrix. 

13.3.18  NUWOFF 

The  first  NSTOP  words  of  this  table  contain  the  numoer  of 
nonzero  off-diagonal  terms  in  each  row  of  the  equation  matrix. 
The  rest  of  the  table  is  used  to  hold  the  linked-list  elements 
(pointed  to  by  the  entries  in  the  ISR  table)  which  indicate  which 
columns  are  nonzero  for  a  given  row  of  the  matrix.  Each  list 
element  contains  two  words:  a  pointer  to  the  next  list  element 
(zero  if  none  exists),  and  the  numoer  of  the  column  containing 
the  nonzero  matrix  element.  Note  that  absolute  subscripts  are 
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used  for  the  linked  list  (the  origin  of  tne  NUMOFF  table  must  not 
change  while  matrix  setup  is  in  progress). 

13.3.19  ISWAP 

This  taole  records  the  matrix  row-swap  performed  in 
subroutine  REORDR.  It  is  indexed  using  the  renumoered  row 
numoer*  and  contains  the  corresponding  original  row  number. 

13.3.20  lEQUA 

This  table  is  the  inverse  of  ISWAP:  it  maps  original  matrix 
row  numoers  into  the  set  of  reordered  row  positions. 

13.3.21  lORDER 

This  table  records  the  simultaneous  swapping  of  rows  and 
columns  (equations)  done  to  maximize  the  sparsity  of  the  equation 
matrix.  It  is  indexed  using  the  renumbered  equation  number*  and 
contains  the  corresponding  original  equation  numoer. 

13.3.22  JMNODE 

This  table  is  the  inverse  of  lORDER:  it  maps  equations  into 
their  reordered  position. 

13.3.23  lUR 

This  taole*  in  conjunction  with  the  lUC  taole*  records  the 
position  of  the  upper-triangular  matrix  elements  in  the  one¬ 
dimensional  storage  used  to  hold  the  matrix  coefficients.  The 
i'th  entry  of  the. lUR  table  indicates  the  starting  offset  within 
the  lUC  table  for  the  nonzero  column  numbers  of  the  i'th  row. 
The  last  entry  for  the  i'th  row  is  indicated  oy  the  contents  of 
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the  {i'th  +  1)  entry  minus  1. 

13.3.24  lUC 

This  table*  indexea  throuyh  the  lUR  table*  records  the 
numbers  or  the  nonzero  columns  for  the  upper-triangular  equation 
matrix  terms. 


13.3.25  ILC 

This  table  is  the  dual  of  the  lUR  table*  for  the  lower- 
triangular  matrix  terms.  It  is  accessed  oy  column  numoer*  and 
indicates  the  starting  position  within  the  ILR  table  where  the 
nonzero  row  positions  are  stored. 


13.3.26  ILR 

This  table  is  the  dual  of  the  lUC  taole*  for  the  lower- 
triangular  part  of  the  equation  matrix.  It  contains  the  numbers 
of  the  nonzero  rows  for  given  columns  of  the  matrix- 


13.3.27  MACINS 

This  table  contains  the  machine  code  (generatec  oy 
subroutine  CODGEN )  which  is  used  to  solve  the  system  of  linear- 
izea  circuit  equations. 


13.3-28  LVNIMl 

This  table  contains  the  solution  to  the  linearizea  circuit 
equations  obtained  from  the  previous  Newton-Raphson  iteration- 
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13.3.29  LVN 

This  taiole  contains  the  system  of  circuit  equations  (the 
equation  coefficients  and  the  right-hand  side  vector).  The  vari- 
aoles  LYNL>  LYU>  and  LYL  contain  the  offsets  in  the  LVN  table  to 
the  beginning  of  the  matrix  diagonals  upper-triangle*  and  lower- 
triangle  terms*  respectively- 

13.3.30  LXO 

This  table  contains  nonlinear  device  operating-point  infor¬ 
mation.  During  transient  analysis*  the  taole  also  contains  in-r 
formation  about  capacitor  charges  and  inductor  fluxes.  A  precise 
description  of  the  taole  contents  for  each  circuit  element  is  in¬ 
cluded  in  Section  13-2.  Note  that  this  taole  always  contains  in¬ 
formation  relative  to  the  current  iteration/sweep  point- 

13.3.31  LXl  -  LX7 

These  taoles  contain  the  previous  contents  of  the  LXO  taole- 
In  particular*  LXl  contains  the  previous  contents  of  LXO*  LX2 
contains  the  previous  contents  of  LXl*  and  so  forth. 

13.3.32  LTD 

This  table  stores  previous  values  of  time-delayed  sources 
within  the  circuit.  The  LXi  tables  are  inadequate  for  this  pur¬ 
pose  oecause  more  than  seven  delayed  values  of  some  source  may  be 
necessary. 
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13.3.33  LOUTPT 

This  tsDle  is  used  to  store  tiie  values  of  requested  output 
variables  (node  voltages  and/or  voltage-source  currents)  during 
analysis.  Each  entry  in  the  taole  consists  of  n+1  words;  the 
first  word  is  the  value  of  the  swept  variable^  and  the  rest  are 
the  corresponding  output  variaole  values.  For  ac  analysis.  the 
size  of  each  entry  in  the  table  is  2* (n+1)  words  (for  the  real 
and  imaginary  parts  of  the  complex  output  variaole  values). 

13.3.34  NDIAGC 

This  table  serves  the  same  purpose  for  the  ac  analysis  as 
does  the  NDIAG  table  for  the  dc/transient  analyses  (see  Section 
13.3.6),  except  that  this  table  is  twice  as  long  since  for  the  ac 
analysis  the  right-hand  side  is  complex-valued. 

13.3.35  LVNC 

This  table  is  the  complex-valued  analog  to  the  LVN  table 
(see  Section  13.3.29). 

13.3.36  LDO  and  LDl 

These  taoles  are  used  in  the  small-signal  distortion 
analysis,  which  is  performed  using  suoroutines  DINIT  (for  ini¬ 
tialization)  and  DISTO.  The  LDO  table  contains  suo-computations 
which  need  oe  evaluated  only  once  for  a  given  operating  point. 
The  LDl  table  holds  copies  of  the  solution  vector  at  different 
frequencies  as  the  distortion  analysis  is  performed. 
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13.3.37  LVWCT 

This  tacle  is  only  allocated  if  ooth  noise  and  small-signal 
distortion  analyses  have  been  requested.  It  is  used  to  store  the 
solution  to  the  current  frequency-point  because  both  the  noise 
and  distortion  analyses  destroy  the  right-hand  side  vector. 

13.3.38  LOCX 

This  taole  contains  the  interpolated  x-axis  values  for  out¬ 
put  variables. 

13.3.33  LOCY 

This  taole  contains  the  interpolated  y-axis  values  of  the 
output  variables.  For  n  sweep  points  and  ra  output  variables,  it 
contains  m*n  words,  with  all  the  values  for  the  first  output 
variable  first,  then  the  values  for  the  second  output  variaole, 
and  so  forth. 


-83- 


SPICE2  REPORT 


LABELED-COHi-lON  VARIABLES 


COMi'lOil  VARIABLES 


13-4  Labeled-Common  Variable  Descriptions 


13.4.1  MEMRY 


This  COWi'lON-olock  contains  variables  used  to  manage  the 
dynamically-allocated  tables. 


name 


description 


LORG 

ICORE 

i'UXCOR 

i4AXUSE 

MEliAVL 

LDVAL 

HUHBLK 

LOCTAB 

LTAB 


NODPLC (LORG+n )  is  word  n  in  memory 
total  number  of  memory  words  in  use 
maximum  value  of  ICORE  for  present  circuit 
maximum  number  of  words  required  oy  circuit 
number  of  words  of  available  memory  (in  other 
words,  (words  allocated)  -  (words  required}) 
NODPLC (LDVAL)  is  last  word  of  dynamically- 
managed  memory 

number  of  allocated  taoles  (does  not  include 
the  "table  entry"  table) 

tlODPLC  (LOCTAB+1 )  is  the  first  Word  of  the 
"table  entry"  table 

WODPLC (LTAB+1 )  is  the  first  word  of  the  entry 
in  the  "table  entry"  table  found  oy  function 
MEMPTR 
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13.4.2  TABINF 


This  COMMON 
and  some  table 
table  pointer  to 


name 

lELMNT 

ISBCKT 

NSBCKT 

I  UNSAT 

NUNS AT 

I TEMPS 

NUMTEM 

ISENS 

NSENS 

I  FOUR 

NFOUR 

IFIELD 

ICODE 

IDELIM 

ICOLUM 

INSIZE 

JUNODE 

LSBKPT 

NUHBKP 

I  ORDER 

JMiiODE 

lUR 

lUC 

ILC 

ILR 

NUMOFF 

ISR 

NMOFFC 

ISEJ 

ISEOl 

NEwN 

NODEVS 

NDIAG 

ISWAP 

lEOUA 

i'lACINS 

LVNIMl 

LXO 

LVN 


-block  contains  taole  information  (table  pointers 
lengths).  (The  notation  "tp ( something ) "  means  "a 
a  table  containing  'something'.") 


tp (input  data  linked-list  storage) 

tp ( suDcircui t  definition  stack) 

size  of  ISBCKT  table 

tp ( unsatisf ied  name  list) 

size  of  lUNSAT  table 

tp (analysis  temperatures) 

size  of  ITEMPS  table 

tp(dc  sensitivity  outputs) 

size  of  ISENS  table 

tp(Fourier  analysis  outputs) 

size  of  IFOUR  table 

tp (scanned  input  line  fields) 

tp(type  of  data  in  corresponaing  IFIELD  entry) 

tp(deliraiter  character  for  IFIELD  entry) 

tp (beginning  column  number  of  IFIELD  entry) 

size  of  IFIELD,  ICODE,  IDELIM,  and  ICOLUM  tables 

tp (nodes  used  in  circuit  description) 

tp (breakpoint  values  for  transient  analysis) 

size  of  LSBKPT  table 

tp( record  of  row  renumbering) 

tp (inverse (lORDER  table)) 

tpdUC  indices) 

tp  (nonzero  colurms  in  upper  triangle) 
tp(ILR  indices) 

tp (nonzero  rows  in  lower  triangle) 
tp(nuniDer  of  nonzero  off-diagonal  entries  in  row) 
tp(ptrs  to  linked-list  of  matrix  structure) 
tp (number  of  nonzero  off-diagonal  entries  in  column) 
tp(ptrs  to  voltage-defined  elements  (L,E,V,H,T)) 
tp  (element  node  flag  for  corresponding  ISEiJ  entry) 
tp (equation  nuraoer  for  corresponding  ISEJ  entry) 
tp (number  of  voltage-defined  oranches 
inciaent  to  node) 

tp(in  REORDR:  whether  diagonal  is  nonzero; 

in  analysis:  copy  of  LVN) 
tp( record  of  equation  swaps) 
tp (inverse (ISWAP) ) 
tp (machine  instructions  (code)) 
tp(previous  solution  (copy  of  LVN)) 
tp (current  circuit  state  information) 
tp (right-hand-side  for  set  of  circuit  equations) 
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LYNL 

LYU 

LYL 

LXl  -  l:<7 

LDO 

LDi 

LTD 

LOUTPT 


LVN  offset;  matrix  oiagonai  terms 

LVi)  offset:  matrix  upper-triangle  terras 

LVN  offset:  matrix  lower-triangle  terras 

tp(past  values  of  LXO  table) 
tp ( distortion  analysis  constants) 
tp (distortion  analysis  work  space) 
tp (previous  timepoint  information  for  elements 
which  involve  delay  (transmission  lines)) 

tp( saved  output  values  to  be  printed/plottea ) 
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13.4.3  MISCEL 

This  COi'L'lON-hlock  contains  miscellaneous  variaoles  contain¬ 
ing  the  program  version,  execution  time,  etc. 


QAme 

description 

APROG (  ) 

program  version  and 

release  date 

AT  I  ME 

real-time 

cloclc  (  hh 

:  mm : s  s  ) 

ABATE 

date  (  dd 

mmm 

yy) 

ATITLE(  ) 

circuit  title 

(first 

input  line) 

RSTATS (  ) 

program  statistics: 

(1)  - 

cpu 

time: 

READ IN 

(2)  - 

cpu 

time : 

SETUP 

(3)  - 

cpu 

time : 

dc  transfer  curves 

(4)  - 

no- 

iter: 

dc  transfer  curves 

(5)  - 

cpu 

time : 

dc  operating  point 

(6)  - 

no- 

iter: 

dc  operating  point 

(7)  - 

cpu 

time : 

ac  analysis 

(8)  - 

no. 

iter: 

ac  analysis 

(9)  - 

cpu 

time : 

transient  analysis 

(10)  - 

no- 

iter: 

transient  analysis 

(11)  - 

cpu 

time : 

print/plot  generation 

(20)  - 

NSTOP 

(21)  - 

no- 

matrix 

terms  before  REORDR 

(22)  - 

no- 

matrix 

terms  after  REORDR 

(23)  - 

no. 

of  fillin  matrix  terms 

(26)  - 

“operation 

count"  for  matrix 

(27)  - 

percent  sparsity  of  matrix 

IWIDTH 

input  line  width 

LWIDTri 

output  line  width 

NOPAGE 

nonzero  = 

>  inhibit  page  ejects  for  subtitle 
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13.4.4  LINE 


This  COt'lMON- block  contains  variaoies  pertaining  to 
input  line  to  be  scanned. 


pa me  description 

ACHAR  input  character  (in  IH  format) 

AFIELD ()  next  input  line  (in  A8  format) 

OLDLIN ( )  previous  input  line  (in  A8  format) 

KNTRC  column  numoer  in  current  input  line 

KNTLIM  last  column  to  scan  in  current  inpu 


he  next 


line 
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13.4.5  CIRDAT 


This  COMMON 


name 

LOCATE (  ) 

JELCNT (  ) 
NUNODS 
NCNODS 
NUMNOD 

NSTOP 

NUT 

NLT 

NXTRW 

NDIST 

NTLIN 

IBR 

NU  MVS 


clock  contains  miscellaneous  circuit  data. 


description 

LOCATE (n)  points  to  the  first  element  with  ID  =  n 
(in  table  IELi4i4T) 

JELCNT (n)  is  the  numoer  of  elements  with  ID  *  n 
number  of  user  nodes  (in  nominal  input  circuit) 
number  of  circuit  nodes  (after  subckt  expansion) 
total  circuit  nodes  (including  nodes  added  for 
extrinsic  resistances  in  device  models) 
number  of  circuit  equations 
number  of  nonzero  terms  in  upper  triangle 
number  of  nonzero  terms  in  lower  triangle 
size  of  LXO  table 
size  of  LDO  table 
size  of  LTD  table 

temporary  variable:  next  equation  number  available 
for  element  currents 

number  of  voltage-oef ined  elements 
(inductors,  ( in ) dependent  voltage  sources, 
and  transmission  lines) 
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13.4.6  MOSARG 

This  COi'ii'iOW-blocK  contains  the  mosret  device  parameters  used 
in  subroutines  WOSEUN  and  MOSCAP. 


name 

description 

GAMi-lA 

gamma 

BETA 

k' 

VTO 

zero-bias  threshold  voltage 

PHI 

inversion  potential 

COX 

thin-oxide  capacitance 

VBI 

built-in  voltage 

XNFS 

fast  surface  state  density 

Xi'lSUB 

effective  substrate  doping 

XD 

depletion-layer-width  constant 

XJ 

metallurgical  function  aepth 

XL 

channel  length 

XLAHDA 

lambda 

UTRA 

mobility  transfieid  factor 

UEXP 

mobility  exponent 

VBP 

voltage  at  which  mooility  begins  to  degrade 

VON 

threshold  voltage 

VD3AT 

saturation  voltage 

Gd 

d{id)/d(vgs) 

GdBS 

d (id)/d (vos ) 

GOS 

d(id)/d (vds) 

CDRAIN 

drain  current  (id) 
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13.4.7  STATUS 


This  COidi’lOW-Dlock  contains  variables  which  determine  the 
program  analysis  status. 


description 


OMEGA 
TIME 
DELTA 
DELOLD( ) 
AG(  ) 

VT 

XNI 

EGFET 

MODE 


MOD EDO 


I  CALC 
INITF 


METHOD 


lORD 

MAXORD 

HOHCOH 

ITERHO 

ITEMHO 

NOSOLV 


2  *  PI  *  FREO 
Simulation  time 
current  timestep 
previous  timesteps 

truncation-error  estimation  coefficients 
thermal  voltage  (=K*T/0  or  BOLTS *TEMP/CHARG£ ) 
intrinsic  carrier  concentration 
energy  gap  for  mosfets 
analysis  mode: 

1  -  dc  analysis  (subtype  MODEDC) 

2  -  transient  analysis 

3  -  ac  analysis 

dc  analysis  type: 

1  -  dc  operating  point 

2  -  initial  transient  point 

3  -  dc  transfer  curve 
number  of  output  sweep  points 

analysis  state  within  iteration  for  given  sweep  point 

1  -  converge  with  'off'  devices  floating 

2  -  initialize  junction  voltages 

3  -  converge  with  'off'  devices  held  off 

4  -  store  small-signal  parameters 

5  -  first  timepoint  i«i  transient  analysis 

6  -  first  iteration  of  sweep  point: 

predict  junction  voltages 
numerical  integration  method  flag: 

1  -  trapezoidal 

2  -  Gear 

integration  order 
maximum  integration  order 
number  of  noncon vergent  oranches 
iteration  nuraoer  for  current  sweep  point 
temperature  numoer 

1  =>  “use  initial  conditions"  (UIC)  specified 
for  transient  analysis 
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13.4.8  FLAGS 

This  COMMOH-dIocjc  contains  miscellaneous  output  and  error 
flags. 


IPRr^TA 
IPRUTL 
IP RUTH 
IPRiNTN 
IPRNTO 
LIMTIM 
LIdPTS 
LVLCOD 


LVLTIi'l 


ITLl 

ITL2 

ITL3 

ITL4 

ITL5 

IGOOF 

NOGO 

KEOF 


nonzero  ■>  print  accounting  information 
nonzero  =>  print  circuit  summary  list 
nonzero  *>  print  device  model  summary 
nonzero  *>  print  node  taole 
nonzero  =>  print  options 

seconds  of  cpu  time  to  reserve  for  output 
maximum  number  of  sweep  points 
machine  coae  flag: 

1  -  no  machine  code 

2  -  generate  and  execute  machine  code 

to  solve  the  circuit  equations 
timestep  determination  algorithm  flag: 

1  -  iteration  count 

2  -  truncation-error  estimation 
iteration  limit:  dc  operating  point 
iteration  limit:  per  point  of  ac  transfer  curve 
lower  transient  analysis  limit  (if  LVLTIrl  *  1) 
upper  transient  analysis  limit 

iteration  limit  for  entire  transient  analysis 
local  error  flag 

global  error  flag  (error  if  nonzero) 
nonzero  =>  eno-of-file  on  input 
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13.4.9  KNSTNT 

This  COMMON-jalock  contains  various  execution-time  constants. 


namg  description 


TWOPI 

XLOG2 

X LOG 10 

ROOT2 

RAD 

BOLTZ 

CHARGE 

CTOK 

Gi'lIN 

RELTOL 

ABSTGL 

VNTOL 

TRTOL 

CHGTOL 

EPSO 

EPSSIL 

EPSOX 


2  *  PI 
ln(2) 

In(lO) 
sqrt (2) 

1  radian*  in  degrees  (360/TWOPI) 

Boltzmann's  constant 
electronic  charge 

Kelvin  temperature  equivalent  to  0  Centigrade 

minimum  branch  conductance 

relative  convergence  tolerance 

absolute  current  convergence  tolerance 

absolute  voltage  convergence  tolerance 

truncation-error  estimation  tolerance 

charge  error  tolerance 

permittivity  of  free  space 

permittivity  of  silicon 

permittivity  of  silicon  dioxide 
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13.4.10  DC 


This  CGHHON-olock  contains  variahles  pertinent 


to  the  dc 


analyses . 


name 


description 


TCSTAR 

TCSTOP 

TCINCR 

ICVFLG 

ITCELM 

KS30P 


KINEL 


KIDIN 

KOVAR 

KIDOUT 


start  value  for  dc  transfer  curve 
final  value  for  dc  transfer  curve 
increment  for  dc  transfer  curve 
nunvoer  of  dc  transfer  curve  sweep  points 
(0  »>  no  dc  transfer  curve  requested) 
pointer  to  source  to  be  swept 

nonzero  =>  dc  operating  point  must  be  computed 
(small-signal  transfer  function  and/or  dc  sensi¬ 
tivity  analyses  requested) 

small-signal  transfer  function  analysis  flag: 

0  -  no  small-signal  transfer  function 

>0  -  pointer  to  input  source 

ID  of  element  pointed  to  by  KlilEL 
pointer  to  output  variable  for  small-signal 
transfer  function  analysis 
<unused> 
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13.4.11  AC 


This  COWMON 
analyses. 


name 

FSTART 

FSTOP 

FINCR 


SKW2 

REFPRL 

SPW2 

JACFLG 

IDFREw 


INOISE 

NOSPRT 

NOSOUT 

IJOSIN 

IDIST 


IDPRT 


block  contains  variables  pertinent  to  the  ac 


description 

starting  frequency  for  ac  analysis 
final  frequency  for  ac  analysis 
if  frequency  sweep  is  linear/  FIUCR  is  an 
additive  constant. 

if  frequency  sweep  is  logarithmic  (octave/decade ) / 
FINCR  is  a  multiplicative  constant, 
ac  distortion  analysis:  skw2 

ac  distortion  analysis:  reference  power  level 
ac  distortion  analysis:  spw2 
numoer  of  ac  sweep  points  (0  =>  no 
ac  analysis  requested] 
type  of  frequency  sweep: 

1  -  decade 

2  -  octave 

3  -  linear 

nonzero  =>  perform  noise  analysis 
noise  summary  interval:  every  NOSPRT  sweep 
points/  a  noise  summary  is  printed 

pointer  to  noise  analysis  output  variable 
pointer  to  noise  analysis  input  source 
distortion  analysis  flag: 

0  -  no  distortion  analysis 

>0  -  pointer  to  distortion  load  resistor 

distortion  summary  interval:  every  IDPRT  sweep 
points,  a  distortion  summary  is  printed 
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13.4.12  TRAIJ 

This  COt'L'lON-olock  contains  variables  pertinent  to  tne  tran¬ 
sient  analysis. 


ngrog  description 


TSTEP 
TSTOP 
T START 
DELMAX 
TDi«lAX 

FORFRE 

JTRFLG 


transient  analysis  print  interval 
final  simulation  time  for  transient  analysis 
time  at  which  to  begin  saving  outputs 
maximum  value  for  DELTA  (time step) 
maximum  delay  for  delayed-value  elements 
(transmission  lines) 

Fourier  analysis  frequency 
number  of  transient  output  points 
(0  =>  no  transient  analysis  requested) 
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13.4.13  OUTINF 


This  COMMON-clock  contains  variables  used  in  the  generation 
of  output. 


name. 


d.eggriptipn 


STRIi>lG(  ) 
y VAR ( ) 
XSTART 
XIWCR 
ITAB(  ) 
ITypE(  ) 
ILOCy {  ) 


WPOINT 

NUMOUT 

KNTR 

IJUMDGT 


scratch  line  buffer  used  in  generating  output 

output  values  for  one  line 

starting  x-value  for  output 

increment  to  get  next  x-value  for  output 

pointers  to  output  variaole  definitions 

typ()  values  for  output  variaoies 

logarithmic  plot  flag: 

1  -  not  logarithmic 

2  -  logarithmic 

number  of  output  points 

numoer  of  output  for  current  analysis  mode 
temporary  output  variable  counter 
number  of  significant  digits  to  use  for 
taoular  output  variable  listing 
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This  CGi'IMON 
30a  environ.uent . 


name 

USRIDl 

U3RID2 

.'lAXTAP 

ITAPE 

I'lAXECS 

lECS 

MAXi-lEM 

I  MEM 

MAXLIN 

I LINES 

MAXPCH 

IPUNCH 

i-lAXTIH 

IT  I  ME 

MAXPPU 

IPPU 

I EFT I M 

ISPTIM 

MAXDLR 

ICOST 

XCJEX( ) 


LABELED -COMMON  VARIABLES 


/CJE/ 


alocK  contains  variaales  descrioin'^  the  current 
Its  contents  are  set  by  suaroutine  GETCJE* 


dassr.ip.tion 

<unused> 

<unused> 

<unused> 

<unused> 

<unused> 

<unused> 

<unused> 

maximum  amount  of  memory  availaale  to  SPICE 
<unused> 

<unuseo> 

number  of  output  lines  generated 
<unused> 

<unused> 

cpu  time  limit  (in  milliseconds) 
elapsed  cpu  time  (in  milliseconds) 

<unused> 

<unusea> 

<unused> 

<unused> 

<unused> 

cost  of  300.  excluding  lines  printed  (in  cents) 
<unused> 
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13.4.15  BLA^^K 

This  CQHHON-i5lock  contains  the  array  into  which  all  aynamic 
memory  allocation  is  made. 


name 

VALUE  ( ) 
NODPLC (  ) 
CVALUE( ) 


glooal  dynamically-managed  allocation  array 
VALUE-equiva lence  tor  integer  values 
VALUE-equivalence  for  complex  values 
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13.5  MOSFET  Equations 


The  MOSFET  model  in  SPICE2  is  based  on  three  independent 
variables:  VGS t  VDS »  and  VBS.  Variables  used  in  the  model  equa¬ 
tions  are  defined  at  the  end  of  this  section- 


Channel  modulation: 

Leff  =  LO* (1-Lambda *VDS ) 

yariatighi 

Ueff  =  UO* {Vbp/(Vgst-UTRA*VDS ) )~UEXP 

for  (VGS-VTO-U’TRA*VDS  )  >  Vbp 

(otherwise,  Ueff  *  UO) 

Two-dimensional'  body  effect: 

Gammad  =  Gamma* (1-Gl) 

G1  =  X3*(sqrt(l+2*W/Xj)-1)/L0 

W  =  Xd*sqrt (PHIs-VBS ) 

Saturation  voltage: 


Vdsat  =  VGS-Vbi+0. 5*Gammad‘'2*Vl 


VI  =  1-sqrt (1+4* (VGS-Vfb-VBS )/Gammad~2 ) ) 

P.rgin  quirant  = 


Id 

VON 

F2 

VFACT 

FI 

VSATON 


Subthreshold  region 

Beta* (Fl-2*Gammad*F2/3 ) *exp (VGS-VON )/VFACT 

Vbi+Gammad  *  sqrt ( PHIs-VBS ) 

+k*T/q* (l+q*Nf s/Cox+Gamma/2*l/sqrt (PHIs-VBS ) ) 

(PHIs+VSATON-VBS )' (3/2 )- (PHIs-VBS (3/2 ) 

VON-Vbi-VSATON/2- (2/3 ) *Gammad/VSATON *F2 

(VON-Vbi-VSATON/2 ) *VSATON 

minimum(VDS,  Vdsat  at  VGS  *  VON) 
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Linear  region 

Id  =  Beta*  (Fl-2*Gamiiiad*F2/3  ) 

Beta  =  Kp* (Uef f /UO ) * (LO/Lef f ) 

Kp  =  U0*Cox*W/L0 

FI  =  (VGS-Vbi-VDS/2 )*VDS 

F2  =  (PHIs+VDS-VBS )~ (3/2 )- (PHIS-VBS )~ (3/2 ) 


Saturation  region 

As  in 

Capacitance 

linear,  with  "VDS"  replaced  by  "VDSAT 

V 

II 

• 

CBD 

=  CBD(0)/(1-VBD/PB)~0.5  for  VBD 

=  linear  extrapolation  for  VBD 

<  FC*PB 
>  FC*PB 

CBS 

=  same  equation  as  CBD 

CGB 

=  overlap  capacitance  (from  geometry) 

+  Cl (Cox 

Cl(Cox)=  constant  W*L*Cox  for  VGS  <  VTO-PHIs 

=  linear  ramp  to  0  for  VTO-PHIs  <  VGS  <  VTO 
=  0  for  VGS  >  VTO 


CGS  =  overlap  capacitance  (from  geometry)  +  C2(Cox) 

C2(Cox)=  constant  0  for  VGS  <  VTO-PHIs/2 

=  linear  ramp  to  (2/3)*Cox  for  VTO-PHIs/2  <  VGS  <  VTO 
=  constant  (2/3)*Cox  for  VTO  <  VGS  <  VTO+VDS 
=  (2/3 ) *Cox* (1- (VGS-VTO-VDS )~2/(2* (VGS-VTO)-VDS )“2  ) 
for  VGS  >  VT+VDS 

CGD  =  constant  0  for  VGS  <  VT+VDS 

=  (2/3 ) *Cox* (1- (VGS-VTO)~2/(2* ( VGS-VTO ) -VDS )“2 ) 
for  VGS  >  VT+VDS 


Gamma  =  model  parameter 

if  not  specified  and  NSUB  is  nonzero. 
Gamma  is  computed  using  the  expression 


Gamma  =  Xd*sqrt  (Vx+sqrt  ( i+Vx''2  )  )  /  (L0*VDS  ) 

where 

Vx  =  (VDS-Vdsat)/4 
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Xd 

PHIS 

Vfb 

Vbi 


q 

k 

T 


=  depletion  constant  (  sqrt (2*EPSSIL/q*NSUB ) ) 
=  surface  potential  =  2*PHIf 
=  PHIms  -  Qss/Cox 
=  Vfb  +  PHIS 
=  electronic  charge 
*  Boltzmann's  constant 
S'  temperature  (in  degrees  Kelvin) 
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NOTATION 


I3«6  Element  Load  Templates 


This  section  describes  the  various  element  "stamps"  or  loca¬ 
tions  in  the  circuit  equation  coefficient  matrix  to  which  dif¬ 
ferent  elements  add>  The  notation  used  below  is  descrioed  in  the 
following  figure: 


ni  the  i'th  element  node 

Kelt)  the  equation  added  to  the  I'lodified  Hodal  equations  to 
solve  for  the  current  in  element  "elt" 

A()  the  Modified  Nodal  coefficient  matrix 

B()  the  right-hand  siae  for  the  set  of  circuit  equations 


13.6.1  Resistor 

A(nl,nl)  =  A(nl,ni)  +  1/R 
A(nl»n2)  »  A(nl»n2)  -  1/R 
A{n2.nl)  =  A(n2fnl)  -  1/R 
A(n2»n2)  =  A(n2,n2)  +  1/R 


13.6.2  Capacitor 

A(nl»nl)  =  A(ni»nl)  +  C/h  (h  =  timestep;  infinity  for  dc ) 

A(nl»n2)  =  A(nl/n2)  -  C/h  (note:  conductance  depends  upon 

A(n2*nl)  =  A(n2»nl)  -  C/h  the  integration  method) 

A(n2fn2)  =  A(n2/n2)  +  C/h 

B(nl)  =  3{nl)  +  C/h*V(C)  (V(C)  =  voltage  across  C) 

B(n2)  =  B(n2)  -  C/h*V(C) 


13.6.3  Inductor 

A(nl  ,i{L) )  -  +1 
A  (n2  »  i  (L )  )  *  -1 

A(i(L) ,  nl)  «  +1 
A ( i (L )  /  n2  )  »  -1 

A(i(l).i(l))  =  -L/h  (h  *  timestep;  infinity  for  dc ) 

B(i(L))  =  B(i(L))  -  L/h*I(L)  (I(L)  “  current  through  L) 
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13.6.4  Voltage-controlled  Current  Source 

A(nl,ncl)  *  A(nl»ncl)  +  gm 

A(nl»nc2)  »  A(nl»nc2)  -  gm 

A(n2/ncl)  «  A(n2.ncl)  -  gm 

A(n2fnc2)  =  A(n2/nc2)  +  gm 

(gm  *  d (source)/d (controlling  nodes  {ncl»nc2)) 


13.6.5  Voltage-control leu  Voltage  Source 

A{  nl , i{E) )  =  +1 
A (  n2 I i ( E ) )  »  -1 
A  (i  (E)  ,  nl  )  =  +1 
A  (i  (E  )  ,  n2  )  =  -1 
A ( i ( E ) »  ncl )  =  -Av 
A (i (E ) /  nc2 )  =  +Av 

(Av  =  d (source)/d (controlling  voltage  (ncl,nc2)}) 


13.6.6  Current-controlled  Current  Source 

A (nl , i ( Vc ) )  »  +Ai 
A ( n2 , i ( Vc ) }  ®  -Ai 

(Vc  =  voltage  source  through  which  controlling 
current  flows; 

Ai  =  d ( source ) /d (controlling  current  i(V))) 


13.6.7  Current-controlled  Voltage  Source 

A(  nl ,  i(H)  )  =  +1 
A(  n2,  i  (H)  )  =  -1 
A  ( i  ( H  )  f  nl )  =  +1 
A  (i  (H  )  .  n2  )  =  -1 
A(i(H) ,i(Vc) )  *  R 

(R  =  d ( source )/d ( current  flowing  in  Vc)) 
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13.6.8  Independent  Voltage  Source 


A (  nl , i ( V)  ) 
A (  n2 . i ( V )  ) 
A(i(V),  nl) 
A(i(V),  n2) 

A{i(V),i(V)) 


A(  nl , i (V)  )  +  1 
A{  n2,i(V)  )  -  1 
A(i (V) ,  nl )  +  1 
A  { i  ( V ) ,  n2  )  -  1 
A{i(V) ,i(V) )  +  1 


B{i(v) )  »  V 


13.6.9  Independent  Current  Source 

B(nl)  =  B{nl)  -  I 
B(n2)  =  3(n2)  +  I 


SPICE2  REPORT  _ CONVERSIOi'J  PROBLEMS  CALL-BY-ADORESS 

13.7  Potential  Conversion  Problems 

13.7.1  Call-by-address 

An  implicit  assumption  in  the  design  of  the  memory  manage¬ 
ment  package  in  SPICE2  is  that  it  is  possible  to  change  the  value 
of  some  variable  given  its  "address**.  The  (U.C.  Computer  Center) 
system  function  LOCF  is  used  for  that  purpose:  it  takes  a  single 
arguments  and  returns  the  address  of  the  argument.  Proolenis  ar¬ 
ise  if  the  FORTRAN  compiler  uses  ca 11-by- va lue  instead  of  call- 
by-addressi  since  in  that  event  calls  to  LOCF  will  always  return 
the  address  of  the  local  copy  of  the  passed  parameter. 

13.7.2  Dynamic  Region  Size 

SPICE2  is  written  so  as  to  take  advantage  of  the  possibility 
of  dynamically  changing  its  region  size  (the  number  of  words  of 
memory  allocated  to  the  300  oy  the  operating  system).  This 
dynamic  region  variation  is  achieveu  oy  using  "illegal"  suo- 
scripts  in  the  array  VALUE  to  reference  memory  above  (higher- 
addressed  than)  the  last  word  of  code  for  the  program.  Thus,  oy 
changing  the  number  of  words  allocated  to  the  ^oostep,  the  size 
of  VALUE  can  be  effectively  varied. 

If  dynamic  region  variation  is  either  not  aesirable  or  im¬ 
possible.  one  may  simply  declare 

DIMENSION  VALUE (50000) 

and  live  with  a  fixed-size  maximum  amount  of  available  space, 
with  corresponding  waste  for  small-sizea  circuits. 


-106- 


SPICE2  REPORT 


COiJVERSIOi^  PROBLEMS 


REGION  SIZE 


13.7.3  Characters-per-word 

SPICE2  assumes  d  characters  per  real-valued  variable.  This 
assumption  is  "wired"  into  the  code  in  many  places;  fortunatelyr 
it  is  a  reasonable  assumption  for  most  large  computers. 

13.7.4  Words-per-real-value 

SPICE2  was  originally  written  for  the  CDC  6400  computer 
which  has  60-bit  words.  As  a  results  only  one  word  is  used  both 
for  integer  and  real  values.  However,  since  more  than  48  bits  of 
floating-point  precision  are  necessary  for  accurate  solutions  of 
medium-sized  circuits.  most  non-CDC  computers  require  that 
floating-point  variaoles  be  typed  "DOUBLE  PRECISION"  (or 
"REAL*8").  This  causes  problems  in  SPICE  oecause  the  memory 
manager  allocates  space  in  terms  of  "words"  instead  of  "some 
numoer  of  integers/reals"  (which  would  oe  more  machine  indepen¬ 
dent).  Changes  in  SPICE2  to  solve  this  problem  for  a  specific 
computer  are  extensive  oecause  modifications  roust  be  made 
throughout  the  program. 

13.7.5  Array  suoscripts 

As  mentioned  above,  SPICE2  achieves  its  dynamic  region  size 
variation,  together  with  dynamic  memory  allocation,  by  using 
“illegal"  (in  a  strict  ANSI  FORTRAN  sense)  subscripts.  There  are 
two  ways  to  resolve  this  problem  on  computers  which  check  array 
subscripts.  One  may  locate  the  origin  of  the  /BLANK/  COi-li-lON- 
block  at  the  end  of  all  executed  code  for  SPICE2.  (in  which  case 
the  dynamic  region-size  capaoility  may  oe  used).  Alternatively, 
one  can  simply  declare  VALUE  to  be  or  size  25000  (or  some  other 


-107- 


SPICE2  REPORT 


COiJVERSION  PROBLEMS 


ARRAY  SUBSCRIPTS 


large  number)  and  dispense  with  dynamic  region-size 
(and  put  up  with  the  wasted  memory  in  which  this  method 


variations 
results ) . 
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13.8  Program  Listing 

Persons  who  wish  to  ootain  the  SPICE2  program  should  request 
the  updated  version  from  the  Electronics  Research  Laboratory# 
University  of  California.  Berkeley.  California  94720.  The  La¬ 
boratory  charges  only  a  nominal  handling  charge  for  any  of  the 
programs  that  it  has  availaole. 
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